【求助】用grads画某一点上某一层比湿随时间变化的图

Nancy1217

各位大神，第一次画这种一维图，想画某一点上某一个层面的比湿随时间变化的图，比如我想画北京上空700hPa比湿随时间变化的图，该怎么画？
这是我写的gs文件：
'reinit'
'open E:\ctl\9m.ctl'
'set grads off'
'set lon 116'
'set lat 40'
'set lev 700'
'set t 3 30'
'define p=lev'
'define es=(6.112*exp(17.67*(tc-273.15)/(tc-29.65)))'
'define e=es*rh/100'
'define q=(0.622*e/(p-0.378*e))'
i=3
while(i<=30)
'set lon 116'
'set lat 40'
'set lev 500'
'set gxout bar'
'set t 'i''
'd q'
i=i+1
endwhile
'printim E:\1\q\1.png white'
;

最后出图什么都显示不出来。怎么破？

男紫汗

不要循环。

Nancy1217

男紫汗 发表于 2017-5-4 08:55
不要循环。

是要set t  3 30 ?

Sue.Chang

你上面已经设置过时间范围是3到30 不需要在循环以及设置经纬度了 你试试改成gxout contour

Nancy1217

Sue.Chang 发表于 2017-5-5 16:26
你上面已经设置过时间范围是3到30 不需要在循环以及设置经纬度了 你试试改成gxout contour

'reinit'
'open E:\ctl\9m.ctl'

'set grads off'
'set lon 116'
'set lat 40'
'set lev 850'
'set t 3 30'

'define p=lev'
'define es=(6.112*exp(17.67*(tc-273.15)/(tc-29.65)))'
'define e=es*rh/100'
'define q=(0.622*e/(p-0.378*e))'
'set gxout contour'
'd q'
'printim E:\1\q\1.png white'
;

这是我新改的，图是出了。。可是不太对。。

Nancy1217

@Sue.Chang  

Nancy1217

出来了，谢谢各位。
'reinit'
'open E:\ctl\9m.ctl'

'set grads off'
'set lon 116'
'set lat 40'
'set lev 850'
'set t 3 30'

'define p=lev'
'define es=6.112*exp((17.67*tc)/(tc+243.5))'
'define e=es*rh/100'
'define q=(0.622*e/(p-0.378*e))'
'set gxout contour'
'd q'
'printim E:\1\q\1.png white'
;
这是最终的gs。不完美，我会继续努力的！

yqy1204

男紫汗 发表于 2017-5-4 08:55
不要循环。

我想画一个区域平均的时间序列，也就是要用到tloop（aave(pp,lon=73,lon=107lat=30,lat=50)）,这是我的ctl文件
dset E:\nnnnn\st11a.grd
dtype station
stnmap E:\nnnnn\sta.map
undef -32766
title Rain Data Sample
xdef 71 linear 70.0 1.0
ydef 41 linear 15.0 1.0
zdef 1 levels  1000
tdef 50 linear jan1961 1yr
vars 1
pp 0 99 Rain Data
endvars!
我尝试写脚本，却总是出错，问了好多人都不会这个函数。能不能请教一下您

男紫汗

yqy1204 发表于 2017-5-7 20:12
我想画一个区域平均的时间序列，也就是要用到tloop（aave(pp,lon=73,lon=107lat=30,lat=50)）,这是我的ct ...

endvars!这里多个感叹号的原因