相关系数画不出来

王兴菊

请各位高手看一下，是不是ctl写的有问题，画不出图来，水平有限，谢谢了。

夏天爱裤衩

直接把ctl考上来呗，帮你解答问题还要花钱？这个划不来啊

夏天爱裤衩

好吧我，我花的更多，遇到穷神了

artmunich

求相关系数的gs只写lev=1000时的数据，ctl怎么会是17层呢？

另外，如楼上所说，直接贴code

王兴菊

很抱歉啊，没有考虑到会浪费各位前辈的积分，谢谢了。那ctl是不是应该改成这样啊
dset D:\ncep\b.dat   
title NMC reanalysis atlas
undef -9.96921e+36
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 levels 1000  
tdef 760 linear 00Z01JAN1948 1mo
vars 1
res 17 -999 Monthly Mean of res
endvars

这个是在网上下载的gs

'sdfopen  D:\ncep\air.mon.mean.nc'
'sdfopen  D:\ncep\rhum.mon.mean.nc'

'set gxout fwrite'
'set fwrite D:\ncep\b.dat'
'define air=0'
'define rhum=0'

lats=0.0

'set lev 1000'

while(lats<=60.0)

lons=120.0

while(lons<=260.0)

'define res=scorr(ave(air.1,t=1,t=396),ave(rhum.2,t=1,t=396),lon='lons',lon='lons+2.5',lat='lats',lat='lats+2.5')'
'd res'
'print lon='lons',lon='lons+2.5',lat='lats',lat='lats+2.5''

lons=lons+2.5

endwhile

lats=lats+2.5


endwhile

'disable fwrite'
'printim D:\ncep\res.gif gif white x1200 y800'
*run D:\ncep\a.gs

这个是我要出图的gs

'reinit'
'open D:\ncep\a.ctl'
'set grads on'
'set grid on'
*'set mpdset cnworld'
'set map 1 1 10'
'set mpdset guizhou guizhou_q'
'define res=0'
'd res'
'printim D:\ncep\res.gif gif white x1200 y800'
*run D:\ncep\a1.gs

请指教

mofangbao

王兴菊 发表于 2012-7-18 20:49

                               
登录/注册后可看大图

很抱歉啊，没有考虑到会浪费各位前辈的积分，谢谢了。那ctl是不是应该改成这样啊
dset D:\ncep\b.dat   
...

你自己先把图画出来吧，画不出来就先定位一下错误，一下子贴这么多怕是没时间看完

王兴菊

我调试的时候都没有出错了，就是不出图

王兴菊

因为以前没有写过ctl ,所以不知道是不是ctl错了，我自己就是找不到原因啊

x233386794

谢谢分享！！