水汽通量散度试了论坛所有去白边的方法还是不行啊，，大神急求

xrpamhf

'sdfopen D:\sandu\uwnd.mon.mean.nc'
'sdfopen D:\sandu\vwnd.mon.mean.nc'
'sdfopen D:\sandu\rhum.mon.mean.nc'
'sdfopen D:\sandu\air.mon.mean.nc'
'set lon 90 140'
'set lat 10 60'
'set lev 600'
i=169

while(i<=170)
'set t 'i
'set grads off'
'set grid off'
'set gxout shaded'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'd hdivg(q*uwnd.1/g,q*vwnd.2/g)'
'run cbarn'
'q w2xy 100.37 38.85'
'set lon 80 130'
'set lat 20 50'
'draw shp D:\xueshuupan\heihe\province'
'print D:\sandu\9.eps'
'c'
i=i+1
endwhile
;

xrpamhf

在线等，特别着急

xrpamhf

1099221723

'reinit'
'sdfopen D:\sandu\uwnd.mon.mean.nc'
'sdfopen D:\sandu\vwnd.mon.mean.nc'
'sdfopen D:\sandu\rhum.mon.mean.nc'
'sdfopen D:\sandu\air.mon.mean.nc'
'set lon 79.99 130.01'
'set lat 19.99 50.01'
'set lev 600'
i=169

while(i<=170)
'set t 'i
'set grads off'
'set grid off'
'set gxout shaded'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'd hdivg(q*uwnd.1/g,q*vwnd.2/g)'
'run cbarn'
'q w2xy 100.37 38.85'
'set lon 79.99 130.01'
'set lat 19.99 50.01'
#'draw shp D:\xueshuupan\heihe\province'
'set mpdset cnworld'
#'print D:\sandu\9.eps'
'printim D:\sandu\1.png white'
#'c'
i=i+1
endwhile
;

xrpamhf

1099221723 发表于 2019-4-3 18:53
'reinit'
'sdfopen D:\sandu%uwnd.mon.mean.nc'
'sdfopen D:\sandu\vwnd.mon.mean.nc'

好的老兄，这会还在做实验，完了第一时间就试试你修改的这个，真的特别特别感谢！！！！！

xrpamhf

1099221723 发表于 2019-4-3 18:53
'reinit'
'sdfopen D:\sandu%uwnd.mon.mean.nc'
'sdfopen D:\sandu\vwnd.mon.mean.nc'

果然可以了灰常感谢

qxlong

请问楼主，你画的这个散度是月平均的？

xrpamhf

qxlong 发表于 2019-4-13 16:07
请问楼主，你画的这个散度是月平均的？

是的，用的月均值