关于500hpa月平均高度距平场的问题

随风一路

想做1949至2015年4月的500hpa月平均高度距平图，可是画出来的不是距平！问下，我的gs哪里出问题了？大神们可以帮我改改吗？附上我的gs文件和绘出来的图，谢谢！
gs
'reinit'
'sdfopen D:\wangpei\hgt.mon.mean.nc'
'set fwrite D:\wangpei\ave4.grd'
'set gxout fwrite'
'define ahgt=ave(hgt,t=1218,t=1938,12)'
tt=1218
while(tt<=1938)
'set lev 500'
'set lon 105 225'
'set lat -66.5 66.5'
'set t  ' tt
' nhgt=hgt-ahgt'
' d nhgt '
tt=tt+12
endwhile
'disable fwrite'
'set mpdset cnworld'
'set gxout contour'
'set cint 4'
'set t 8'
'd nhgt'
'enable print D:\wangpei\fugaotai.gmf '
'print'
'disable print'
*'reinit'*'reinit'

午夜狂魂

感谢楼主分享

午夜狂魂

感谢楼主分享

oucdaqi

你在定义ahgt之前就应该先对经纬度时间维和z维进行设置，而且去除气候趋势还有个'modify *** seasonal'的一个命令可以不用像LZ那么麻烦

随风一路

oucdaqi 发表于 2016-4-25 19:10
你在定义ahgt之前就应该先对经纬度时间维和z维进行设置，而且去除气候趋势还有个'modify *** seasonal'的一 ...

可以说说具体怎么改吗？这些设置的顺序应该不会影响吧！

随风一路

oucdaqi 发表于 2016-4-25 19:10
你在定义ahgt之前就应该先对经纬度时间维和z维进行设置，而且去除气候趋势还有个'modify *** seasonal'的一 ...

可以说说具体怎么改吗？这些设置的顺序应该不会影响吧！

oucdaqi

你把ctl粘出来，要不然没法具体。我认为可能有影响，毕竟没试过，lz你可以试试嘛

随风一路

oucdaqi 发表于 2016-4-27 00:10
你把ctl粘出来，要不然没法具体。我认为可能有影响，毕竟没试过，lz你可以试试嘛

试过了，不行呀！这个我直接用gs画的图，然后有关文件的信息直接都是使用sdfopen命令打开获取的，时间开始为1850年。

oucdaqi

随风一路 发表于 2016-4-27 22:54
试过了，不行呀！这个我直接用gs画的图，然后有关文件的信息直接都是使用sdfopen命令打开获取的，时间开 ...

这样吧，LZ，我贴一个我以前用的差不多类似的你瞅瞅，有用最好没用LZ再想办法
'reinit'
'set grads off'
'sdfopen f:\finshfortran\seaice\icec.mnmean.nc'
'set x 1 360'
'set y 1 180'
'set t 1 12'
'define clicec=tloop(ave(icec,t+0,t=412,12))'
'modify clicec seasonal'
'set t 1 412'
'define iceca=icec-clicec'

随风一路

oucdaqi 发表于 2016-4-28 21:48
这样吧，LZ，我贴一个我以前用的差不多类似的你瞅瞅，有用最好没用LZ再想办法
'reinit'
'set grads off ...

先谢谢啦！