求问如何画同期与多年平均高度场的距平

x1n3r

***本人新手***      
情况是这样的，
资料是2.5×2.5°的NC逐日资料和NC的ltm多年平均资料，
要做的是四天的平均高度场与同期多年平均高度场的距平。
以下是我的gs：

'reinit'

'sdfopen hgt.day.1981-2010.ltm.nc'
'sdfopen hgt.2012.nc'
'enable print anomaly-500hPa-height.gmf'

'set parea 1 10 1 8'
'set frame circle'
'set mpdset cnworld'
'set mproj nps'
'set grads off'
'set map 1 1 5'

'set lon -90 270'
'set lat 20 90'

'set lev 500'
'define aveH1=ave(hgt.1,t=335,t=338)'
'define aveH2=ave(hgt.2,t=336,t=339)'
'define anom=aveH2-aveH1'

'set gxout contour'
'set ccolor 4'
'set cthick 5'
'd anom/10'

'q w2xy 90 18'
x1=subwrd(result,3);y1=subwrd(result,6)
'set string 1 l 5 0'
'draw string 'x1' 'y1' 90E'
'q w2xy 180 16'
x2=subwrd(result,3);y2=subwrd(result,6)
'set string 1 c 5 0'
'draw string 'x2' 'y2' 180E'
'q w2xy 270 18'
x3=subwrd(result,3);y3=subwrd(result,6)
'set string 1 r 5 0'
'draw string 'x3' 'y3' 90W'
'q w2xy 360 16'
x4=subwrd(result,3);y4=subwrd(result,6)
'set string 1 c 5 0'
'draw string 'x4' 'y4' 180W'
'print'
'printim anomaly-500hPa-height.png white'
'disable print'
'reinit'

运行了之后总是提醒错误，好像意思是需求超出了文件限制。

请大家帮忙看看吧~{:soso_e154:}  

kongfeng0824

请问楼主2.5×2.5°的NC逐日资料是哪里下载呢？时间跨度是多长呢？

river

具体说下NC的ltm多年平均资料

x1n3r

kongfeng0824 发表于 2013-5-22 01:26

                               
登录/注册后可看大图

请问楼主2.5×2.5°的NC逐日资料是哪里下载呢？时间跨度是多长呢？

就是一年的逐日资料，我用的是2012年的资料，是闰年的，366天的。ltm是365天的多年平均资料。

x1n3r

river 发表于 2013-5-22 08:30

                               
登录/注册后可看大图

具体说下NC的ltm多年平均资料

资料是1981-2010年多年高度场逐日平均资料，分辨率2.5×2.5°的。

river

x1n3r 发表于 2013-5-22 21:30

                               
登录/注册后可看大图

资料是1981-2010年多年高度场逐日平均资料，分辨率2.5×2.5°的。

在set lev 500后面，加上set t 1试试

x1n3r

river 发表于 2013-5-22 22:54

                               
登录/注册后可看大图

在set lev 500后面，加上set t 1试试

我是做四天平均的，如果在set lev 500后面加set t 1 之后的平均不就已经覆盖了嘛？

river

x1n3r 发表于 2013-5-22 23:00

                               
登录/注册后可看大图

我是做四天平均的，如果在set lev 500后面加set t 1 之后的平均不就已经覆盖了嘛？

这个和你的四天平均没什么关系啊，你求完平均之后不就相当于只剩下一个时次了么，这时候的set t 1就是代表这么个意思。其实那个不一定是1，也可以换成其他数字。你先试试能不能出图，然后在讨论这个，否则都是没有意义的

x1n3r

river 发表于 2013-5-23 07:09

                               
登录/注册后可看大图

这个和你的四天平均没什么关系啊，你求完平均之后不就相当于只剩下一个时次了么，这时候的set t 1就是代表 ...

恩呢，谢谢，我试过了呢，还是一样的错误。

x1n3r

river 发表于 2013-5-23 07:09

                               
登录/注册后可看大图

这个和你的四天平均没什么关系啊，你求完平均之后不就相当于只剩下一个时次了么，这时候的set t 1就是代表 ...

'reinit'

'sdfopen hgt.day.1981-2010.ltm.nc'
'sdfopen hgt.2012.nc'
'enable print anomaly-500hPaltm-height.gmf'

'set parea 1 10 1 8'
'set frame circle'
'set mpdset cnworld'
'set mproj nps'
'set grads off'
'set map 1 1 5'

'set lon -90 270'
'set lat 20 90'
'set lev 500'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set clopts 1 1 0.10'

'set gxout contour'
'set cterp on'
'set csmooth on'
'set lon -90 270'
'set lat 20 90'
'set z 1'
'set t 1'
'define aveH1=ave(hgt.1,t=335,t=338)'
'define aveH2=ave(hgt.2,t=336,t=339)'
'define anom=aveH2-aveH1'
'set ccolor 4'
'set cthick 5'
'd anom/10'

'q w2xy 90 18'
x1=subwrd(result,3);y1=subwrd(result,6)
'set string 1 l 5 0'
'draw string 'x1' 'y1' 90E'
'q w2xy 180 16'
x2=subwrd(result,3);y2=subwrd(result,6)
'set string 1 c 5 0'
'draw string 'x2' 'y2' 180E'
'q w2xy 270 18'
x3=subwrd(result,3);y3=subwrd(result,6)
'set string 1 r 5 0'
'draw string 'x3' 'y3' 90W'
'q w2xy 360 16'
x4=subwrd(result,3);y4=subwrd(result,6)
'set string 1 c 5 0'
'draw string 'x4' 'y4' 180W'

'print'

'printim anomaly-500hPa-height.png white'

'disable print'

'reinit'


受你刚才所说的启发，我又稍改了下程序，可是还有错误。