关于NCL的垂直积分函数Vibeta的困惑

青青草心

问题描述：在计算青藏高原地区整层水汽通量时，需要从地表气压Ps（t,y,x）向上垂直积分，以前使用了grads的vint函数，这个函数可以从Ps开始积分（函数使用方法：http://cola.gmu.edu/grads/gadoc/gradfuncvint.html）。NCL中垂直积分函数为vibeta，这个函数给出的例子都是具体的等压面P1到P2积分，虽然输入了地表气压数据，但是不知道高原地表以下的数据有没有参与积分过程？[size=13.3333px]比如vint = vibeta[size=13.3333px](p,x(time|:,lat|:,lon|:,lev|:),1,psfc,1000,300)表示从1000积分到300hPa,那么高原地区的水汽通量是从Ps还是1000向上积分的？（使用的欧洲中心再分析数据在高原地区的850hPa也是有数值的）

[size=13.3333px]翻看了NCL官网也没有理解，所以发帖请各位坛友来指点迷津！

卢小傻子

可以用另一个函数，官网上也指出了  实现从地面开始积分吧
This matches an alternative approach
                                             ; dp = layer thickness
  dp   = dpres_plevel(p, psfc, ptop, 0)   
                                             ; weighted sum/sum_of_layer_thickness
  fint = dim_sum_n(x*dp, 0)/dim_sum_n(dp,0)

  print(dp)  ; (/50.5, 75, 112.5, 125, 100, 100, 100, 75, 50, 50, 40, 25, 20, 15, 10, 5/)
  print("dim_sum_n(dp,0)="+dim_sum_n(dp,0))  ; 1003
  print("fint="+fint)                        ; fint ===> 254.605 K

cici破晓碎寒

请问楼主解决了吗

ANNA0217

vibeta提示：If psfc is not a constant, then its dimension sizes must be equal to all but the rightmost dimension sizes of x。这个怎么解决啊

lt5201314

ANNA0217 发表于 2020-12-22 21:17
vibeta提示：If psfc is not a constant, then its dimension sizes must be equal to all but the rightmo ...

请问解决了吗，我也遇到了这个问题~

阿吉木柯

lt5201314 发表于 2022-9-19 15:11
请问解决了吗，我也遇到了这个问题~

看到了我就回答一下吧。因为你的x最右纬度不是level。把level这个纬度搞到最后就行了，例如x(time|:,lat|:,lon|;,level|:)

旺旺小小酥

卢小傻子 发表于 2020-5-1 08:38
可以用另一个函数，官网上也指出了  实现从地面开始积分吧
This matches an alternative approach
       ...

请问这个这个x是什么呀？