- 积分
- 2310
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-8-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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行就是程序的最后一行(使用插值函数那一行)
谢谢各位大佬!
|
-
错误提示
-
读取的数据文件
|