水汽通量散度的图画出来不大对

shuling216

水汽通量散度的图画出来不大对

独奏古筝到天明

'reinit'
'open e:\nc\fnl_20160420_00_00.ctl'
*'set mpdset cnworld guangdong_q'
'set lon 100 140'
'set lat 20 50'
'set lev 850'
'define_colors'
tt=1
while(tt<=10)
'set t 'tt''
'set grads off'
'set grid off'
'define p=lev'
'define e=exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))*(RHprs/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=UGRDprs*q/9.8'
'define q2=VGRDprs*q/9.8'
'set clevs  0  2 4 6 8  10  12  14  16 18 20  '
'set ccols  0    24 23 22 21 43 44 45 46 47 48 49'
'set gxout shaded'
*'set cmin 5'
'd mag(q1,q2)'
*'d q1;q2'
'cbarn'
'set gxout barb'
'd skip(UGRDprs,2);skip(VGRDprs,2)'

'set gxout contour'
'define p=lev'
'define e=(6.112*exp((17.67*(tmpprs-273.16))/(tmpprs-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhprs/100'
'define qdiv=hdivg(q*ugrdprs/9.8,q*vgrdprs/9.8)'
'd qdiv*1e8'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
' draw title 850hPa vaporflux`b(shaded)`n and qdivergence `b(contour)`n 'recs''
'printim e:\nc\850hPa_vaporflux_'tt'_'recs'UTC.png white'
'c'
tt=tt+1
endwhile

'reinit'

独奏古筝到天明

'reinit'
'open e:\nc\fnl_20160420_00_00.ctl'
*'set mpdset cnworld guangdong_q'
'set lon 100 140'
'set lat 20 50'
'set lev 850'
'define_colors'
tt=1
while(tt<=10)
'set t 'tt''
'set grads off'
'set grid off'
'define p=lev'
'define e=exp(17.67*(TMPprs-273.15)/(TMPprs-29.65))*(RHprs/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=UGRDprs*q/9.8'
'define q2=VGRDprs*q/9.8'
'set clevs  0  2 4 6 8  10  12  14  16 18 20  '
'set ccols  0    24 23 22 21 43 44 45 46 47 48 49'
'set gxout shaded'
*'set cmin 5'
'd mag(q1,q2)'
*'d q1;q2'
'cbarn'
'set gxout barb'
'd skip(UGRDprs,2);skip(VGRDprs,2)'

'set gxout contour'
'define p=lev'
'define e=(6.112*exp((17.67*(tmpprs-273.16))/(tmpprs-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhprs/100'
'define qdiv=hdivg(q*ugrdprs/9.8,q*vgrdprs/9.8)'
'd qdiv*1e8'
'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
' draw title 850hPa vaporflux`b(shaded)`n and qdivergence `b(contour)`n 'recs''
'printim e:\nc\850hPa_vaporflux_'tt'_'recs'UTC.png white'
'c'
tt=tt+1
endwhile

'reinit'