求助：绘制整层水汽通量和散度却只能画出一小部分图像

WOT

'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\pres.mon.mean.nc'
'set z 1 8'
'set lon 70 150'
'set lat 10 70'


*'set t 161'
*'set t  801'
'define qud5=ave(uwnd*shum.3/9.8,t=1,t=859)'
'define qvd5=ave(vwnd.2*shum.3/9.8,t=1,t=859)'
*statr integral from pres
'set z 1'
'define qus=vint(pres.4(lev=0),qud5,300)'
'define qvs=vint(pres.4(lev=0),qvd5,300)'
'define div=hdivg(qus,qvs)'


'set mpdset cnworld cnriver'
'set map 1 1 3'
'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 10'
'set ylint 10'
'set lon 60 150'
'set lat 10 60'
'set gxout shaded'
'd div'
'set gxout vector'
'd qus;qvs'
*adjust the colorbar location
'cbarn 1.2 0'
'draw title 1000-300hpa vapor flux and divergence'
'printim E:\GrADS\test\9.png white x2000 y1800'
'disable print'

river

看脚本也没看出啥大问题，你把两个经纬度的设置统一一下，要不就把最前面那个经纬度设置删掉

WOT

river 发表于 2023-1-7 13:38
看脚本也没看出啥大问题，你把两个经纬度的设置统一一下，要不就把最前面那个经纬度设置删掉

还是和原来一样

栤唥de嶶笑

这个好像有点问题，pres只有地面一层吗，pres.4(lev=0)是指从高空0hPa积分到300hPa？从1000hPa积分应该是pres.4(lev=1000)？
'define qus=vint(pres.4(lev=0),qud5,300)'
'define qvs=vint(pres.4(lev=0),qvd5,300)'

另外，你这个公式'define qvd5=ave(vwnd.2*shum.3/9.8,t=1,t=859)'应该存在一些问题，具体看看我的帖子 http://bbs.06climate.com/forum.p ... 2709&extra=page%3D1

最后，计算时间平均的垂直积分，建议逐月计算，再取平均，这样不会出错，尤其对于里面有指数等复杂运算的

WOT

栤唥de嶶笑 发表于 2023-1-9 17:07
这个好像有点问题，pres只有地面一层吗，pres.4(lev=0)是指从高空0hPa积分到300hPa？从1000hPa积分应该是pr ...

您好，我重新在论坛里找了一个脚本，按照您之前说的pres.4(lev=0)更改为pres.4(lev=1000)，但是没法出图了，保持原来的pres.4(lev=0)是可以出图的，您再帮我看看是我脚本有问题吗？
'reinit'
'sdfopen E:\GrADS\testwww\uwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\vwnd.mon.mean.nc'
'sdfopen E:\GrADS\testwww\shum.mon.mean.nc'
'sdfopen E:\GrADS\testwww\pres.mon.mean.nc'

'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 3'
'set mpdset cnworld cnriver'

'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 10'
'set ylint 10'

'set t 413 869'
'set z 1 8 '
'define qu5=uwnd*shum.3/1000'
'define qv5=vwnd.2*shum.3/1000'
'set t 413'
'set z 1 8 '
'define qud5=ave(qu5,t+0,t=864,12)'
'define qvd5=ave(qv5,t+0,t=864,12)'
'define qdiv5=hdivg(qud5,qvd5)'

'set t 414 870'
'set z 1 8 '
'define qu6=uwnd*shum.3/1000'
'define qv6=vwnd.2*shum.3/1000'
'set t 414'
'set z 1 8 '
'define qud6=ave(qu6,t+0,t=864,12)'
'define qvd6=ave(qv6,t+0,t=864,12)'
'define qdiv6=hdivg(qud6,qvd6)'

'set t 415 871'
'set z 1 8 '
'define qu7=uwnd*shum.3/1000'
'define qv7=vwnd.2*shum.3/1000'
'set t 415'
'set z 1 8 '
'define qud7=ave(qu7,t+0,t=864,12)'
'define qvd7=ave(qv7,t+0,t=864,12)'
'define qdiv7=hdivg(qud7,qvd7)'

'set t 416 872'
'set z 1 8 '
'define qu8=uwnd*shum.3/1000'
'define qv8=vwnd.2*shum.3/1000'
'set t 416'
'set z 1 8 '
'define qud8=ave(qu8,t+0,t=864,12)'
'define qvd8=ave(qv8,t+0,t=864,12)'
'define qdiv8=hdivg(qud8,qvd8)'

'set z 1 8 '
'define qud=(qud5+qud6+qud7+qud8)/4'
'define qvd=(qvd5+qvd6+qvd7+qvd8)/4'
'define qdiv=(qdiv5+qdiv6+qdiv7+qdiv8)/4'

'set z 1'
'define qus=vint(pres.4(lev=1000),qud,300)'
'define qvs=vint(pres.4(lev=1000),qvd,300)'

'define shuiqi=mag(qus,qvs)'

*'define sandu=hdivg(qus,qvs)'


'set lon 40 120'
'set lat 20 60'

'set gxout shaded'
'd shuiqi'
'set gxout barb'
'set ccolor 1'
'set gxout vector'
'd qus;qvs'
'cbarn 1.2 0'

'draw title 1000-300hpa vapor flux and divergence'

'printim E:\GrADS\test\44.png  white x2000 y1800'
'disable print'