求助！关于grads画气压经度剖面图叠加地形问题

wdd1992

要素资料是necp1*1，地形资料也是从论坛下载的，计算44N比湿高度剖面图
先贴出gs，
'reinit'
'open E:/ice317/fnl_20180317_06_00.ctl'
'open E:/ice317/dx1/orog1.ctl'
'set grads off'
'set grid off'
'set lon 75 95'
'set lat 44'
'set t 1'
'set lev 1000 300'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
*计算比湿
'define es=(6.112*exp((17.67*(TMPprs-273.15))/(TMPprs-29.65)))'
'define qs=(0.62197*es/(lev-0.378*es))'
'define qv=qs*RHprs/100'
*输出地形
'set dfile 2'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set lat 44 '
'set lon 75 95'
'set t 1'
'set z 1 8'
'set gxout shaded'
'set csmooth on'
'set clab off'
'set clevs 0'
'set ccols 0 5 0'
'd orog.2'

*设置 等值线标记的颜色，粗细，大小，输出图形
'set gxout contour'
'set clab on'
'd maskout(qv,-orog.2)'
*图注
'draw ylab hPa'
'draw title 44N QV'
'printim E:/ice317/QV.png white'


grads报错
大家看看哪里的问题啊？

wdd1992

之前也尝试了不用maskout函数，先用shaded画地形，再画出比湿的等值线，但是画出的图形在地形里还有等值线的数值。好头疼啊

river

可否先计算出三维的比湿数据，再来画剖面？
'set lat 44'
就像你切瓜或者其他什么的东西，你想一下剖面是一个三维立体的东西的一个切面，你从最开始都没有这个三维立体的东西，你从哪开始切？

wdd1992

river 发表于 2018-5-11 13:19
可否先计算出三维的比湿数据，再来画剖面？
'set lat 44'
就像你切瓜或者其他什么的东西，你想一下剖面是 ...

感谢回复！学习了。看了另外一个你的帖子，最后画出来了。地形最后用了shade2语句。

284329681

wdd1992 发表于 2018-5-14 21:24
感谢回复！学习了。看了另外一个你的帖子，最后画出来了。地形最后用了shade2语句。

shade2语句是什么啊？能否推荐该贴