求大佬帮我看看我求整层水汽通量的gs

LNYZ

这是我的gs
'reinit'
'sdfopen E:\shum.mon.mean.nc'
'sdfopen E:\uwnd.mon.mean.nc'
'sdfopen E:\vwnd.mon.mean.nc'
'sdfopen E:\pres.mon.mean.nc'
'set grads off'
'set grid off'
'set lon  20 90'
'set lat  0  70'
'set z 1 8'
'define shum1=ave(shum.1,t=450,t=816,12)'
'define shum2=ave(shum.1,t=451,t=816,12)'
'define shum3=ave(shum.1,t=452,t=816,12)'
'define q=(shum1+shum2+shum3)/3'
'define uwnd1=ave(uwnd.2,t=450,t=816,12)'
'define uwnd2=ave(uwnd.2,t=451,t=816,12)'
'define uwnd3=ave(uwnd.2,t=452,t=816,12)'
'define u=(uwnd1+uwnd2+uwnd3)/3'
'define vwnd1=ave(vwnd.3,t=450,t=816,12)'
'define vwnd2=ave(vwnd.3,t=451,t=816,12)'
'define vwnd3=ave(vwnd.3,t=452,t=816,12)'
'define v=(vwnd1+vwnd2+vwnd3)/3'
'define pres1=ave(pres.4,t=450,t=816,12)'
'define pres2=ave(pres.4,t=451,t=816,12)'
'define pres3=ave(pres.4,t=452,t=816,12)'
'define p=(pres1+pres2+pres3)/3'
'define qu=u*q*1000'
'define qv=v*q*1000'
'define qus=vint(p(lev=0)/100,qu,300)/9.8'
'define qvs=vint(p(lev=0)/100,qv,300)/9.8'
'set z 1'
'define qum=ave(qus,t=1,t=31,1)'
'define qvm=ave(qvs,t=1,t=31,1)'
'difine sq=mag(qum,qvm) '
'set gxout vector'
'd qum;qvm'
'set gxout contour'
'd sq'
'printim E:\div2\quv.png'
'reinit'
;

LNYZ

{:5_275:}{:5_275:}

18393810639

好吧 。我是为了金币的。