一段时间平均的整层水汽通量，图出来了，不知道正不正确

太阳雨

在论坛上搜集了大神们关于Grada画整层水汽通量的好多脚本，自己写了一个，是计算一段降雨时间的平均的整层水汽通量，把图和gs脚本贴出来给大家看一看，在计算每一层的水汽通量的时候以及再算平均的时候没用do while之类的循环语句，所以不确定我的图是不是一定正确，希望有大神能指点一下，如果是正确的话大家可以借鉴一下，我自己也还在验证这个图的正确性，谢谢大家！！！！！

16天平均的整层水汽通量

太阳雨

有错别字，是Grads
我直接吧脚本贴出来吧，方便浏览
'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 -20 80'
'set lon 40 160'
'set z 1 8'
'set t 163 178'
'set grid off'
'set xlint 20'
'set ylint 10'
'set grads off'
*'set timelab off'
*'set mpdset mres'
*'set map 1 1 6'
'set mpdset cnworld'
'define qu=uwnd*shum.3'
'define qv=vwnd.2*shum.3'
'define qus=vint(pres.4(lev=0)/100,qu,700)'
'define qvs=vint(pres.4(lev=0)/100,qv,700)'
'set t 163'
'define qum=ave(qus,t=163,t=178,1)'
'define qvm=ave(qvs,t=163,t=178,1)'
'set z 1'
'define a=sqrt(qum*qum+qvm*qvm)'
'set gxout shaded'
'set cmin 150'
'd a'

'run cbarn.gs'
'set arrscl 0.5 700'
'set gxout stream'
'set strmden 4'

'd qum;qvm'
'printim D:\9801\1jieko\11.png white'
*'disable print'
*'disable fwrite'
;

river

太阳雨 发表于 2013-7-4 09:34

                               
登录/注册后可看大图

有错别字，是Grads
我直接吧脚本贴出来吧，方便浏览
'reinit'

看着的话没有什么问题。不过有一个疑问
'define qus=vint(pres.4(lev=0)/100,qu,700)'
'define qvs=vint(pres.4(lev=0)/100,qv,700)'
你的pres资料只有一层还是有多层啊？一般来说，用z 代替高度层的时候使用的是1 2 3 ····
用lev的时候要写他实际的那层气压，比如1000  925  850   700·····

太阳雨

river 发表于 2013-7-4 17:47

                               
登录/注册后可看大图

看着的话没有什么问题。不过有一个疑问
'define qus=vint(pres.4(lev=0)/100,qu,700)'
'define qvs=vin ...

pres只有一层的

ctlinfo这样写的
z 1 linear 0 0

pres 0 -999 mean daily presssure at surface

我担心的主要是循环有木有问题，还有这个应该是lev=0还是lev=1的问题

river

太阳雨 发表于 2013-7-5 11:19

                               
登录/注册后可看大图

pres只有一层的

ctlinfo这样写的

你这个不需要循环，看着没有什么问题。就是你那个lev的问题，用lev的话就要具体写气压的值，比如你这个只有地面一层，应该是lev=1000，也可以写成z=1的形式。

DQ153LX

我运行了一下  出不了图   麻烦问一下   这是什么情况？