mvc_horiz_mod.F

路径

LMDZ.MARS\libf\dynphy_lonlat\phymars\mvc_horiz_mod.F

所属目录/模块

libf/dynphy_lonlat/phymars

文件定位

最大值复合（Maximum Value Composite）水平插值模块。将高分辨率地表数据（如 MOLA 地形、summit、base 等亚网格地形参数）映射到 GCM 目标网格，对每个目标格点取其覆盖范围内所有源格点的最大值。由 datareadnc 在处理亚网格地形变量时调用。

作者：Chao Wang。

定义的符号

依赖的模块

调用的关键例程

输入

输出

共享状态与副作用

无。纯计算例程，不读写任何 module 变量、common block 或文件。

核心逻辑

1. 初始化：iik/jjk 清零，outdata 初始化为 -999999（远低于火星地表最低点 ~8200 m）
2. 构建目标网格边界：
   • tmplon(1)=rlonv(1)，tmplon(2..imn+1) 为相邻经度中点，tmplon(imn+2)=rlonv(imn+1)
   • tmplat(1)=rlatu(1)，tmplat(2..jmn+1) 为相邻纬度中点，tmplat(jmn+2)=rlatu(jmn+1)
3. 源→目标映射：遍历每个源格点 (i,j)，在目标网格中找到其落入的单元格 (m,n)，记录到 iik(i,j)=m、jjk(i,j)=n；使用 goto 1000 找到即跳出内层循环
4. 最大值复合：遍历目标网格每个单元格 (m,n)，收集所有映射到该单元的源格点，取 max(outdata(m,n), zdata(i+imo*(j-1)))

伪代码

Initialize outdata = -999999 (sentinel)
Build template grid boundaries (tmplon, tmplat) at midpoints

FOR each source point (i,j):
    FOR each target cell (m,n):
        IF rlonuo(i) in [tmplon(m), tmplon(m+1))
           AND rlatvo(j) in [tmplat(n+1), tmplat(n)):
            Record iik(i,j)=m, jjk(i,j)=n
            GOTO next source point

FOR each target cell (m,n):
    FOR each source point (i,j):
        IF iik(i,j)==m AND jjk(i,j)==n:
            outdata(m,n) = max(outdata(m,n), zdata(i+imo*(j-1)))

参与的主题流程

写法特点

• 纯 Fortran 90 模块：无 USE/INCLUDE/CALL，完全自包含
• goto 1000 跳转：找到目标单元格后立即跳出，避免 O(imojmoimn*jmn) 全遍历
• -999999 哨兵值：注释说明 "lowest point on Mars ~8200m"，确保任何有效数据都能替换初始值
• 一维→二维索引：源数据 zdata(i+imo*(j-1)) 使用列优先一维展开，与 NetCDF 读取格式对应
• 四重循环复杂度：映射步骤和复合步骤均为 O(N²·M²)，适用于本场景（源/目标网格均较小）

复现要点

• 输入经纬度单位为弧度（pi=2*ASIN(1.) 虽声明但未使用）
• rlonv/rlatu 维度为 (imn+1)/(jmn+1)，即目标网格边界坐标
• outdata 中保持 -999999 的格点表示无源数据落入（调用方 datareadnc 后续替换为 0）

待确认

• pi 变量已声明并赋值但未使用，可能为遗留代码
• 纬度方向 rlatvo(j) .lt. tmplat(n) .and. rlatvo(j) .ge. tmplat(n+1) 假设纬度递减（北→南），对递增排列不适用
• 对于跨 0° 经度的源格点，周期性处理由调用方负责，本例程不处理

相关页面

• datareadnc — 唯一调用方，在 k=5,6,7（hmons/summit/base）时使用 mvc_horiz
• common-dependencies — interp_horiz 的 COMMON 水平插值边界；与本文件最大值复合插值形成对照。
• dynphy_lonlat/index — 所属目录概览