求助利用10年NCEP逐日资料求10年平均高度场

gooo_s

如题，想用10年NCEP逐日位势高度再分析资料，画出10年的梅雨期平均500hPa高度场图，编写的GS如下：
'reinit'
'sdfopen D:\2\hgt\hgt.2001.nc'
'sdfopen D:\2\hgt\hgt.2002.nc'
'sdfopen D:\2\hgt\hgt.2003.nc'
'sdfopen D:\2\hgt\hgt.2004.nc'
'sdfopen D:\2\hgt\hgt.2005.nc'
'sdfopen D:\2\hgt\hgt.2006.nc'
'sdfopen D:\2\hgt\hgt.2007.nc'
'sdfopen D:\2\hgt\hgt.2008.nc'
'sdfopen D:\2\hgt\hgt.2009.nc'
'sdfopen D:\2\hgt\hgt.2010.nc'
'enable print D:\2\hgt\hgt10.gmf'
'set lon 70 140'
'set lat 15 75'
'set lev 500'
'set t 168'
'have1=ave(hgt.1,t+0,t=178,1)/9.8'
'have2=ave(hgt.2,t+1,t=189,1)/9.8'
'have3=ave(hgt.3,t+4,t=193,1)/9.8'
'have4=ave(hgt.4,t-2,t=197,1)/9.8'
'have5=ave(hgt.5,t+9,t=180,1)/9.8'
'have6=ave(hgt.6,t+4,t=193,1)/9.8'
'have7=ave(hgt.7,t+1,t=209,1)/9.8'
'have8=ave(hgt.8,t-7,t=186,1)/9.8'
'have9=ave(hgt.9,t+10,t=196,1)/9.8'
'have10=ave(hgt.10,t+0,t=199,1)/9.8'
'define have=(have1+have2+have3+have4+have5+have6+have7+have8+have9+have10)/10.0'
'set cint 4'
'set xlint 5'
'set ylint 5'
'set grib off'
'd have'
'set clevs 588'
'set ccols 2'
'set cthick 9'
'd have'
'print'
'disable print'
;
但是运行出来就显示这个 ，求助大神看下哪里错了。。还是这种算法不可以。。。

G-mimzy

我觉得是这一句 'set t 168' ，打开了多个数据文件，但是 t=168 是固定的某一时次，之后你再取到几年后的时间段，但还只是 t+0 ，比如 'have10=ave(hgt.10,t+0,t=199,1)/9.8' ，肯定是不在 hgt.10 这个量的时间范围内，所以报错。 'have10=ave(hgt.10,t+t1,t=t2,1)/9.8' ,t1,t2要取相应的更大的数，不然取不到后面几年的数据吧……（好混乱，试试吧）

gooo_s

求大神。。

残月

在定义新变量的时候怎么着也得有 define 关键字吧..建议用循环做，篇幅少很多。。

易小凯

打开文件的时候用
'sdfopen d:\data\ncep\vwnd.2001.nc d:\data\ncep\vwnd.%y4.nc 3653'这种形式
比较方便
计算have的时候用define语句

mofangbao

可能是你的
ave(hgt.1,t+0,t=178,1)/9.8  这里，使得你的时间设置超出了资料范围，你检查一下

mofangbao

先逐语句过一遍，看在哪一条的时候出现了问题

其微

同求教，我也有这个问题><或者有谁直接写一个针对楼主问题的可以画图的gs文件也行啊！好像不能同时打开多个hgt.****.nc似的><

carol

我也遇到同样的问题，单独的把hgt.2后面的.2去掉就行，加上就不行，向大神求教这是什么原因

太阳雨

同问，求指教