求水汽通量散度，大神们请赐教

右手年华

'reinit'
'sdfopen E:\rhum.2015.nc'
'sdfopen E:\uwnd.2015.nc'
'sdfopen E:\vwnd.2015.nc'
'sdfopen E:\air.2015.nc'

*绘图区域
'set vpage 0 11 0 8.5'
'set parea 1 10.5 3.5 8'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set lon 60 160'
'set lat 0 70'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set ylevs 10N 20N 30N 40N 50N 60N 70N 80N'
'set xlevs 40E 50E 60E 70E 80E 90E 100E 110E 120E 130E 140E 150E 160E'
'set clopts 1 1 0.2'
'set ccolor 1'

'set lev 850'
i=158mpdset
while i<163 do
'set t 'i
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.1/100'
'define qu=q*uwnd.2/g*1000'
'define qv=q*vwnd.3/g*1000'
*'define sqd=mag(qu,qv)'
'define sqdiv=hdivg(qu,qv)*1e6'
i=i+1
endwhile

'define vap1=ave(sqdiv,t=158,t=162)'
'define qu1=ave(qu,t=158,t=162)'
'define qv1=ave(qv,t=158,t=162)'

'set lon 70 140'
'set lat 5 50'

'set gxout shaded'
'setrgb'
'set ccols  43 44 45 46 47 48 49 0'
'set clevs -13 -10 -9 -8 -7 -6 -5'
'd vap1'
'set gxout vector'
'd qu1;qv1'
'enable print E:\850sq\850sq.mean.gmf 1024 768 white'
'print'
'disable print'
'c'

右手年华

过年了，人都没了

beiruoyiqiao

我也在想画水汽通量图，求指教

Akbar

这程序要改正吧

nunu18

咦~~~~楼主的问题呢？

深山老怪

这个程序不能运行啊  

深山老怪

这个程序不能运行啊