leapfrog 动力计算和 tracer 输送段

输入范围

LMDZ.COMMON-6.3\LMDZ.COMMON\libf\dyn3d\leapfrog.F
LMDZ.MARS\libf\dynphy_lonlat\phymars\leapfrog_nogcm.F

Mars 运行参与度：条件经过。COMMON 完整 GCM 入口会执行本页的动力/输送段；MARS nogcm 变体保留接口但把动力 tendency 清零，不调用 caldyn/caladvtrac/integrd。

完整 GCM 调用顺序

每个 label 2 步中，完整 leapfrog.F 的动力段顺序是：

1. tpot2t 把 teta 和 pk 转成温度 temp。
2. 构造 tsurpk=cpp*temp/pk。
3. geopot(ip1jmp1,tsurpk,pk,pks,phis,phi) 计算 geopotential。
4. caldyn(...) 计算 du/dv/dteta/dp/w/pbaru/pbarv。
5. 若 planet_type=="generic" 且 ok_guide，把 uforc nudging 加到 du。
6. 若 forward .OR. leapf，调用 caladvtrac(q,pbaru,pbarv,p,masse,dq,teta,flxw,pk)。
7. 若 offline 且启用 IOIPSL，调用 fluxstokenc 写 mass flux。
8. 调用 integrd(nqtot,...,dv,du,dteta,dq,dp,...,q,ps,masse,phis) 写回状态。
9. Titan tidal 分支可额外加 tidal_forces。

tendency 输入输出

integrd 接收 nqtot，因此 tracer 维度由 infotrac contract 控制。

NODYN 和 nogcm

完整 leapfrog.F 有 #ifdef NODYN 分支：将 dv/du/dteta/dq/dp/flxw 置零，绕过 caldyn/caladvtrac/integrd 前的动力计算。

MARS phymars/leapfrog_nogcm.F 则直接把：

dv=0
du=0
dteta=0
dq=0

并跳过完整动力段和 integrd。它后续通过 physics 后的 pressure/CO2/temperature mixing 来替代部分再分布效果。

与已有 helper 页的关系

• caldyn 串并行对照：解释 du/dv/dteta/dp/w/pbaru/pbarv 的动力来源。
• advtrac 串并行对照：解释 iadv 如何决定 tracer 输送路径。
• integrd 串并行对照：解释 tendency 写回和非法压力检查。
• fluxstokenc 串并行对照：解释 offline flux storage。

复现要点

• caladvtrac 使用步初状态和 caldyn 给出的 mass flux；不要把 physics 后的 masse 回填到这一段。
• 完整 GCM 中 caladvtrac 只在 forward 或 leapf 步调用，不在 backward Matsuno 步调用。
• nogcm 页面中的 q 仍会被 physics 和 CO2 mixing 改变，但不是由 caladvtrac 改变。

待确认

• NODYN 预处理构建与 nogcm 运行的实际使用差异，需要 build/arch 页面确认。