如何用NCL从地表气压开始进行垂直积分

奚子惠

使用       p_file=addfile("pres.nc","r")
      pres=p_file->pres
      p       = (/ 1000.,925.,850.,700.,600.,500., \
               400.,300.,250.,200.,150.,100., \
                70.,50.,30.,20.,10. /)
     linlog = 1
     pbot   = 1000.                                             
     ptop   = 300.

     u = vibeta (p,qu_new(time|:,lat|:,lon|:,level|:),linlog,pres,pbot,ptop)/9.8 进行垂直积分，发现这个函数是从pbot积到ptop，但是我想从地表气压pres往上积分，请教各位大神们如何解决？

风往北吹

vibeta
Performs vertical integration using beta factors.

Prototype

        function vibeta (
                p          : numeric,  
                x          : numeric,  
                linlog [1] : integer,  
                psfc       : numeric,  
                pbot   [1] : numeric,  
                ptop   [1] : numeric   
        )

        return_val [dimsizes(psfc)] :  numeric
Arguments

p
A multi-dimensional array of pressure levels. The rightmost dimension must be at least of length 3. The order is bottom-to-top.

x
A multi-dimensional numeric array to be integrated. If p is a 1-dimensional array, then the rightmost dimension of x must be same size as p. Otherwise, x and p must be same size. The order of x's level dimension must be bottom-to-top.

linlog
A scalar integer, 1 = linear interpolation, 2 = log interpolation.

psfc
A multi-dimensional numerical array of surface pressures. Must be same size as x, minus the rightmost dimension.

pbot
A scalar numeric equal to the lower limits of integration.

ptop
A scalar numeric equal to the upper limits of integration.

风往北吹

psfc
A multi-dimensional numerical array of surface pressures. Must be same size as x, minus the rightmost dimension.

niuda

  为什么要 “ /9.8” ?

talkd

风往北吹 发表于 2015-6-6 07:42
psfc
A multi-dimensional numerical array of surface pressures. Must be same size as x, minus the ri ...

这个是不是只是为了判断哪里是missingvalue用的。实质上就是从pbot开始积分的？

talkd

niuda 发表于 2015-6-6 08:19
为什么要 “ /9.8” ?

z坐标和p坐标换算

niuda

talkd 发表于 2015-6-6 10:40
z坐标和p坐标换算

  一般不是  /(pbot-ptop) 的吗？

另外如果想从地表气压pres往上积分的话，是否可以设置 pbot的值很大，这样是否会默认从地表气压开始积分？

talkd

niuda 发表于 2015-6-6 11:25
一般不是  /(pbot-ptop) 的吗？

另外如果想从地表气压pres往上积分的话，是否可以设置 pbot的值很 ...

好主意，我觉得值得一试

奚子惠

niuda 发表于 2015-6-6 08:19
为什么要 “ /9.8” ?

因为我要算大气质量流，公式里面有要/g的哦

奚子惠

talkd 发表于 2015-6-6 13:08
好主意，我觉得值得一试

我来试试~~