850hPa单层水汽通量散度

小贝壳

'reinit'
'sdfopen C:\。。。\201605_u.nc'
'sdfopen C:\。。。\201605_v.nc'
'sdfopen C:\。。。\201605_q.nc'
'set mpdset cnworld'
'set ylint 5'
'set xlint 5'
'set map 15'
'set grads off'
'set grid off'
'set lev 850'
i=137
while(i<158)
'set t 'i
'set lat 20 35'
'set lon 100 125'
'set gxout shaded'
*light blue to dark blue
'set rgb 41 225 255 255'
'set rgb 42 180 240 250'
'set rgb 43 150 210 250'
'set rgb 44 120 185 250'
'set rgb 45  80 165 245'
'set rgb 46  60 150 245'
'set rgb 47  40 130 240'
'set rgb 48  30 110 235'
'set rgb 49  20 100 210'
'set clevs -20 -15 -10 -5 0'
'set ccols 49 47 45 43 42 0'
'define qdiv=hdivg(q.3*u.1/9.8,q.3*v.2/9.8)'
'define sqdiv=qdiv*100000000'
'd sqdiv'
'cbarn'
'set gxout contour'
'set clevs  -20 -15 -10 -5 '
'set cthick 6'
'set ccolor 1'
'set cstyle 1'
'd sqdiv'
'set gxout barb'
'set  strmden 1'
'set digsiz 0.07'
'set cthick 6'
'set ccolor 1'
'd skip(u.1,9,9);v.2'
'printim C:\。。。\'i'.png white'
i=i+1
endwhile
;

小贝壳

单位是10-8e /(g*cm2*hPa*s)