【求助】如何用grads做高度场距平的经度—时间剖面

洛喵喵 · 发表于 2018-4-28 16:09:25

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

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

x

想画这种2011 年5 月21 日至6 月30 日20°—25°N 平均的200 hPa 位势高度距平的经度—时间剖面图：

资料是

NCEP的

2.5x2.5的2011年高度场逐日资料hgt2.2011.nc，和1921-2010年30年平均的逐日资料hgt.day.1981-2010.ltm.nc。

下面是我的gs文件：
'reinit'
'sdfopen C:\ziliao\hgt2.2011.nc'
'set dfile 1'
'set lev 200'
'set lon 30 130'
'set lat 20'
'set t 145 181'
'define wp=ave(hgt.1,lat=20,lat=25)/10'
'sdfopen C:\ziliao\hgt.day.1981-2010.ltm.nc'
'set dfile 2'
'set lev 200'
'set lon 30 130'
'set lat 20'
'set t 145 181'
'define wp30yr=ave(hgt.2,lat=20,lat=25)/10'
'define jp=wp-wp30yr'
'd jp'

;

报错：

还有一种方法是两个nc文件中的变量直接减：
'reinit'
'sdfopen C:\ziliao\hgt2.2011.nc'
'sdfopen C:\ziliao\hgt.day.1981-2010.ltm.nc'
'set lev 200'
'set lon 30 130'
'set lat 20 25'
'set t 145 181'
'define jp=hgt.1-hgt.2'
'set lon 30'
'define wp=ave(jp,lat=20,lat=25)/10'
'd wp'

;
报错：

好像因为30年平均的日期是1:1:1:0，而2011年日期是2011:1:1:0，所以不能直接相减，在论坛上搜索方法，说要define转移变量，可是也不太会，尝试多次都报错了。

因为这幅图已经卡了很多天，论坛也没有找到合适解决方式，请求大神指教~小白第一次发帖，如有做的不合适的地方欢迎指出，谢谢！