OLCI script AMT paper

Analysis code for "Boundary layer turbulence characterization from OLCI-derived integrated water vapor structure functions: a global multi-site analysis" (Kermarrec, 2026, AMT).

**Contents:** `olci_amt_analysis.m` (main script, reproduces all figures and tables), `seasonal_count.m` (monthly/seasonal product counts).

**Requirements:** MATLAB R2022b+, Statistics and Machine Learning Toolbox, Mapping Toolbox (optional for world map).

**Input data:** (1) OLCI SF results as `.mat` files with fields `heff_ref` (L₀), `a1_ref` (α₁), `a2_ref` (α₂), `C_ref` (σ²), `sf_ref` (60-bin SF), `centers`, `heff_detection_rate`, `iwv_mean`, produced by `olci_batch_v4b_generic.m`. (2) ERA5 `.nc` files: hourly BLH (`era5_{site}_hourly_*.nc`) and sub-daily wind/TCWV (`era5_{site}_single_*.nc`) from the Copernicus Climate Data Store.

**Usage:** Set `MAT_DIR`, `ERA5_DIR`, `ERA5_HOURLY_DIR` at the top of `olci_amt_analysis.m`, then run. Figures are saved to `./figures/`. Run `seasonal_count.m` afterwards for product distribution tables.

**Sites:** Highveld (n=42), India (n=115, 2016–2024), Pampas (n=58), US Plains (n=177, 2016–2024), Yangtze (n=31).

**Key parameters:** `RD_MIN=0.85`, `HEFF_MIN/MAX=350/3000 m`, `MAD_THRESH=3.5`, `SMOOTH_WIN=5` months.

**Contact:** Gaël Kermarrec, [e-mail], Institute of Meteorology and Climatology, Leibniz Universität Hannover.

本代码对应论文《基于海洋与陆地颜色仪器（OLCI）反演的积分水汽结构函数的边界层湍流特征表征：全球多站点分析》（Kermarrec, 2026, AMT）。

**内容：** `olci_amt_analysis.m`（主脚本，可复现所有图表与表格）、`seasonal_count.m`（用于统计月度/季度产品数量）。

**运行要求：** MATLAB R2022b及以上版本，需安装统计与机器学习工具箱（Statistics and Machine Learning Toolbox），制图工具箱（Mapping Toolbox）为可选依赖，用于绘制全球地图。

**输入数据：** (1) OLCI结构函数结果，存储为`.mat`格式文件，包含字段`heff_ref`（特征高度L₀）、`a1_ref`（α₁）、`a2_ref`（α₂）、`C_ref`（方差σ²）、`sf_ref`（60-bin结构函数）、`centers`、`heff_detection_rate`、`iwv_mean`，上述结果由`olci_batch_v4b_generic.m`生成。(2) ERA5 `.nc`格式文件：来自哥白尼气候变化数据服务中心的逐小时边界层高度（Boundary Layer Height, BLH）文件`era5_{site}_hourly_*.nc`，以及亚日尺度风和总柱水汽（Total Column Water Vapor, TCWV）逐次观测文件`era5_{site}_single_*.nc`。

**使用方法：** 在`olci_amt_analysis.m`文件顶部设置`MAT_DIR`、`ERA5_DIR`、`ERA5_HOURLY_DIR`三个路径参数后运行脚本，生成的图表将保存至`./figures/`目录。后续可运行`seasonal_count.m`以生成产品分布统计表。

**研究站点：** 海菲尔德（Highveld, n=42）、印度（India, n=115, 2016–2024）、潘帕斯草原（Pampas, n=58）、美国大平原（US Plains, n=177, 2016–2024）、长江流域（Yangtze, n=31）。

**关键参数设置：** `RD_MIN=0.85`、`HEFF_MIN/MAX=350/3000 m`、`MAD_THRESH=3.5`、`SMOOTH_WIN=5`个月。

**联系方式：** Gaël Kermarrec，邮箱：[e-mail]，所属机构：汉诺威莱布尼茨大学气象与气候学研究所（Institute of Meteorology and Climatology, Leibniz Universität Hannover）。