夏季整层平均水汽通量出图的问题

Raina

我是非气象专业，研究GrADS一段时间了，我想做一个2015年夏季(6-8月)整层（surface to 300hpa）水汽通量图，不需要阴影图，这是我的gs文件，以及出的图，我感觉方向和单位都不对。请各位高手帮我看看哪里出了问题。不胜感激。
'reinit'
'sdfopen E:\wyr\uwnd.2015.nc'
'sdfopen E:\wyr\vwnd.2015.nc'
'sdfopen E:\wyr\shum.2015.nc'
'sdfopen E:\wyr\pres.sfc.2015.nc'
'set lat 0 50'
'set lon 50 140'
'set z 1 8'
'set t 180 270'
'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'
'define qus=vint(pres.4(lev=0)/100,qu,300)/9.8'
'define qvs=vint(pres.4(lev=0)/100,qv,300)/9.8'
'define qum=ave(qus,t=180,t=270,12)'
'define qvm=ave(qvs,t=180,t=270,12)'
'set z 1'
'set t 270'
'd qum;qvm'
'printim E:\wyr1115.png'
;

Raina

我觉得从地面到300hpa的积分我好像没有求对 我用q dims 得到的是lve是1000

Raina

但是6月-8月的时间应该没有问题 是180-270

D.B.L

http://bbs.06climate.com/forum.p ... p;extra=&page=1
你看一下这个就明白了

Raina

D.B.L 发表于 2015-11-16 08:24
http://bbs.06climate.com/forum.php?mod=viewthread&tid=6383&extra=&page=1
你看一下这个就明白了

这个帖子我看了好多次 是不用ave么？

D.B.L

你最后用ave没问题，前面的计算不对啊
你自己琢磨一下

D.B.L

http://bbs.06climate.com/forum.php?mod=viewthread&tid=4920
你在参考一下这个

Raina

'reinit'
'sdfopen E:\wyr\uwnd.2015.nc'
'sdfopen E:\wyr\vwnd.2015.nc'
'sdfopen E:\wyr\shum.2015.nc'
'sdfopen E:\wyr\pres.sfc.2015.nc'
'set lat 0 40'
'set lon 80 140'
'set map 1 1 1'
'draw map'
'set cthick 4'
'set xlopts 1 4 0.12'
'set ylopts 1 4 0.12'
'set z 1 8'
'set t 152 242'
'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'
'define qus=vint(pres.4(lev=0)/100,qu,300)/9.8'
'define qvs=vint(pres.4(lev=0)/100,qv,300)/9.8'
'set z 1'
'set t 152'
'd qus;qvs'
'printim E:\wyr2.png x1000y800 white'
;

Raina

D.B.L 发表于 2015-11-16 09:18
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4920
你在参考一下这个

谢谢！我还是不太清楚前面哪里出了问题 你能帮我看看么 8楼我把新的脚本和图贴出来了

D.B.L

Raina 发表于 2015-11-16 09:25
谢谢！我还是不太清楚前面哪里出了问题 你能帮我看看么 8楼我把新的脚本和图贴出来了

'reinit'
'sdfopen E:\wyr\uwnd.2015.nc'
'sdfopen E:\wyr\vwnd.2015.nc'
'sdfopen E:\wyr\shum.2015.nc'
'sdfopen E:\wyr\pres.sfc.2015.nc'
'set lat 0 40'
'set lon 80 140'
'set map 1 1 1'
'draw map'
'set cthick 4'
'set xlopts 1 4 0.12'
'set ylopts 1 4 0.12'
'set z 1 8'
'set t 152 242'
'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'

'set z 1'
'set t 152'
'define qus=vint(pres.4(lev=0)/100,qu,300)/9.8'
'define qvs=vint(pres.4(lev=0)/100,qv,300)/9.8'

'd qus;qvs'
'printim E:\wyr2.png x1000y800 white'
;
你试试