关于水汽通量图

★冷霜子★

水汽通量叠加风场，又加了个588线，这里588会不会有点多余呢

wenwen952724299

'reinit'
'sdfopen  f:\data2\air.mon.mean.nc'
'sdfopen  f:\data2\rhum.mon.mean.nc'
'sdfopen  f:\data2\uwnd.mon.mean.nc'
'sdfopen  f:\data2\vwnd.mon.mean.nc'
'sdfopen  f:\data2\hgt.mon.mean.nc'
*'enable print e:\chengxu\1.gmf'
'set parea 0.5 10.5 1.5 8.2'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set t 520'
'set lev 850'
'define es1=6.112*exp(17.67*air/(air+243.51))'
'define qs1=0.62197*es1/(lev-0.378*es1)'
'define q1=qs1*rhum.2/100'
'define qx1=q1*uwnd.3/9.8'
'define qy1=q1*vwnd.4/9.8'
'define uvq1=1000*mag(qx1,qy1)'
'define u1=uwnd.3'
'define v1=vwnd.4'
'set lev 500'
'define h1=hgt.5'

'set t 760'
'set lev 850'
'define es2=6.112*exp(17.67*air/(air+243.51))'
'define qs2=0.62197*es2/(lev-0.378*es2)'
'define q2=qs2*rhum.2/100'
'define qx2=q2*uwnd.3/9.8'
'define qy2=q2*vwnd.4/9.8'
'define uvq2=1000*mag(qx2,qy2)'
'define u2=uwnd.3'
'define v2=vwnd.4'
'set lev 500'
'define h2=hgt.5'

'set t 628'
'set lev 850'
'define es3=6.112*exp(17.67*air/(air+243.51))'
'define qs3=0.62197*es3/(lev-0.378*es3)'
'define q3=qs3*rhum.2/100'
'define qx3=q3*uwnd.3/9.8'
'define qy3=q3*vwnd.4/9.8'
'define uvq3=1000*mag(qx3,qy3)'
'define u3=uwnd.3'
'define v3=vwnd.4'
'set lev 500'
'define h3=hgt.5'

'set t 388'
'set lev 850'
'define es4=6.112*exp(17.67*air/(air+243.51))'
'define qs4=0.62197*es4/(lev-0.378*es4)'
'define q4=qs4*rhum.2/100'
'define qx4=q4*uwnd.3/9.8'
'define qy4=q4*vwnd.4/9.8'
'define uvq4=1000*mag(qx4,qy4)'
'define u4=uwnd.3'
'define v4=vwnd.4'
'set lev 500'
'define h4=hgt.5'

*'define uvq=(uvq2+uvq1)/2'
'define uvq=(uvq4+uvq3)/2'
*'define h=(h1+h2)/20'
'define h=(h3+h4)/20'

'set lon 20 180'
'set lat -20 70'
'set gxout shaded'
'set cmin 0'
'set clevs 1 2 3 4 5 6 7 8 9 10'
*'set ccols 21 22 23 24 25 26 27 28 29 30 31 32 33'
'e:\opengrads\output.gs'
'd uvq'
'cbarn'

'set gxout vector'
'set arrowhead 0.06'
'd skip((u1+u2)/2,2.4,1.8);skip((v1+v2)/2,2.4,1.8)'
*'drawvc skip((u3+u4)/2,2.4,1.8) skip((v3+v4)/2,2.4,1.8) 9.5 0.5 0.5 50'
*'print'
*'disable print'
'set gxout contour'
'set clevs 588'
'set cthick 8'
'set ccolor 1'
*'set clab off'
'd h3'
'printim e:\chengxu\shuiqi.png white'
;

clare

请教lz，你的水汽通量gs方便贴出来吗？个人grads出图的水汽通量散度和micpas有差异

星星

588线是怎么加啊

星星

588线多余不多余，要看你研究的内容是什么了，要说明什么，如果与你研究的内容没什么关系，就不要加了

★冷霜子★

星星 发表于 2014-5-27 11:15
588线多余不多余，要看你研究的内容是什么了，要说明什么，如果与你研究的内容没什么关系，就不要加了

设置只显示588就可以了，我是不知道加了588这个图的名字怎么搞

weiguangjian

什么天气过程需要研究500hPa的水汽状况

nunu18

楼主水汽一般研究的是850hpa，和700hpa，越到高层水汽研究的意义不大，你画副高线是想研究什么？

lightxxia

nunu18 发表于 2014-6-1 21:32
楼主水汽一般研究的是850hpa，和700hpa，越到高层水汽研究的意义不大，你画副高线是想研究什么？

个人觉得850的风场和水汽叠加588hpa副高相对来说是较好的配置吧，一般还是研究较低层的水汽输送

梦之蓝

这图是怎么做出来的啊

sweet33

求楼主贴出gs，目前正在学习，想要在论文中加一个这个，跪求~~~