WRF新手必看：Single Domain Case模拟全流程详解（附常见错误排查）
WRF单域模拟实战指南：从零开始到结果可视化
1. WRF单域模拟概述
天气研究与预报模型（WRF）作为当前最先进的中尺度气象模拟系统，已成为科研机构和业务部门的核心工具。单域案例（Single Domain Case）是WRF建模中最基础的配置方案，特别适合初学者掌握WRF工作流程和参数设置逻辑。
为什么选择单域模拟作为入门？
- 配置复杂度低：无需处理嵌套网格间的交互
- 计算资源需求小：适合在个人工作站或小型服务器上运行
- 调试方便：问题定位更直接，适合学习错误排查方法
- 核心功能完整：包含WRF所有关键流程环节
提示：虽然单域模拟相对简单，但它包含了WRF所有核心组件（WPS、real.exe、wrf.exe）的完整工作流程，是掌握更复杂嵌套模拟的必要基础。
2. 环境准备与数据获取
2.1 系统要求与软件依赖
WRF运行需要以下基础环境：
- Linux操作系统（推荐Ubuntu/CentOS）
- 至少16GB内存（对于100x100网格的测试案例）
- 50GB以上磁盘空间（用于存储地理数据和输出结果）
- 必备编译器：gcc/gfortran（建议版本9+）
- 并行计算支持：MPI库（如OpenMPI）
- 辅助工具：csh、tar、curl等基础命令
关键软件组件版本建议：
| 组件名称 | 推荐版本 | 备注 |
|---|---|---|
| WPS | 4.1+ | 地理数据处理系统 |
| WRF | 4.1.2+ | 主模型核心 |
| netCDF | 4.7.0+ | 数据格式支持 |
| Jasper | 2.0.0+ | GRIB解码库 |
2.2 静态地理数据下载
地理数据是WRF模拟的基础，需通过WPS的geogrid程序处理：
# 创建数据目录 mkdir -p ~/Build_WRF/WPS_GEOG cd ~/Build_WRF # 下载基础地理数据集（约1.5GB） wget https://www2.mmm.ucar.edu/wrf/src/wps_files/geog_10m.tar.gz tar -zxvf geog_10m.tar.gz mv geog WPS_GEOG常见缺失数据解决方案：
- 运行geogrid.exe时若报错缺少特定数据集（如orogwd_10m）
- 单独下载缺失包并解压到WPS_GEOG目录
- 重新运行geogrid.exe直至所有数据完备
3. WPS配置详解
3.1 namelist.wps关键参数
&share wrf_core = 'ARW', max_dom = 1, start_date = '2000-01-24_12:00:00', end_date = '2000-01-26_00:00:00', interval_seconds = 21600 / &geogrid parent_id = 1, parent_grid_ratio = 1, i_parent_start = 1, j_parent_start = 1, e_we = 100, e_sn = 100, geog_data_res = 'default', dx = 27000, dy = 27000, map_proj = 'mercator', ref_lat = 41.0, ref_lon = 86.0, truelat1 = 30.0, truelat2 = 60.0, stand_lon = 86.0, geog_data_path = '/path/to/WPS_GEOG' /参数解析：
e_we/e_sn：东西/南北向网格点数（决定分辨率）dx/dy：网格间距（米），与地图投影相关map_proj：投影方式（mercator/lambert/polar等）ref_lat/ref_lon：区域中心坐标
3.2 地理数据处理流程
生成网格定义：
./geogrid.exe >& geogrid.log检查生成的geo_em.d01.nc文件：
ncview geo_em.d01.nc气象数据预处理：
- 下载GRIB格式的初始场数据（如GFS/FNL）
- 使用ungrib提取气象要素：
ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable ./link_grib.csh /path/to/GRIB/files ./ungrib.exe
数据插值到网格：
./metgrid.exe生成met_em.d01.*文件供WRF使用
4. WRF主模型运行
4.1 namelist.input配置要点
&time_control run_days = 0, run_hours = 36, start_year = 2000, start_month = 01, start_day = 24, start_hour = 12, end_year = 2000, end_month = 01, end_day = 26, end_hour = 00, interval_seconds = 21600, history_interval = 180, / &domains time_step = 150, time_step_fract_num = 0, time_step_fract_den = 1, max_dom = 1, e_we = 100, e_sn = 100, e_vert = 45, p_top_requested = 5000, num_metgrid_levels = 27, dx = 27000, dy = 27000, grid_id = 1, parent_id = 1, parent_grid_ratio = 1, /物理过程参数选择：
mp_physics：微物理方案（8=Thompson方案）cu_physics：积云参数化（3=Grell-Freitas）bl_pbl_physics：边界层方案（2=YSU）
4.2 模型执行流程
初始化真实数据：
./real.exe生成wrfinput_d01（初始场）和wrfbdy_d01（边界条件）
主模型运行：
mpirun -np 4 ./wrf.exe根据硬件配置调整进程数（-np参数）
运行监控：
- 查看rsl.error.0000日志文件
- 监控输出文件生成进度：
ls -lh wrfout_d01_*
5. 常见错误排查指南
5.1 地理数据相关问题
错误现象：
ERROR: Could not open /path/WPS_GEOG/orogwd_10m/con/index解决方案：
- 确认WPS_GEOG路径正确
- 检查缺失的具体数据集名称
- 从UCAR官网下载补充数据包
- 解压到WPS_GEOG对应目录
5.2 时间设置错误
典型报错：
Timing for input data is not valid排查步骤：
- 检查namelist.wps中：
- start_date/end_date格式
- interval_seconds与输入数据频率一致
- 确认met_em文件时间序列完整
- 验证计算机系统时区设置
5.3 数值不稳定问题
表现特征：
- 运行中出现NaN值
- 模型崩溃报"floating point exception"
调试方法：
- 减小时间步长（time_step）
- 增加扩散系数（diff_opt/km_opt）
- 尝试不同的物理参数化组合
- 检查初始场数据质量
6. 结果分析与可视化
6.1 输出文件结构
WRF生成的主要结果文件：
wrfout_d01_*：主输出文件（netCDF格式）wrfxtrm_d01_*：极值统计文件wrfrst_d01_*：重启文件
关键变量说明：
| 变量名 | 描述 | 单位 |
|---|---|---|
| T2 | 2米气温 | K |
| U10/V10 | 10米风场 | m/s |
| RAINNC | 累积降水 | mm |
| LH | 潜热通量 | W/m² |
6.2 基础可视化方法
使用NCL快速查看结果：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" a = addfile("wrfout_d01_2000-01-24_12:00:00","r") times = wrf_user_getvar(a,"times",-1) do i = 0,dimsizes(times)-1 slp = wrf_user_getvar(a,"slp",i) ; 海平面气压 wrf_smooth_2d(slp, 3) ; 平滑处理 contour = wrf_contour(a,slp,False) plot = wrf_map_overlays(a,contour,True,True) end doPython替代方案（使用xarray+wrf-python）：
import xarray as xr from wrf import getvar, to_np ds = xr.open_dataset("wrfout_d01_2000-01-24_12:00:00") t2 = getvar(ds, "T2") - 273.15 # 转换为摄氏度 lat = to_np(getvar(ds, "XLAT"))[0,:] lon = to_np(getvar(ds, "XLONG"))[:,0]7. 性能优化技巧
7.1 计算资源调配
MPI任务分配建议：
| 网格规模 | 推荐CPU核心数 | 内存需求 |
|---|---|---|
| 100x100 | 4-8 | 16GB |
| 200x200 | 8-16 | 32GB |
| 500x500 | 16-32 | 64GB+ |
优化运行命令示例：
mpirun --map-by core --bind-to core -np 16 ./wrf.exe7.2 I/O性能提升
- 使用并行NetCDF输出：
io_form_history = 2 io_form_restart = 2 nio_tasks_per_group = 4 nio_groups = 2 - 减少输出频率和变量
- 考虑使用时间平均替代瞬时输出
8. 扩展应用方向
掌握基础单域模拟后，可进一步探索：
- 嵌套网格配置（多域模拟）
- 数据同化（WRFDA）应用
- 集合预报系统构建
- 耦合模式开发（如与海洋模型耦合）
- 定制化物理参数化方案
实际项目中遇到的一个典型问题是地理数据路径设置错误，导致geogrid.exe无法找到高程数据。通过添加环境变量GEOG_DATA_PATH可以避免每次修改namelist.wps的麻烦：
export GEOG_DATA_PATH="/path/to/WPS_GEOG"