季度距平求教啊

曼联现在没$

278803532 发表于 2014-5-17 20:51
可以先把数字写出来再画图，在GRADS写数据的时候用AAVE函数就可以实现，时间一个循环，层次一个循环，层次 ...

大概懂了你的思路，但是自己做起来还是有困难，用了你的方法之后，显示不出图，是一片黑

278803532

曼联现在没$ 发表于 2014-5-18 10:04
大概懂了你的思路，但是自己做起来还是有困难，用了你的方法之后，显示不出图，是一片黑

你写数据没写错吧，
写数据的文件，类似

it=1   
while(it<=6)   
'set lon 90 120'
'set lat 10 40'
*----------第一次循环---------------------
ih=4   *----------1948年4月对应时间
ik=    *----------2012年4月对应时间

iy1=ih+it *----------1948年5月对应时间
iy2=ik+it *----------2012年5月对应时间
*----------第一次循环---------------------

iz=1
while(iz<=21)
'set z 'iz''

'set lon 90 120'
'set lat 10 40'

'define rh1=rh(t='iy2')-ave(rh,t='iy1',t=iy2'，12）
'define gg=aave(rh1,lon=97,lon=106,lat=21,lat=29)'
'd gg'

iz=iz+1
endwhile
it=it+1
endwhile

**************
描述文件

dset        D:\zdz\rh.grd
title  gyleftht
undef 9.999e+20
XDEF  1 LINEAR     100   1
YDEF  1 LINEAR     21  1
zdef 21 levels  1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200

150 100
tdef 6 linear may2012 1mo
vars   1
n  21 99 temp
endvars

**************
画图
'set zlog on'
'set ylevs 1000 850 700 600 500 400 300 200 100'

'set gxout contour'
'set grads off'
'set grid off'
'set mpdset mres'
'set csmooth on'
'set tlsupp year'

'set gxout shaded'
'set lev 1000 100'
'set t 1 6'
'set cint 0.000005'
'set cmax 0'
'd n'
'print'

'set gxout contour'
'set cint 0.000005'
'd n'
'print'

只能帮你到这里了，
再不会我也没有办法了，都不知道你哪里错了。。。

曼联现在没$

278803532 发表于 2014-5-18 10:33
你写数据没写错吧，
写数据的文件，类似

好的，我再试试，十分感谢你

曼联现在没$

278803532 发表于 2014-5-18 10:33
你写数据没写错吧，
写数据的文件，类似

d gg不是显示gg变量吗，这没看懂，为什么循环一次要显示一次gg

278803532

曼联现在没$ 发表于 2014-5-18 11:32
d gg不是显示gg变量吗，这没看懂，为什么循环一次要显示一次gg

，程序都看不懂木有办法了，第一次循环是写5月的，第二次循环是写6月的，7,8,9,10一直这样写出来，你不写出来时间轴怎么会有变量。。。。

书道阁

不错哈！。。。

曼联现在没$

278803532 发表于 2014-5-18 15:35
，程序都看不懂木有办法了，第一次循环是写5月的，第二次循环是写6月的，7,8,9,10一直这样写出 ...

额，我是初学者请见谅，但是d不是绘图指令吗，每次d不会要显示一张图？那循环里有d不是有好多图？还是d有其他功能

曼联现在没$

书道阁 发表于 2014-5-18 17:33
不错哈！。。。

啥不错。。。。。。。。。。。。。。？

曼联现在没$

排山倒海 发表于 2014-5-12 10:12
用grads读取nc资料，转成grd格式的，再用fortran处理，如果直接用grads好像不太容易

好麻烦啊。。。。。。。

曼联现在没$

月是故乡明 发表于 2014-5-13 11:34
多个层次是可以输出在一个grd中的，只需要在fwrite过程中使用do while循环语句即可
且亦可以直接使用gra ...

我试了，出不来图，一片黑，应该是维数设置的不对