垂直离散构造

输入范围

dyn3d_common\disvert.F90
dyn3d_common\disvert_noterre.F

Mars 运行参与度：必经/条件经过。动力初始化必须构造垂直坐标；Mars 非地球路径常依赖 disvert_noterre.F 或相同接口。

例程定位

disvert 和 disvert_noterre 构造 ap/bp/presnivs 等垂直坐标数组。前者是较新的 Fortran 90 实现，读取多种 getin 参数；后者是非地球/传统路径，调用 sig_hybrid 构造 hybrid sigma 坐标。

对照表

核心流程

1. 读取垂直层数和坐标控制参数。
2. 根据采样方式构造层界面和层中点。
3. 计算 hybrid 坐标系数 ap/bp 与参考压力层。
4. 写入 comvert_mod 供 pressure、Exner、history 垂直轴使用。

复现要点

• 垂直坐标是后续 pressure、mass、temperature 转换的共同输入，不能只记录层数。
• vert_sampling 等配置会改变层厚分布；Mars 复现实验必须保存 run.def 中相关键。
• disvert_noterre 的非地球路径和 Mars 动力-物理共享状态约定见 vertical-coordinate-mars-contract。

相关页面

• dyn3d_common 总览
• 质量压力能量总览

已确认

• COMMON iniconst 在 planet_type!="earth" 时默认 disvert_type=2，即 Mars/Generic 默认调用 disvert_noterre；仍可由 run.def 的 disvert_type 覆盖。