WRF-Chem基本流程（1）

一、确定研究区域、时间范围和嵌套层数

先去WRF Domain Wizard（WRF Domain Wizard (jiririchter.github.io)）绘制研究区大致范围，左侧会自动生成网格信息，之后会在namelist.wps中填入这些信息。

二、修改namelist.wps

&share
 wrf_core = 'ARW',
 max_dom = 2,
 start_date = '2021-06-20_00:00:00','2021-06-20_00:00:00',
 end_date   = '2021-07-10_00:00:00','2021-07-10_00:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

wrf_core是选择核心，WRF有两种核，分别是：

1. ARW（Advanced Research WRF）核心：

   • 由NCAR（美国国家大气研究中心）主导开发。
   • 适用于研究和高级应用。
   • 提供更多的物理过程和数值算法选项，允许用户进行更复杂和定制化的模拟。
   • 基于有限差分方法，适合高分辨率模拟。

2. NMM（Nonhydrostatic Mesoscale Model）核心：

   • 由NOAA（美国国家海洋和大气管理局）的NCEP（国家环境预报中心）开发。
   • 主要用于操作性天气预报。
   • 采用有限体积方法，具有更强的数值稳定性，适合大规模并行计算。
   • 设计上更为简化，适合实际预报需求。

max_dom指的是嵌套层数，这里设置2就是两层嵌套。

start_date和end_date分别是指模拟的起始和终止时间，因为嵌套层数是2，所以要写两列。

interval_seconds指的是输入数据（气象数据，如全球分析数据或再分析数据）的时间间隔，单位是秒。这里使用的是FNL再分析数据，时间间隔为6个小时，因此写21600秒。FNL数据可以在这里下载（WRF - Free Data (ucar.edu)）

io_form_geogrid指的是地理数据文件的输出格式，就是在运行完geogrid.exe后的文件格式。1表示二进制格式（bin），2表示网格格式（nc）。

&geogrid
 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  53,
 j_parent_start    =   1,  25,
 e_we              =  60, 220,
 e_sn              =  44, 214,
 geog_data_res = 'default','default',
 dx = 3000,
 dy = 3000,
 map_proj = 'lambert',
 ref_lat   =  34.242,
 ref_lon   = 117.186,
 truelat1  =  34.242,
 truelat2  =  34.242,
 stand_lon = 117.186,
 geog_data_path = '/home/user/wrfchem4/geog/WPS_GEOG'
/

parent_grid_ratio表示子网格相对父网格分辨率的比率，第一列设置1，第二列设置3表示子网格的分辨率要比父网格的分辨率大3倍，父网格设置为3000，则子网格为1000。

geog_data_res用于指定地理数据的分辨率，详细信息见（http://t.csdnimg.cn/KcTAl）

geog_data_path表示静态地理数据的存放位置，静态地理数据可以在这里下载：Static Data Downloads (ucar.edu)

其他的参数根据WRF Domain Wizard设置即可。

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/
&metgrid
 fg_name = 'FILE'
/

这里是一些输出文件和中间文件格式，无需修改。

至此namelist修改完成，可以使用ncview查看设置的范围，这里不做演示。