观潮 · 模式预报与气象 AI 观察
scale_region.py 把网格从 46—12 公里加密到 15—4 公里:静态场看着都对,飓风却在几小时内消散
> 一位用户把 Helene 个例的区域网格用 scale_region.py 从 46—12 公里加密到 15—4 公里,除按分辨率调整积分步长外,namelist、streams、graph.info 全部与原网格运行保持一致。结果原网格上的飓风在第 1、5、10 小时的海平面气压与 850 hPa 风场都正常,加密网格上飓风迅速消散,日志没有任何报错、积分也没有提前终止。NCAR 收集了静态场与 namelist 后,回复延迟,并明确表示「还没有找到解决你这个问题的方法」。用户排除了 Grell-Freitas 尺度敏感项这一嫌疑:关掉对流之后,飓风依然迅速消散。
症状:换了分辨率,气旋没了
操作流程是这样的:先用 grid-rotate 与 MPAS-Limited-Area-master 在 46km—12km 的变分辨率网格上做区域化,把区域中心放在墨西哥湾,得到 Helene.grid.nc。然后跑两组积分:
网格来源 · 分辨率 · 结果
第一组 · Helene.grid.nc · 46—12 公里 · 初始时刻、第 5 小时、第 10 小时的 MSLP 与 850 hPa 风场看起来正常
第二组 · Helene.rgrid.nc · 15—4 公里 · 飓风迅速消散
第二组的网格由加密脚本生成:
```
./scale_region.py Helene.grid.nc Helene.rgrid.nc 3 22 -88
```
除按分辨率调整 config_dt 之外,两组运行的 namelist 与 streams 完全一致,graph.info 也相同。用户检查后认为网格本身没问题,grid、init、static 文件看起来都是正常的,问题只在积分阶段出现——更麻烦的是日志没有任何报错,模拟也不提前终止,看起来像一次「正常跑完」的失败。
第一条线索:对流有没有打开
NCAR 方面的第一反应是从 namelist 里找问题,注意到对流被关掉了:
```
config_physics_suite = 'convection_permitting'
```
对方提出两个问题:为什么不想开对流?46—12 公里那次是否也关掉了?
用户的答复澄清了这一点:两次运行都用默认的 convection_permitting 套件;此外也用 mesoscale_reference 套件跑过,结果在性质上相似。至于 namelist 里出现的关对流配置,来自后续的计算——他特意这样设置,是为了检验 Grell-Freitas 的尺度敏感特性是不是问题根源,因为气旋中心及其附近区域表现出明显异常。
这条线索随即被否掉了:关掉对流的那组试验里,飓风同样迅速消散。
不是个案
用户补充,另一个西太平洋的热带气旋个例出现了类似问题。这让问题从「Helene 这一个例的配置失误」扩大到「局部加密网格上 TC 强度的共性现象」。
走到一半的官方跟进
NCAR 随后要求用户上传一批文件,以便复现排查:
编号 · 要求上传的内容
3 · namelist.init_atmosphere 与 streams.init_atmosphere(静态、初始与边界条件)
4 · namelist.atmosphere 与 streams.atmosphere(模式积分用)
6 · Helene.static.nc 与 Helene.rstatic.nc(对应 46—12 公里与 15—4 公里两套网格的静态资料)
用户按要求上传到对方的 Box 存储并告知。此后他又两次追问进展,NCAR 回复解释了延迟原因——假期与刚刚结束的一场 tutorial,并明确表示尚未找到解决办法,一旦弄清楚问题出在哪会再来答复。帖子最后一次更新仍是用户在问有没有后续进展。
这条帖子的边界在哪
- 没有归因,也没有解法。 官方到最后只完成了资料收集,明确说没找到问题所在。「关键在各通胀变量的季节性分量」这类同步的结论在这里不存在。
- "I don't believe there are issues here" 是用户自己的目视判断。 他检查 grid、init、static 后认为没问题,但没有给出量化验证。
- 对流这一项是被排除的假设,不是结论。 关对流组同样消散,只能说明它不是这次差异的来源。
- 两张对比图的形态无法在这条帖子里核验。 用户提到的 MSLP 与 850 hPa 风场图是附件形式,正文中没有给出数值化的强度对比。
一句话判断
如果要在 MPAS 上用 scale_region.py 做 TC 个例的局地加密,这条帖子值得先读:它给出了可复现的操作序列、明确的输入对比(除 config_dt 外全同),以及一条已经排除掉的错误假设;但截至目前,NCAR 尚未给出成因与修复方案。