求助：绘制整层水汽通量与文献中图像不符

WOT

这个脚本是绘制1982-2020年5~8月的整层（1000-300hpa）水汽通量脚本，这个脚本可以绘制出图像，图像如图一所示，但是无法与文献中的图像对应上，请各位帮我看一下，哪里出错了数据时间范围是1948年1月到2022年11月

'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=869,12)'
'define qvd5=ave(qv5,t+0,t=869,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=870,12)'
'define qvd6=ave(qv6,t+0,t=870,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=871,12)'
'define qvd7=ave(qv7,t+0,t=871,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=872,12)'
'define qvd8=ave(qv8,t+0,t=872,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=0),qud,300)'
'define qvs=vint(pres.4(lev=0),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'

tulalang

感觉大体轮廓也还是能对的上，你要不用调色盘弄个和文献一样的渐变色试试看？

tbag

确实大体上对得上但是分辨率感觉很不一样

伽蓝鸟

计算水汽通量应该用逐日或者6小时数据，用月平均数据算出来的当然不准确啦

WOT

tulalang 发表于 2023-1-12 22:50
感觉大体轮廓也还是能对的上，你要不用调色盘弄个和文献一样的渐变色试试看？

你好，怎么弄调色盘呢

WOT

伽蓝鸟 发表于 2023-1-13 09:02
计算水汽通量应该用逐日或者6小时数据，用月平均数据算出来的当然不准确啦

气压也要用逐日数据吗？

tulalang

具体学习可参考这个网址http://bbs.06climate.com/forum.p ... &extra=page%3D1
以下是我看着你发的文献色标设置的代码，grads很久没用，可能设置的细节有点问题
'set rgb 20 0 4 247'
'set rgb 21 1 16 231'
'set rgb 22 1 27 216'
'set rgb 23 2 39 200'
'set rgb 24 3 51 184'
'set rgb 25 3 63 168'
'set rgb 26 4 75 153'
'set rgb 27 5 86 138'
'set rgb 28 5 98 122'
'set rgb 29 6 110 107'
'set rgb 30 13 119 97'
'set rgb 31 20 128 87'
'set rgb 32 27 137 77'
'set rgb 33 33 146 67'
'set rgb 34 42 155 59'
'set rgb 35 52 162 52'
'set rgb 36 61 170 45'
'set rgb 37 71 177 38'
'set rgb 38 81 185 31'
'set rgb 39 92 191 26'
'set rgb 40 104 197 21'
'set rgb 41 116 203 16'
'set rgb 42 127 209 11'
'set rgb 43 140 214 8'
'set rgb 44 153 217 7'
'set rgb 45 167 221 6'
'set rgb 46 180 224 5'
'set rgb 47 194 228 4'
'set rgb 48 206 230 3'
'set rgb 49 218 232 2'
'set rgb 50 230 234 1'
'set rgb 51 242 236 0'
'set rgb 52 249 234 0'
'set rgb 53 250 229 0'
'set rgb 54 251 223 0'
'set rgb 55 253 218 0'
'set rgb 56 254 212 0'
'set rgb 57 254 206 0'
'set rgb 58 254 199 0'
'set rgb 59 253 193 0'
'set rgb 60 253 186 0'
'set rgb 61 253 179 0'
'set rgb 62 253 170 0'
'set rgb 63 253 161 0'
'set rgb 64 253 152 0'
'set rgb 65 253 143 0'
'set rgb 66 253 133 0'
'set rgb 67 253 122 0'
'set rgb 68 252 112 0'
'set rgb 69 252 101 0'
'set rgb 70 251 92 1'
'set rgb 71 249 83 3'
'set rgb 72 246 75 4'
'set rgb 73 244 66 6'
'set rgb 74 242 58 8'
'set rgb 75 238 51 12'
'set rgb 76 234 44 15'
'set rgb 77 230 37 19'
'set rgb 78 226 31 22'
'set rgb 79 221 24 26'
'set rgb 80 215 19 31'
'set rgb 81 208 14 36'
'set rgb 82 202 9 41'
'set rgb 83 196 4 46'
'set rgb 84 186 3 54'
'set rgb 85 176 2 62'
'set rgb 86 166 1 71'
'set rgb 87 156 0 79'
'set rgb 88 146 0 87'
'set rgb 89 134 0 96'
'set rgb 90 123 0 105'
'set rgb 91 112 0 114'

'set ccols 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 '

'set clevs 0   5  10  15  20  25  30  35  40  45  50  55  60  65  70  75  80  85  90  95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185 190 195 200 205 210 215 220 225 230 235 240 245 250 255 260 265 270 275 280 285 290 295 300 305 310 315 320 325 330 335 340 345 350'

WOT

tulalang 发表于 2023-1-16 09:41
具体学习可参考这个网址http://bbs.06climate.com/forum.php?mod=viewthread&tid=9294&extra=page%3D1
以 ...

万分感谢！！！