关于水汽通量图

sxdmxw

很好，有意义。

楚狂

谢谢楼主 下载看看

ITCZ

588线怎么加上去的呀？

★冷霜子★

ITCZ 发表于 2017-5-3 16:05
588线怎么加上去的呀？

设置只显示588等高线啊，高度场里找到588，只显示它。

signor9

感谢楼主，学习了一个

唯妮

学习了。图很漂亮。反映水汽路径很明显

初翠翠

图很漂亮，但是没有gs文件

15174078109

初翠翠 发表于 2017-9-25 18:38
图很漂亮，但是没有gs文件

人家贴出来了,在16楼

juliauh

要看副高是不是必须要说明的， 文章中应该有相应的描述

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'
;