利用探空站数据（怀俄明和IGRA）和ERA5计算ZTD、ZHD和ZWD
1.有关 Matlab 获取代码关注咸鱼获取：
怀俄明探空站数据解算PWV和Tm
怀俄明多线程下载探空站数据（包括检查和下载遗漏数据的代码）
对IGRAv2进行质量控制得到PWV和Tm的 matlab 代码
算 IGRAv2 探空站的 Tm 和 PWV
提取探空站 IGRAv2 全部数据：
ERA5 解算合集（温度、气压、PWV、水汽压和 Tm）代码获取
2.计算ZWD
计算ZWD使用温度、水汽压、高程，利用下面公式进行计算。
水汽压的计算方法使用ECMWF里面的公式进行计算(ECMWF, 2007. IFS Documentation CY31R1 - Part II: Data Assimilation. https://doi.org/10.21957/m46uhsu4q.)，它是把温度分为三个区间：
3.计算ZHD
计算ZHD使用气压、纬度和海拔，利用Saastamoinen模型得到。
参考文献：Saastamoinen, J., 1972. Atmospheric Correction for the Troposphere and Stratosphere in Radio Ranging Satellites. In: Henriksen, S.W., Mancini, A., Chovitz, B.H. (Eds.), Geophysical Monograph Series. American Geophysical Union, Washington, D.C., pp. 247–251. https://doi.org/10.1029/GM015p0247.
4.计算ZTD
使用ZHD+ZWD=ZTD。
5.代码获取
关注上面的咸鱼，里面有三种方法计算，怀俄明探空站点、IGRAv2探空站点和ERA5格网计算