【求助】关于垂直积分

Nico_凌风

对多层物理量进行垂直积分，浏览许多相关帖子没有得到想要的答案后，有几个问题请各位大神指点：
1、GrADS的vint函数，不清楚其“质量加权”是否为普通的垂直积分，但用vint和普通积分得出的结果相差很大；
2、vint的第一个和第三个参数是否为另下载的表层气压（第一个参数）、且顶层自设的固定值（第三个参数）在所积分物理量ctl中存在即可？
3、一般的积分算法理论上是否可以无限的分割得到中间值最后相加？
4、接着上面的问题，对于1000-300hPa的资料，NCEP仅有8层，而ERA5有20层，我分别垂直积分后，ERA5的值明显比NCEP的大，这样是否正常？是否需要层次数目相同？

非常感谢！

苏拉苏拉

Nico_凌风 发表于 2020-7-24 14:42
也就是说，表层气压数据pres的单位要为hPa(mb)吗？，代入vint中后×其自带的100变成Pa?

是的，vint里带入的气压需要hPa/mb为单位。这是介绍，你看看http://cola.gmu.edu/grads/gadoc/gadocindex.html

苏拉苏拉

3、我认为是的
4、一方面资料本身有差异，另一方面函数计算也会有误差，我感觉层数越密积分应该越准确。你说ERA5的值明显比NCEP大，大多少，如果不是有量级程度的差异，应该也可以接受吧，可以再用第三个数据去验证一下谁准确些。或者你可以把ERA5里对应的那8层挑出来，单独积分一次，结果比较一下，看是不是层数多少造成的差异。

苏拉苏拉

1、vint和普通积分不同，大概差一个量级。
2、是的、以及是的。

苏拉苏拉

1、补充。官网介绍：The calculation is a sum of the mass-weighted layers: f/g * sum(expr * Delta(level))。The scale factors are f=100 and g=9.8  

Nico_凌风

苏拉苏拉 发表于 2020-7-23 18:05
1、补充。官网介绍：The calculation is a sum of the mass-weighted layers: f/g * sum(expr * Delta(leve ...

也就是说，表层气压数据pres的单位要为hPa(mb)吗？，代入vint中后×其自带的100变成Pa?