如何画这种某个站的比湿的时间-高度垂直分布图

爱喝柠檬汁

图片是文献中的图，我也想要画出这样的图，下面是我的gs,请各位路过的大神帮我看看呗，哪里的问题，非常感谢！'reinit'
'sdfopen E:/grads/air.2017.nc'
'sdfopen E:/grads/rhum.2017.nc'
'set lon 0 180'
'set lat 0 90'
'set lev 1000 100'
'set t 583 591'
'set grads off'
'set grid off'
'set mpdset hires'
'set map 15 1 1'
'define p=lev'
'define e=exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)*6.112'
'define q=622*e/(p-0.378*e)'
'set lon 82'
'set lat 44'
'set t 583 591'
'set lev 1000 100'
'set cint 1'
'set ccolor 1'
'set cstyle 1'
'set cthick 2'
'set gxout contour'
'd q'
'printim E:/grads/pou.jpg white'
；

Masterpiece

'set lon 531'
你再检查这个吧...

爱喝柠檬汁

Masterpiece 发表于 2017-11-3 18:00
'set lon 531'
你再检查这个吧...

不好意思，那是我编辑时候不小心犯的错误，我已经改过来了，我实际的gs没有这个问题，还是画不出来呢，

爱喝柠檬汁

这个是grads提示的错误，请高手指教

river

画图之前先了解一下你自己用的资料，好好看看你的相对湿度有多少层，上限是哪一层？

Masterpiece

'sdfopen E:/grads/air.2017.nc'
'set dfile 1'
'sdfopen E:/grads/rhum.2017.nc'
'set dfile 2'

'define e=exp(17.67*(air.1-273.15)/(air.1-29.65))*(rhum.2/100)*6.112'

你试试这样改

爱喝柠檬汁

Masterpiece 发表于 2017-11-3 22:21
'sdfopen E:/grads/air.2017.nc'
'set dfile 1'
'sdfopen E:/grads/rhum.2017.nc'

试了，和我之前的一样，应该不是这个毛病，感谢亲

爱喝柠檬汁

river 发表于 2017-11-3 21:15
画图之前先了解一下你自己用的资料，好好看看你的相对湿度有多少层，上限是哪一层？

嘻嘻，我查了一下，只有lev=1000这一层，这样就画不了高度剖面了吧，我想画文献里那个（如图）固定经纬度的比湿的时间—高度图，应该用啥资料呢，咋设置命令呢，求大神帮帮我呀，

爱喝柠檬汁

想画文献里这种图呀！，@river

爱喝柠檬汁

文献里说的是，850hpa某一站的比湿的时间-高度剖面，我不懂850hpa就是一层的量啊，怎么能做出垂直分布呢、请大神赐教，感谢，@river