对IGRAv2进行质量控制，得到IGRA-PWV的Matlab代码

1. Matlab 获取代码回复关键词：

怀俄明探空站数据解算PWV和Tm：怀俄明探空站数据解算PWV和Tm

怀俄明多线程下载探空站数据（包括检查和下载遗漏数据的代码）：怀俄明多线程下载

对IGRAv2进行质量控制得到PWV和Tm的 matlab 代码：IGRAv2进行质量控制得到PWV和Tm

算 IGRAv2 探空站的 Tm 和 PWV：IGRAv2计算Tm和PWV

提取探空站 IGRAv2 全部数据：

ERA5 解算合集（温度、气压、PWV、水汽压和 Tm）代码获取：ERA5合集

2. IGRA气压层对PWV的影响

  在论文（DOI: 10.1175/2007JCLI1944.1）中利用IGRA数据计算出的全球和年度平均温度和湿度曲线，分析了对流层顶以上、300、500 和 700 hPa 缺失数据对PWV的影响程度；以及对地表、1000 和 850 hPa 缺失数据对PWV的影响。

   分析表明，PW 对地表、1000和850 hPa 的缺失数据很敏感；300 hPa 以上的缺失数据只会给 PW 带来很少的干偏差。

3.IGRAv2数据质量控制

 根据上面的结论进行质量控制。

 1.顶层至少达到300 hPa

2.地面有数据可用（IGRAv2第一层即是地面数据，满足温度、湿度、气压有数据即可）

3.在低于(高于)1000hPa的台站，地面以上至少有五(四个)个标准压力水平

4.连续记录温度或湿度之间的压力间差值大于200 hPa的剔除。

在武大张卫星的论文（DOI: 10.1175/JCLI-D-16-0591.1）里面有介绍：

4. 进行质量控制得到IGRA-PWV的Matlab代码

如下图代码所示，压缩包的内容。

5. 下面是部分代码展示：

% 读取IGRA的文件
% 基于matlab2022写的
% IGRA下载地址https://www1.ncdc.noaa.gov/pub/data/igra/derived/derived-por/
% 里面有很多站点已经弃用，并且该网址包括建站以后的所有数据集,
% 因此我们需要提取对应年份的数据,然后计算得到探空站得PWV

%% -------------------------质量控制-----------------------------
% 我们实施的数据质量控制策略来源：DOI: 10.1175/2007JCLI1944.1。
% 1.顶层至少达到300 hPa，
% 2.地面有数据可用（第一层即是地面数据，满足温度、湿度、气压有数据即可）
% 3.在低于(高于)1000hPa的台站，地面以上至少有五(四个)个标准压力水平。
% 4.连续记录温度或湿度之间的压力间隙不能大于200 hPa。

%% ----------------------------------------------------------------
% 提取数据存到"v2data"，按照时间进行排列，最后存到mat文件中
% 每个站点的每年存一个mat文件，名称为：站点名字-年份.mat
% --------------数据的含义如下所示---------------

%     V2data.pwv       PWV(包括时间);
%     V2data.lon       经度;
%     V2data.lat       纬度;
%     V2data.Tm        加权平均温度(包括时间);
%     V2data.ID        站点名称；
%
% ------------------------------------------
clc,clear;