斜剖面不能画计算量大点的物理量吗？

酒国

根据river   http://bbs.06climate.com/forum.php?mod=viewthread&tid=10966中的程序，画了风速，可以出图：

而要画水汽通量的斜剖面，没有画出来，是因为水汽通量计算步骤多吗？下面是相关程序，都是根据上面链接中做的。
'reinit'
'open e:/wrf/yuand01/yuan01.ctl'
'set grads off'
'set grid off'
'set zlog on'
'set x 1'
'set y 1'
'set lev 1000 100'
'set t 15'
'set parea 6 10.4 0.5 7'
'define pre=lev'
'define e=(6.112*exp((17.67*(tk-273.16))/(tk-29.65)))'
'define qs=(0.622*e/(pre-0.378*e))'
'define q=qs*rh/100'
'define qu=(u*q/9.8)'
'define qv=(v*q/9.8)'
'define uvq=1000*mag(qu,qv)'
lon1 = 105
lon2 = 113
lat1 = 17
lat2 = 25
lon = lon1
'collect 1 free'
while (lon <= lon2)
  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
  'collect 1 gr2stn(uvq,'lon','lat')'
  lon = lon + 1
endwhile
'set csmooth on'
'set lon 105 113'
'set lev 1000 100'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'
'set xlabs 17N,105E|21N,109E|25N,113E'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set gxout contour'
'set csmooth on'
'set cthick 6'
'set ccolor 1'
'set clopts 0 6 0.15'
'd coll2gr(1,-u)'
'close 1'

男紫汗

斜面图真是溜。。。帮顶

kevin520kaka

看不懂啊，这些太高深了

xuwuwanyou

看不懂，都是大牛，来学习的

river

你自己都知道需要先计算出相应的物理量才能画斜剖面，那你不想想出错的原因在哪？是不是因为所需的物理量没有计算出来？找到原因，再解决不是更容易吗

酒国

river 发表于 2016-9-14 20:37
你自己都知道需要先计算出相应的物理量才能画斜剖面，那你不想想出错的原因在哪？说不是因为所需的物理量没 ...

嗯嗯，好的

deemo7

楼主 最后解决了吗？我也是能画出风速 画不出剖面

酒国

deemo7 发表于 2016-9-23 07:14
楼主 最后解决了吗？我也是能画出风速 画不出剖面

还没。。不过大体懂了。。我画过类似的：用gs提出去水汽通量散度放在grd文件，然后其ctl，然后再把ctl放在画斜剖面的gs里面就可以了

deemo7

酒国 发表于 2016-9-23 16:02
还没。。不过大体懂了。。我画过类似的：用gs提出去水汽通量散度放在grd文件，然后其ctl，然后再把ctl放 ...

好的 谢谢{:5_257:}

micolian

谢谢楼主分享，先果断拷贝，然后自己画一下