从零开始掌握CESM2集合预报：pertlim参数实战指南

刚接触气候模式的研究者常会困惑：为什么同样的模式设置，每次运行结果都有细微差异？这背后隐藏着大气系统的混沌本质。1963年气象学家洛伦兹用"蝴蝶效应"形象说明，即使初始条件存在极其微小的差异，经过足够长时间也会导致完全不同的天气状态。集合预报正是应对这种不确定性的科学方法——不是做一次模拟，而是同时进行多组略有差异的模拟，通过统计规律提高预测可靠性。

作为地球系统模式中的佼佼者，CESM2（Community Earth System Model version 2）提供了便捷的集合实验功能。本文将聚焦pertlim这个关键参数，带您完成首个集合预报实验。不同于单纯的技术文档，我们会用烹饪来类比：就像在蛋糕配方中微调糖量观察成品变化，pertlim就是对大气初始温度的"调味料"，1.d-14这样的数值就是我们的"调料计量单位"。

1. 集合预报的核心概念与价值

1.1 为什么需要多组模拟？

想象预测某地未来两周的天气。如果只运行一次模式，就像只抛一次硬币就断定正反面概率。集合预报相当于同时抛几十枚硬币，通过正反面比例计算可靠概率。具体来说：

• 误差过滤：各成员结果的共性部分更可能反映真实趋势，随机误差会被平均弱化
• 可靠性评估：成员间差异程度直接反映预测可信度（差异越小越可靠）
• 概率预报：可以计算某天气现象出现的百分比概率

提示：典型的集合规模在20-50个成员之间，具体数量取决于计算资源和研究目的

1.2 pertlim的物理意义解析

pertlim参数的单位看起来令人费解（如1.d-14），其实这是科学计数法的Fortran表示方式：

fortran复制1.d-14 = 1 × 10⁻¹⁴

这个量级相当于给初始大气温度场施加约0.00000000000001K的扰动——比蝴蝶扇动翅膀引起的空气温度变化还要小多个数量级。但经过模式积分，这种微观扰动会被非线性过程放大，最终形成可观测的差异。

2. 实验环境准备

2.1 基础环境检查

确保已经完成以下前置步骤：

1. CESM2成功安装并通过基础测试案例
2. 拥有至少50GB的临时存储空间（用于多个成员输出）
3. 掌握基本的case创建和提交命令

验证环境是否就绪：

bash复制# 检查CESM根目录设置
echo $CESM_ROOT_DIR

# 测试简单案例运行
cd $CESM_ROOT_DIR/cime/scripts
./create_newcase --case test_case --compset F2000climo --res f19_g17
cd test_case
./case.setup
./case.build
./case.submit

2.2 创建集合实验案例

建议使用以下组合配置，平衡计算成本与教学效果：

创建命令示例：

bash复制./create_newcase --case ensemble_base --compset F2000climo --res f19_g17

3. pertlim参数实战配置

3.1 单成员参数设置

进入案例目录后，修改user_nl_cam文件：

1. 用文本编辑器打开配置文件
2. 在文件末尾添加扰动参数
3. 保存并关闭文件

典型配置示例：

fortran复制! 初始场扰动参数设置
&cam_initfiles_nl
   pertlim = 1.d-14
/

注意：Fortran的namelist格式要求严格，&符号后不能有空格，/必须单独成行

3.2 批量生成多成员脚本

手动修改每个成员效率低下，推荐使用脚本自动化：

bash复制#!/bin/bash

for i in {1..5}  # 生成5个成员
do
   # 创建成员目录
   cp -r ensemble_base ensemble_$i
   
   # 进入目录修改参数
   cd ensemble_$i
   echo "&cam_initfiles_nl" > user_nl_cam
   echo "   pertlim = $i.d-14" >> user_nl_cam
   echo "/" >> user_nl_cam
   
   # 提交作业
   ./case.submit
   cd ..
done

将上述脚本保存为generate_ensemble.sh后执行：

bash复制chmod +x generate_ensemble.sh
./generate_ensemble.sh

4. 结果分析与常见问题

4.1 基础验证方法

运行结束后，检查各成员是否产生预期差异：

bash复制# 查看表面温度输出差异
ncdiff ensemble_1/run/ensemble_1.cam.h0.0001-01.nc \
       ensemble_2/run/ensemble_2.cam.h0.0001-01.nc \
       diff_1_2.nc

预期应该看到：

• 初始时间步差异极小（≈1e-14量级）
• 随时间推移差异逐渐扩大
• 热带地区差异发展更快（混沌特性更强）

4.2 典型报错解决方案

问题1：ERROR: pertlim not recognized

• 检查拼写是否正确（应为pertlim不是pert_limit）
• 确认参数放在&cam_initfiles_nl而非其他namelist块

问题2：成员结果完全相同

• 确认各成员的user_nl_cam中pertlim值确实不同
• 检查是否意外复制了相同的初始条件文件

问题3：运算中途崩溃

• 尝试减小pertlim值（如改为1.d-13）
• 检查计算节点内存是否充足

4.3 高级技巧：扰动量级选择

不同研究目的适用的扰动强度：

在最近一次海温异常研究中，我发现当pertlim大于1.d-10时，部分成员会出现非物理震荡。这提醒我们参数设置需要结合具体compset测试调整。