calfis.F

路径

LMDZ.COMMON-6.3\LMDZ.COMMON\libf\dynphy_lonlat\calfis.F

calfis.F 定义串行动力调用物理的接口例程 calfis。它把动力网格上的压力、温度、风、tracer、位势和垂直质量通量整理到物理列数组，调用 call_physiq 获得物理 tendency，再把 tendency 映射回动力网格。

条件经过：COMMON 串行动力启用物理时经过；Mars 标准物理入口多在 phymars 侧，但本文件保留 COMMON lon-lat 串行动力-物理接口逻辑。

符号	类型	行号	作用
`calfis`	subroutine	6	串行动力-物理调用接口。

依赖	用途
`infotrac:nqtot,niadv,tname`	tracer 数量、物理 tracer 名称和 advection index 映射。
`control_mod:planet_type,nsplit_phys`	行星类型和物理步分裂次数。
`callphysiq_mod:call_physiq`	行星物理统一调用入口。
`cpdet_mod:tpot2t,t2tpot`	位温和温度互转，支持 Cp(T) 适配。
`comconst_mod`	`daysec/dtvr/dtphys/kappa/cpp/g/rad/pi` 等常量。
`gr_dyn_fi`、`gr_fi_dyn`	动力格点与物理列互转。

参数组	说明
动力状态	`pps/pp/ppk/pteta/pq/pphi/pphis` 等，从动力网格转为物理列。
风场	`pucov/pvcov` 转为物理侧 `zufi/zvfi`，极点用经纬方向分解。
物理 tendency 输出	`pdufi/pdvfi/pdhfi/pdqfi/pdpsfi`，单位按注释为每秒 tendency。
时间标志	`debut` 保存首次调用，`lafin` 表示最后一次 physics 调用。

首次调用检查 ngridmx，分配 zonal mean 工作数组和纬度索引 klat。
将动力表面压力、压力层、Exner/位温和 tracer 从 iip1/jjp1 网格压缩为 ngridmx 物理列。
调用 tpot2t 得到物理温度 ztfi。
调用 gr_dyn_fi 转换位势和地表位势；计算物理列中的风、相对涡度和垂直质量通量 flxwfi。
在 CPP_PHYS 下按 nsplit_phys 循环调用 call_physiq，每个 split 使用 dtphys/nsplit_phys 和分裂后的 jH_cur_split。
每个 split 后用物理 tendency 更新 zufi/zvfi/ztfi/zqfi，并累计 zdufic/zdvfic/zdtfic/zdqfic。
split 结束后把累计 tendency 除以 nsplit_phys。
调用 gr_fi_dyn 或显式循环把压力、热力、tracer、u/v tendency 写回动力网格；tracer 写回使用 iiq=niadv(iq)。