grads计算水汽通量求助！！！！（已解决，并更新）

太阳雨

wet510 发表于 2012-10-20 12:09

                               
登录/注册后可看大图

这个可以用ave函数来实现啊，我在14楼贴的程序时间维是变化的，从1998年6月1日12时到1998年8月31日12时， ...

十分感谢楼主的指教，我现在把修改过后的gs脚本和图贴出来，还希望楼主帮我看看是否正确，还有就是右下角图标单位的量级有点没搞清楚，我现在的单位是不是g/(m*s)?麻烦楼主再给指点一下，感谢~
'reinit'
'sdfopen D:\9801\uwnd.1998.nc'
'sdfopen D:\9801\vwnd.1998.nc'
'sdfopen D:\9801\shum.1998.nc'
'sdfopen D:\9801\pres.sfc.1998.nc'
'set lat -10 70'
'set lon 40 160'
'set z 1 8'
'set t 163 178'
'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'
'define qus=vint(pres.4(lev=0)/100,qu,700)/9.8'
'define qvs=vint(pres.4(lev=0)/100,qv,700)/9.8'
'define qum=ave(qus,t=163,t=178,1)'
'define qvm=ave(qvs,t=163,t=178,1)'
'set z 1'
'set t 163'
'd qum;qvm'
'printim D:\9801\1jieko.png'
;

wet510

太阳雨 发表于 2012-10-20 20:06

                               
登录/注册后可看大图

十分感谢楼主的指教，我现在把修改过后的gs脚本和图贴出来，还希望楼主帮我看看是否正确，还有就是右下角 ...

脚本没什么问题，水汽通量的单位也是对的

zdf2075

太阳雨 发表于 2012-10-20 20:06

                               
登录/注册后可看大图

十分感谢楼主的指教，我现在把修改过后的gs脚本和图贴出来，还希望楼主帮我看看是否正确，还有就是右下角 ...

'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'

这里为什么要乘以1000？ncep里面的比湿单位不是g/kg吗？

kong369

学习学习哦！！~~~~~

wet510

zdf2075 发表于 2013-3-14 20:18

                               
登录/注册后可看大图

'define qu=uwnd*shum.3*1000'
'define qv=vwnd.2*shum.3*1000'

关于单位可以参考这篇帖子，具体情况看个人

整层水汽通量散度计算及其单位
http://bbs.06climate.com/forum.p ... 20&fromuid=6323

coco12345

赞一个 正好碰到这个问题~~~

coco12345

好在有高手指点啊，太感谢啦

星河~~

好贴~~学到了

煎豆

刚好要画这张图。。太好了

phys

好贴！！我也试试