爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6167|回复: 1

关于Cressman插值

[复制链接]
发表于 2019-7-28 09:43:39 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
begin
        f=addfile("C:/Users/Pam/Desktop/precip.cma.daily.195101-201903.nc/precip.cma.daily.195101-201903.nc","r")
        lat=f->stnlat
        lon=f->stnlon
        ;time=f->time(24773:24804) ;20180801-20180831
        ;Time=cd_calendar(time,2)
        ;asciiwrite("C:/Users/Pam/Desktop/precip.cma.daily.195101-201903.nc/2.txt",Time)
        iStation=ind(lat.ge.25.and.lat.le.35.and.lon.ge.95.and.lon.le.110)
        ;asciiwrite("C:/Users/Pam/Desktop/precip.cma.daily.195101-201903.nc/iStation.txt",iStation)
        pre=f->pretotal(24773:24804,iStation)
        StationPre=dim_sum_n(pre,0)
        Lat=f->stnlat(iStation)
        Lon=f->stnlon(iStation)
        oLat=fspan(25,35,21)
        oLat!0="Lat"
        oLat@long_name="Latitude"
        oLat@units="degrees_east"
        oLat&Lat=oLat
        oLon=fspan(95,110,31)
        oLon!0="Lon"
        oLon@long_name="Longitude"
        oLon@units="Longitude"
        oLon&Lon=oLon
        GridPre=new((/dimsizes(oLon),dimsizes(oLat)/),"float")
        rscan = (/.5 , .3 , .15 , .05/)
        StationPre@_FillValue=32766
        GridPre=obj_anal_ic_deprecated(Lon,Lat,StationPre,oLon,oLat,rscan,False)
end
运行后为什么总会出现如图所示的错误提示,第30行就是程序的最后一行(使用插值函数那一行)
谢谢各位大佬!

错误提示

错误提示

读取的数据文件

读取的数据文件
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2019-7-28 10:14:36 | 显示全部楼层
各位我知道啦!我的站点数据为int型,把它改成float后就好啦
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表