Plotting data and scripts for "The Altitudinal Profile of the Equatorial Ionization Anomaly Response to the intense X9.3 Solar Flare on 6 September 2017"

# Reviewer-Response Figure Reproduction Package This folder contains the simplified and desensitized plotting package prepared for the manuscript revision of the 6 September 2017 X9.3 solar flare study. The reorganization follows the reviewer request to make the figure generation procedure clearer and easier for readers to run, without requiring access to the full original TIEGCM, Swarm, GOES, or intermediate diagnostic files. ## Purpose The reviewer specifically asked for clearer support of the figure interpretation and for a practical description of how the TIEGCM-based results can be used by readers. To address this, the original complex scripts were separated into two stages: 1. The full physical calculation workflow exports only the core plotting variables into `Fig1_Data.mat` to `Fig7_Data.mat`.2. The simplified scripts in this folder read those exported `.mat` files and reproduce Figures 1-7 directly. This keeps the physical results unchanged while removing dependence on large raw model files, satellite CDF files, NetCDF products, and local intermediate variables. ## Files | File | Role || --- | --- || `Fig1_Data.mat` - `Fig7_Data.mat` | Desensitized exported plotting data for Figures 1-7. || `figure1.m` - `figure7.m` | Minimal MATLAB plotting scripts. Each script generates one figure. || `figutil.m` | Shared helper for loading exported data, compatible field lookup, matrix orientation, and the common diverging colormap. || `run_all_figures.m` | One-command entry point to regenerate all figures. || `Figure1.png` - `Figure7.png` | Output images generated by the simplified scripts. | ## How To Run Open MATLAB, set the current folder to this directory, and run: ```matlabrun_all_figures``` This regenerates: ```textFigure1.pngFigure2.pngFigure3.pngFigure4.pngFigure5.pngFigure6.pngFigure7.png``` Each figure can also be generated independently, for example: ```matlabfigure5``` All scripts save figures using the same format: ```matlabprint(fig, 'FigureX', '-dpng', '-r600');``` ## Figure Mapping | Figure | Scientific Content | Reviewer-Response Role || --- | --- | --- || Figure 1 | IMF, Dst, GOES soft and hard X-ray fluxes | Clarifies the geomagnetic quiet condition and separates soft/hard X-ray flare signatures. || Figure 2 | Swarm A/B/C Ne observations versus TIEGCM simulations | Documents model-observation consistency using the exported comparison data. || Figure 3 | Global altitude-dependent Ne perturbations at 300, 400, and 500 km | Supports the revised discussion of spatial and vertical EIA response. || Figure 4 | Altitude-MLat Ne perturbation profile at GLon = 0 deg | Highlights bottomside enhancement and topside depletion in the revised interpretation. || Figure 5 | O+ continuity-term diagnostics | Provides the revised mechanism separation for photochemical production, ambipolar diffusion, winds, and E x B transport. || Figure 6 | Wind-driven vertical drift components | Supports the revised discussion of zonal and meridional wind contributions. Panel (a) uses a ±10 m/s scale; panels (b) and (c) use ±40 m/s to reflect the dominant meridional wind forcing. || Figure 7 | Perturbed vertical E x B drift and horizontal drift vectors | Clarifies the reviewer concern about downward drift versus horizontal arrow directions. | ## Data Desensitization The plotting package is physically desensitized: the scripts do not read the original TIEGCM NetCDF outputs, GOES NetCDF files, Swarm CDF files, magnetic-coordinate files, or diagnostic intermediate structures. Instead, each `.mat` file stores the final arrays needed for plotting. This makes the reproduction step lightweight while preserving the visual appearance and quantitative values of the original figures. The exported data include, for example: - Figure 1: IMF parameters and GOES soft/hard X-ray fluxes.- Figure 2: Swarm/TIEGCM sampled Ne tracks and correlation coefficients.- Figure 3: longitude, latitude, magnetic latitude, SZA or equivalent solar-terminator information, and altitude-layer Ne perturbations.- Figure 4: latitude-altitude slices of absolute and relative Ne perturbations.- Figure 5: four O+ diagnostic-term matrices and the latitude/height axes.- Figure 6: wind-driven vertical drift profiles.- Figure 7: gridded vertical E x B drift and horizontal drift vectors. ## Notes For Reviewers And Readers - The simplified scripts are intended for figure reproduction and inspection, not for rerunning the full TIEGCM simulation.- The full physical workflow remains the source of the exported `.mat` data; this package provides the reader-facing plotting layer requested during revision.- The visual settings, including panel positions, font sizes, color limits, contour lines, vector overlays, and 600 dpi PNG export, were kept consistent with the manuscript figures.- Figure 1 uses explicit logarithmic tick labels for the X-ray panels to avoid MATLAB automatic exponent formatting issues.- Key quantitative values reported in the manuscript (peak ΔNe, relative enhancement percentages, wind drift magnitudes, and thermospheric stability statistics) are reproducible directly from the exported .mat files using the diagnostic print statements included in each figure script.# 审稿人回复专用图片复现套件 本文件夹包含为2017年9月6日X9.3级太阳耀斑研究的稿件修订工作所准备的简化脱敏绘图套件。本次重构遵循审稿人要求，旨在使图片生成流程更清晰、便于读者复现，且无需获取完整原始的热层-电离层-电动力学通用环流模型（TIEGCM）、Swarm卫星群（Swarm）、地球静止运行环境卫星（GOES）数据或中间诊断文件。 ## 设计目的 审稿人特别要求更清晰地支持图片解读，并提供基于TIEGCM的结果如何被读者使用的实操说明。为此，我们将原复杂脚本拆分为两个阶段： 1. 完整物理计算流程仅将核心绘图变量导出至`Fig1_Data.mat`至`Fig7_Data.mat`。 2. 本文件夹内的简化脚本可读取这些导出的.mat文件，直接复现图1至图7。 此举在保留物理结果不变的同时，移除了对大型原始模型文件、卫星公共数据格式（CDF）文件、网络通用数据格式（NetCDF）产品及本地中间变量的依赖。 ## 文件说明 | 文件 | 作用 | | --- | --- | | `Fig1_Data.mat` 至 `Fig7_Data.mat` | 用于图1至图7的脱敏导出绘图数据。 | | `figure1.m` 至 `figure7.m` | 极简MATLAB绘图脚本，每个脚本生成对应单张图片。 | | `figutil.m` | 共享辅助工具，用于加载导出数据、兼容字段查找、矩阵方向调整及通用发散色图。 | | `run_all_figures.m` | 一键生成所有图片的入口脚本。 | | `Figure1.png` 至 `Figure7.png` | 简化脚本生成的输出图片。 | ## 运行方法 打开MATLAB，将当前工作目录设置为本文件夹，执行以下命令： matlab run_all_figures 该命令将复现以下图片： text Figure1.png Figure2.png Figure3.png Figure4.png Figure5.png Figure6.png Figure7.png 也可独立生成单张图片，例如： matlab figure5 所有脚本均采用统一格式保存图片： matlab print(fig, 'FigureX', '-dpng', '-r600'); ## 图片对应关系 | 图片编号 | 科学内容 | 审稿人回复作用 | | --- | --- | --- | | 图1 | 行星际磁场（IMF）、地磁Dst指数（Dst）、GOES软X射线和硬X射线通量 | 阐明地磁静日条件，并区分软/硬X射线耀斑特征。 | | 图2 | Swarm A/B/C卫星的电子数密度（Ne）观测值与TIEGCM模拟结果对比 | 通过导出的对比数据验证模型与观测的一致性。 | | 图3 | 300km、400km、500km高度处的全球高度依赖型Ne扰动 | 支持关于空间与赤道电离层异常（Equatorial Ionization Anomaly，EIA）响应的修订讨论。 | | 图4 | 地理经度0°处的高度-磁纬度Ne扰动剖面 | 在修订后的解释中突出底部侧增强与顶部侧耗竭现象。 | | 图5 | 氧离子（O+）连续项诊断结果 | 为光化学产生、双极扩散、风场及电场×磁场（E×B）漂移提供修订后的机制分离说明。 | | 图6 | 风驱动的垂直漂移分量 | 支持关于纬向与经向风贡献的修订讨论。子图(a)采用±10 m/s的色标；子图(b)与(c)采用±40 m/s的色标，以体现主导的经向风强迫效应。 | | 图7 | 扰动垂直E×B漂移与水平漂移矢量 | 解答审稿人关于向下漂移与水平箭头方向的疑问。 | ## 数据脱敏处理 本绘图套件已完成物理层面的脱敏处理：脚本不会读取原始TIEGCM NetCDF输出文件、GOES NetCDF文件、Swarm卫星CDF文件、磁坐标文件或诊断中间结构。取而代之的是，每个.mat文件存储了绘图所需的最终数组。此举使复现步骤轻量化，同时保留了原始图片的视觉效果与定量数值。 导出的数据示例包括： - 图1：IMF参数与GOES软/硬X射线通量。 - 图2：Swarm卫星与TIEGCM采样的Ne轨迹及相关系数。 - 图3：经度、纬度、磁纬度、太阳天顶角（Solar Zenith Angle，SZA）或等效的日界线信息，以及各高度层的Ne扰动。 - 图4：绝对和相对Ne扰动的纬度-高度切片。 - 图5：四个O+诊断项矩阵及纬度/高度坐标轴。 - 图6：风驱动的垂直漂移剖面。 - 图7：网格化的垂直E×B漂移与水平漂移矢量。 ## 致审稿人与读者的说明 - 本简化脚本仅用于图片复现与检视，不可用于重新运行完整的TIEGCM模拟。 - 完整物理流程仍是导出.mat数据的唯一来源；本套件仅提供修订阶段收到的面向读者的绘图层。 - 视觉设置（包括面板布局、字体大小、色限、等高线、矢量叠加及600dpi PNG导出）均与稿件图片保持一致。 - 图1的X射线面板采用显式对数刻度标签，以避免MATLAB自动指数格式带来的显示问题。 - 稿件中报告的关键定量数值（峰值ΔNe、相对增强百分比、风漂移幅值及热层稳定性统计量）均可通过每张图片脚本中包含的诊断打印语句，直接从导出的.mat文件中复现。