急急急急！！！！多年平均水通量汽

WOT

本人想计算1961~1986年500hpa的多年平均水汽通量，我先计算1961-2020年12个月的与平均水汽通量，然后加和计算，画出的图里面的箭头很短，求大神指点一下，这究竟是怎么回事啊，还有就是hdivg算的是水汽通量还是水汽通量散度（改脚本改的有点懵了）？
改脚本改的有点懵了
数据使用为ncep月平均数据，其中1961年1月t=157,1986年1月t=457,
画出后的脚本、图像如下

'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'

'set parea 1.5 10.0 1.5 8.0'
'set map 1 1 3'
'set mpdset cnworld cnriver'
'set lev 500'
'set xlopts -10'
'set ylopts -10'
'set xlint 10'
'set ylint 10'

'set t 157 468'
'define qu1=uwnd*shum.3/1000'
'define qv1=vwnd.2*shum.3/1000'
'set t 157'
'define qud1=ave(qu1,t+0,t=468,12)'
'define qvd1=ave(qv1,t+0,t=468,12)'
'define qdiv1=hdivg(qud1,qvd1)'

'set t 158 468'
'define qu2=uwnd*shum.3/1000'
'define qv2=vwnd.2*shum.3/1000'
'set t 158'
'define qud2=ave(qu2,t+0,t=468,12)'
'define qvd2=ave(qv2,t+0,t=468,12)'
'define qdiv2=hdivg(qud2,qvd2)'

'set t 159 468'
'define qu3=uwnd*shum.3/1000'
'define qv3=vwnd.2*shum.3/1000'
'set t 159'
'define qud3=ave(qu3,t+0,t=468,12)'
'define qvd3=ave(qv3,t+0,t=468,12)'
'define qdiv3=hdivg(qud3,qvd3)'

'set t 160 468'
'define qu4=uwnd*shum.3/1000'
'define qv4=vwnd.2*shum.3/1000'
'set t 160'
'define qud4=ave(qu4,t+0,t=468,12)'
'define qvd4=ave(qv4,t+0,t=468,12)'
'define qdiv4=hdivg(qud4,qvd4)'

'set t 161 468'
'define qu5=uwnd*shum.3/1000'
'define qv5=vwnd.2*shum.3/1000'
'set t 161'
'define qud5=ave(qu5,t+0,t=468,12)'
'define qvd5=ave(qv5,t+0,t=468,12)'
'define qdiv5=hdivg(qud5,qvd5)'

'set t 162 468'
'define qu6=uwnd*shum.3/1000'
'define qv6=vwnd.2*shum.3/1000'
'set t 162'
'define qud6=ave(qu6,t+0,t=468,12)'
'define qvd6=ave(qv6,t+0,t=468,12)'
'define qdiv6=hdivg(qud6,qvd6)'

'set t 163 468'
'define qu7=uwnd*shum.3/1000'
'define qv7=vwnd.2*shum.3/1000'
'set t 163'
'define qud7=ave(qu7,t+0,t=468,12)'
'define qvd7=ave(qv7,t+0,t=468,12)'
'define qdiv7=hdivg(qud7,qvd7)'

'set t 164 468'
'define qu8=uwnd*shum.3/1000'
'define qv8=vwnd.2*shum.3/1000'
'set t 164'
'define qud8=ave(qu8,t+0,t=468,12)'
'define qvd8=ave(qv8,t+0,t=468,12)'
'define qdiv8=hdivg(qud8,qvd8)'

'set t 165 468'
'define qu9=uwnd*shum.3/1000'
'define qv9=vwnd.2*shum.3/1000'
'set t 165'
'define qud9=ave(qu9,t+0,t=468,12)'
'define qvd9=ave(qv9,t+0,t=468,12)'
'define qdiv9=hdivg(qud9,qvd9)'

'set t 166 468'
'define qu10=uwnd*shum.3/1000'
'define qv10=vwnd.2*shum.3/1000'
'set t 166'
'define qud10=ave(qu10,t+0,t=468,12)'
'define qvd10=ave(qv10,t+0,t=468,12)'
'define qdiv10=hdivg(qud10,qvd10)'

'set t 167 468'
'define qu11=uwnd*shum.3/1000'
'define qv11=vwnd.2*shum.3/1000'
'set t 167'
'define qud11=ave(qu11,t+0,t=468,12)'
'define qvd11=ave(qv11,t+0,t=468,12)'
'define qdiv11=hdivg(qud11,qvd11)'

'set t 168 468'
'define qu12=uwnd*shum.3/1000'
'define qv12=vwnd.2*shum.3/1000'
'set t 168'
'define qud12=ave(qu12,t+0,t=468,12)'
'define qvd12=ave(qv12,t+0,t=468,12)'
'define qdiv12=hdivg(qud12,qvd12)'


'define qud=qud1+qud2+qud3+qud4+qud5+qud6+qud7+qud8+qud9+qud10+qud11+qud12'
'define qvd=qvd1+qvd2+qvd3+qvd4+qvd5+qvd6+qvd7+qvd8+qvd9+qvd10+qvd11+qvd12'

'define qdiv=qdiv1+qdiv2+qdiv3+qdiv4+qdiv5+qdiv6+qdiv7+qdiv8+qdiv9+qdiv10+qdiv11+qdiv12'

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

'set gxout shaded'

*'E:\GrADS\test\sebiao.gs'
'd qdiv'
'run cbarn'
*'E:\GrADS\test\cbar_matlab 10 1 0'

'set gxout barb'
'set ccolor 1'
'set gxout vector'
'set arrscl 0.5 200'
'set arrowhead 0.04'

'd qud;qvd'
*'cbarn 1.2 0'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type b -interval 20 -position o -suffix `3.'
'run C:\OpenGrADS\Classic\scripts\axis.gs -type L -interval 10 -position o -suffix `3.'

'draw title 1961-1986 500hpa vapor flux '

'printim E:\GrADS\test\88.png  white x2500 y2000'
'disable print'

pefuen

调一下箭头大小

pefuen

hdivg算的是水汽通量散度，mag（q*u/g,q*v/g）才是水汽通量，写代码的时候命名就要标好画的哪个要素，这样才不会乱

pefuen

给你看看我之前画的吧，你的色标也太长了，建议慢慢调整一下

WOT

pefuen 发表于 2023-2-23 15:40
给你看看我之前画的吧，你的色标也太长了，建议慢慢调整一下

大神求一个脚本行吗？

pefuen

WOT 发表于 2023-2-23 15:49
大神求一个脚本行吗？

论坛有很多水汽通量的脚本，我当时也是照着改的，先把东西画出来，哪里不对调哪里。我看了你的代码，太复杂了，我的脚本不适合你，自己摸索会好一点哈

WOT

pefuen 发表于 2023-2-23 16:02
论坛有很多水汽通量的脚本，我当时也是照着改的，先把东西画出来，哪里不对调哪里。我看了你的代码，太复 ...

好的好的，那水汽通量是用hdivg算的吗？