垂直积分出错

dengxiaohua_121

各位前辈，我使用
vint_qu=vibeta(p_new,qu_record,linlog,ps,pbot,ptop)/g*100  
提示的错误还是：warning:vibeta: there must be at least three levels with data above the surface
fatal:vibeta: ier = -999
可我一步一步检查了一下，觉得应该没有问题啊，会的前辈能帮我看一下是什么原因吗？
Variable: p_new
Type: float
Total Size: 32 bytes
            8 values
Number of Dimensions: 1
Dimensions and sizes: [level | 8]
Coordinates:
            level: [1000..100]
Number Of Attributes: 3
  long_name : pressure_level
  units : millibars
  _FillValue : 9.96921e+36
(0) 1000
(1) 975


Variable: qu_record
Type: float
Total Size: 17911552 bytes
            4477888 values
Number of Dimensions: 4
Dimensions and sizes: [time | 31] x [lat | 122] x [lon | 148] x [leve
l | 8]
Coordinates:
            time: [2014071206..2014071918]
            lat: [50.25..-40.5]
            lon: [39.75..150]
            level: [1000..100]
Number Of Attributes: 1
  _FillValue : 9.96921e+36


Variable: ps
Type: float
Total Size: 2238944 bytes
            559736 values
Number of Dimensions: 3
Dimensions and sizes: [31] x [122] x [148]
Coordinates:
Number Of Attributes: 1
  _FillValue : -32767
(0,0,0) 1014.408
(0,0,1) 1014.991
(0,0,2) 1015.525
(0,0,3) 1016.011
(0,0,4) 1016.441

Variable: pbot
Type: integer
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 1000


Variable: ptop
Type: integer
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 100

talkd

第一，ps最好给属性，第二，你计算的地方高度多少？ps以上有3层吗？

dengxiaohua_121

谢谢楼上的回复，ps里头的值是海平面气压值，计算的区域为-40至50，东经40至150，是因为含了青藏高原吗，所以ps以上有的区域不足3层？

longlivehj

dengxiaohua_121 发表于 2015-9-18 11:13
谢谢楼上的回复，ps里头的值是海平面气压值，计算的区域为-40至50，东经40至150，是因为含了青藏高原吗，所 ...

正解

dengxiaohua_121

谢谢longlivehj。我试试，论坛里头好像看到有怎么往上加几层的例子。

气象妹子

我算的区域是30S--30N，60E-160E，，没有缺测值，，为什么还是会提示这个错误………………{:eb303:}{:eb303:}

日不不不落

楼主 你好 你的问题解决了吗  可否告知一下解决方法呢

极地冷高压

longlivehj 发表于 2015-9-18 12:55
正解

真的是这样吗

servicewangyuqi

可是我的数据没有ps这一项要怎么办，就是海平面的气压值，是模拟结果

zys1022

我也正有这样的困惑，pbtop=1000,那1000hpa以下的高度层就表示不出来了？