WZZHHH的博客

上图中红线画出的地方就是两个通用变量名称的改变：level 变为 pressure_leveltime 变为 valid_time。

小白就用简单的把，别纠结那么多，哪个代码容易敲还简单就用哪个，对结果的影响也不是很大。

使用水汽压和温度计算 Tm和 表示第 i 层大气的水汽压和温度, 是第 i 层的厚度。这种方法使用探空站数据计算得到的 Tm 一般为真值与其他数据进行比较。

GNSS ERA5、ERA5 PWV和探空站PWV(RS PWV)进行相互比较。ZTD我用gamit计算得到，ZHD需要气压用Saastamoinen模型计算，该模型需要压力和站位信息。理想情况下，根据积分法，可以通过探空站的大气剖面得到Tm或者ERA5数据积分。因此，压力、温度和Tm是GNSS PWV检索的关键参数。本文使用ERA5数据集通过水平和垂直方向的补偿得到GNSS站点垂直方向的ERA5数据集，得到GNSS站点处的温度、气压和Tm，并积分得到该位置的ERA5 PWV。代码结果如图所示ERA5数据下

我们在获取GNSS站点位置的PWV有很多方法，其中一个是利用ERA5数据进行积分得到。由于ERA5格网位置和GNSS不并址，需要进行水平方向和垂直方向的补偿的补偿获得。这里讲的是垂直方向补偿，水平方向补偿已经写过反距离加权水平内插，附matlab代码(ERA5和GNSS站点不并址的处理方法之水平补偿)-优快云博客，这里介绍垂直方向的补偿。水平方向补偿方法ERA5多层温度、位势、相对湿度的下载方法在这面其中相对湿度取最接近层数据即可，温度和气压用下面方法使用内插估计站点参数即可。假定温度递减率为-6.5 K

我在学习过程，内插方法为反距离加权水平内插，分享我的方法和公式，以及matlab代码。

0. 码字不易，点赞加关注（公众号：WZZHHH，部分资料在公众号可以下载），使用请注明出处（根据我的研究方向，我会不断更新代码）。1.计算PWV的方法一般采用有三种，GNSS ERA5、ERA5 PWV和探空站PWV(RS PWV)进行相互比较。2.本博文使用ERA5数据集通过水平和垂直方向的补偿得到GNSS站点垂直方向的ERA5数据集，计算的GNSS站点上方的ERA5 PWV。3.ERA5的下载方法在前面的博文，自行学习下载4.下面是我的部分代码，具体获取请关注咸鱼：WZZHHH123

6.在下载的页面，区域大以及下载数据较多需要等待较长时间。记得检查下载数据的参数和年月日等是否正确。等待完成后，点击下载就行。4.选择自己的所需的时间段，比如2020年1月1日的24个小时的数据。3.选择合适的参数，比如10m东风、10m北风、和2m温度。5.选择研究区域的经纬度区间，下载格式为NC，然后点击提交。2.选择下载数据，和再分析数据。1.首先进入ERA5数据集。

温度气压温度。