请教一下如何修改-水汽通量垂直剖面图

flyflyfishna

我需要画3天4时次34°N，100-120°E水汽通量垂直剖面图。现在只写了单时次的程序，运行出来X轴是对的，Y轴显示的是数值，我希望是1000-200hpa。程序如下：

'reinit'
'open d:\2011091300.ctl'
'define Es=6.1078*exp(17.27*(tmpprs-273.16)/(tmpprs-35.86))'
'define E=rhprs*Es/100.0'
'define mixr=622*E/(lev-0.378*E)'
*'define moistdiv = hdivg(ugrdprs*mixr,vgrdprs*mixr)*1e4'
'define fluxU=ugrdprs*mixr/9.8'
'define fluxV=vgrdprs*mixr/9.8'
'set gxout vector'
'set csmooth on'
'set grid off'
'set grads off'
'set clab forced'
'set cint 3'
'set lon 100 120'
'set lat 34'
'set lev 1000 200'
'd fluxU;fluxV;mag(fluxU,fluxV)'
'd mag(fluxU,fluxV)'
'enable print d:\20111300.gmf'
'print'
'disable print'
;

天目神眉

set ylabs *|* 就可以Y轴显示你想要的数值了

flyflyfishna

我试着把语句顺序调换了一下出来了，但是不知道数算得对不对。
'reinit'

'open d:\2011091300.ctl'
'set lev 1000 200'
'define Es=6.1078*exp(17.27*(tmpprs-273.16)/(tmpprs-35.86))'
'define E=rhprs*Es/100.0'
'define mixr=622*E/(lev-0.378*E)'
*'define moistdiv = hdivg(ugrdprs*mixr,vgrdprs*mixr)*1e4'

'define fluxU=ugrdprs*mixr/9.8'
'define fluxV=vgrdprs*mixr/9.8'

'set gxout vector'
'set csmooth on'
'set grid off'
'set grads off'
'set clab forced'
'set cint 3'

'set lon 100 120'
'set lat 34'

'set ylab on'
'draw ylab `2p`1/hPa'
'set strsiz 0.25 0.25'
'set string 1 bc'
'd fluxU;fluxV;mag(fluxU,fluxV)'
'd mag(fluxU,fluxV)'

'enable print d:\sq1300.gmf'
'print'
'disable print'
;

飞雪残剑

好东东，收起啦！

水云公子

楼主，小弟不才，这久i就纠结于剖面图，想问一下你的ctl文件里面的nc文件是哪一个？或者说这个ctl是关于什么的啊！？

xiayu1990426

flyflyfishna 发表于 2011-11-12 15:52

                               
登录/注册后可看大图

我试着把语句顺序调换了一下出来了，但是不知道数算得对不对。
'reinit'

你好，我也在画这个图，但是我画出来的数值与你相差好大，我的最大也就2.5，我就不知道是怎么回事了，我没有计算q，因为我有现成的shum，

kanbaoyun

能否把ctl文件也分享一下，感谢了！

sunempire22

剖面要插值,很容易错,这个要学习下

云南小黑哥

我也要画这种图，不过要叠加上地形

转音2

很棒呀~谢谢分享~~~~~~~~~~~~~~