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等基础命令关键软件组件版本建议组件名称推荐版本备注WPS4.1地理数据处理系统WRF4.1.2主模型核心netCDF4.7.0数据格式支持Jasper2.0.0GRIB解码库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微物理方案8Thompson方案cu_physics积云参数化3Grell-Freitasbl_pbl_physics边界层方案2YSU4.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_*重启文件关键变量说明变量名描述单位T22米气温KU10/V1010米风场m/sRAINNC累积降水mmLH潜热通量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替代方案使用xarraywrf-pythonimport 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核心数内存需求100x1004-816GB200x2008-1632GB500x50016-3264GB优化运行命令示例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