求问怎么画时间轴是月份的剖面图~~~~

双喜丸子

有个小错误，ctl里是 tdef 12 linear 00z011988 1mo。
改好了还是读不出。。

river

双喜丸子 发表于 2013-8-21 10:53

                               
登录/注册后可看大图

谢谢好心的river~~~~全是0的问题已经解决了~~~换了1988——1999和2000——20111这两段来做~
现在的gs长这 ...

你可以看一下正确的ctl是怎么写的啊，提示都说的很明白了airp有问题。你除了告诉grads文件里变量名是什么，还得告诉它有几层啊。
airp 0 99


还有问一个问题，你运行以上gs有没有什么错误提示之类的

双喜丸子

river 发表于 2013-8-21 11:05

                               
登录/注册后可看大图

你可以看一下正确的ctl是怎么写的啊，提示都说的很明白了airp有问题。你除了告诉grads文件里变量名是什 ...

刚刚的问题解决啦~就是木有写airp 0 99的缘故，应该细心一点的。。。。
我现在想把17层数据全部写进去，可是为什么画出来的图都是一样的捏。。不管set t多少，d出来的airp都是一个样
'reinit'
'sdfopen e:\nc\air.mon.mean.nc'



'set gxout fwrite'
'set fwrite e:\nc\airpm.grd'

'set x 1 144'
'set y 1 73'
l=1
while(l<=12)
'set t 'l''

h=1
while(h<=17)
'set z 'h''
a=625+l-1
b=769+l-1
c=481+l-1
'define airp=(ave(air,t='a',t='b',12)-ave(air,t='c',t='a',12))'
'd airp'
h=h+1
endwhile

l=l+1
endwhile

'disable fwrite'
;

ctl是这样

dset  airpm.grd
undef -9.99e+33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef 12 linear 00Z01jan1988 1mo
vars 1
airp 0 99
ENDVARS

river

双喜丸子 发表于 2013-8-21 16:37

                               
登录/注册后可看大图

刚刚的问题解决啦~就是木有写airp 0 99的缘故，应该细心一点的。。。。
我现在想把17层数据全部写进去， ...

你到底明不明白，变量后面为什么要加0 99 而0和99又是什么意思？grads最基础的东西，能不能自己先看书弄一下呢。既然你都是17层了，还写0，那不就是告诉grads只有地面一层啊，能有什么变化呢

双喜丸子

river 发表于 2013-8-21 19:31

                               
登录/注册后可看大图

你到底明不明白，变量后面为什么要加0 99 而0和99又是什么意思？grads最基础的东西，能不能自己先看书弄一 ...

哦，谢谢river

river

双喜丸子 发表于 2013-8-21 19:35

                               
登录/注册后可看大图

哦，谢谢river

搞懂了就好，下次就会了

双喜丸子

river 发表于 2013-8-21 21:50

                               
登录/注册后可看大图

搞懂了就好，下次就会了

恩恩，还要多学习