整层水汽通量的计算（根据自己的1961-2014夏季平均）necp资料

choityyy

啊 不好意思我不太会发图片 您方便留个联系方式么？邮箱或者其他都可以  

帅帅f

choityyy 发表于 2019-12-24 19:51
啊 不好意思我不太会发图片 您方便留个联系方式么？邮箱或者其他都可以

你再用用我新上传的第二个脚本

choityyy

帅帅f 发表于 2019-12-24 22:02
你再用用我新上传的第二个脚本

楼主，你好！我试用了您新上传的danceng(ture).gs 现在箭头的方向是对的。但是量级还是不太对，文献上的箭头是150，我画出来箭头是20，我看您的gs是算单层850hPa的，文献里是算整层的，会不会和这个有关系？您知道整层的应该怎么修改gs吗？

帅帅f

choityyy 发表于 2019-12-25 20:25
楼主，你好！我试用了您新上传的danceng(ture).gs 现在箭头的方向是对的。但是量级还是不太对，文献上的 ...

肯定会不一样的 单层和整层的  整层的脚本我上传到了 第三个链接  你要注意的是我算的时间是一整个的时间段，如果你只是单单计算夏季平均  时间循环可能要改变一下

choityyy

帅帅f 发表于 2019-12-26 14:32
肯定会不一样的 单层和整层的  整层的脚本我上传到了 第三个链接  你要注意的是我算的时间是一整个的时间 ...

楼主，你好！我刚刚试了第三个文件，整层的gs，整体没太大毛病。我想要挑出夏季的数据进行计算，先用您最早发过的gs文件里的循环，图画出来的画风场方向还是反的，我又试了您单层水汽通量里的循环方法，gs代码如下：
*--------1948.6-2006.12-----------
'set t 6 708'
'define qu5=uwnd*shum.3/9.8'
'define qv5=vwnd.2*shum.3/9.8'
'set t 6'
'define qud5=ave(qu5,t+0,t=708,12)'
'define qvd5=ave(qv5,t+0,t=708,12)'
'define qus5=vint(pres.4(lev=0),qud5,300)'
'define qvs5=vint(pres.4(lev=0),qvd5,300)'
*'define qdiv5=hdivg(qud5,qvd5)'
'define sandu5=hdivg(qus5,qvs5)'

'set t 7 708'
'define qu6=uwnd*shum.3/9.8'
'define qv6=vwnd.2*shum.3/9.8'
'set t 7'
'define qud6=ave(qu6,t+0,t=708,12)'
'define qvd6=ave(qv6,t+0,t=708,12)'
'define qus6=vint(pres.4(lev=0),qud6,300)'
'define qvs6=vint(pres.4(lev=0),qvd6,300)'
*'define qdiv6=hdivg(qud6,qvd6)'
'define sandu6=hdivg(qus6,qvs6)'

'set t 8 708'
'define qu7=uwnd*shum.3/9.8'
'define qv7=vwnd.2*shum.3/9.8'
'set t 8'
'define qud7=ave(qu7,t+0,t=708,12)'
'define qvd7=ave(qv7,t+0,t=708,12)'
'define qus7=vint(pres.4(lev=0),qud7,300)'
'define qvs7=vint(pres.4(lev=0),qvd7,300)'
*'define qdiv7=hdivg(qud7,qvd7)'
'define sandu7=hdivg(qus7,qvs7)'

'define qus=(qus5+qus6+qus7)/3'
'define qvs=(qvs5+qvs6+qvs7)/3'
'define qdiv=(sandu5+sandu6+sandu7)/3'

很奇怪的一点是 这个gs一运行opengrads就会闪退。我试了下，如果只循环6，7月的话，就不会闪退。请楼主指教闪退原因，或者能否告诉我怎样编写挑出夏季数据的循环，万分感谢！

choityyy

choityyy 发表于 2019-12-26 20:32
楼主，你好！我刚刚试了第三个文件，整层的gs，整体没太大毛病。我想要挑出夏季的数据进行计算，先用您最 ...

再补充一点，如果删掉6月份的循环也是可以的，好像这个gs只能承担两个月的计算，多一个月就会闪退

Donald

学习一下，谢谢分享！！！

duruyi

大佬，会不会改nc文件里面的变量名

yqf13414898540

楼主你好，想请问你1961-2014年夏季平均值是怎么求的？

居一桠

帅帅f 发表于 2019-8-29 22:04
这个是ctl文件的写法dset F:\grads_vapor_flux2%uwnd.grd
undef  -9.99e+8
xdef           144 linear  0 ...

楼主你好，问一下如果是冬季的话，gs这里t='i',t='i+4'要怎么改呀，怎么改成12，1，2月份呀