如果要用上一次化学结果驱动本次模拟,需要将上一次模拟结果命名成 wrf_chem_input_01,
并放在相同文件夹下,namelist设置如下
! io_form_auxinput12 = 2,
! auxinput12_inname = 'wrf_chem_input',
./anthro_emis <input_sap_agr
anthro_emis直接生成wrfchemi,无需convert_emiss这一步了
anthro_emis中需要调整data_yrs_offset = 0, 例如模拟年是2021, 清单年是2019, 这里需要设成2
wrfchem里需要real两次:
第一次,纯气象real, 需要关闭chemopt(=0),
第二次,在生成wrfchemi之后,把wrfchemi链接到run的路径下,打开chemopt,然后real, 这里得到气象化学耦合场
马上要毕业了,毕业前将以前的一些简要经验发出来。由于水平不高,基本都是直接翻译的用户手册,也算是抛砖引玉,欢迎论坛各位批评指正。
**[资料](http://www.06climate.com/)下载**
> SST ftp://polar.ncep.noaa.gov/pub/history/sst
> FNL
**WRF 前处理**
* **运行 geogrid**
> 与一般情况相同
* **运行 ungrib**
> 分别对**FNL** 资料与**SST** 资料进行解压。
> * **首先 **使用 **link_grib.csh** 链接 **fnl **资料
> ./link_grib.csh fnl*
> * 链接合适的Vtable
> ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
> * 在**namelist.wps**中 **&share** 部分作如下改动
> &share
>
> start_date = '2006-09-23_00:00:00',
>
> end_date = '2006-09-24_00:00:00',
>
> interval_seconds = 21600,
> * 执行 **ungrid.exe**
> ./ungrib.exe
>
> 产生
>
> FILE:2006-09-23_00
>
> FILE:2006-09-23_06
>
> FILE:2006-09-23_12
>
> FILE:2006-09-23_18
>
> FILE:2006-09-24_00
> * **然后**,使用 **link_grib.csh** 链接**SST **数据
> > ./link_grib.csh rtg_sst*
>
> * 链接合适的 Vtable(选择 Vtable.SST)
> ln -sf ungrib/Variable_Tables/Vtable.SST Vtable
> * 编辑 **namelist.wps**中的 **&share** 和 **&ungrib **部分
> 注意:SST数据的间隔时间以天为计。所以需要将SST插值至6小时间隔,以匹配FNL资料。通过将interval_seconds=21600,可以将SST数据插值至6小时。通过要注意改变prefix,这样才不会覆盖掉之前解压的fnl。
>
> &share
>
> start_date = '2003-07-02_00:00:00',
>
> end_date = '2003-07-03_00:00:00',
>
> interval_seconds = 21600,
>
> &ungrib
>
> prefix = 'SST',
> * 执行 **ungrib.exe**
> ./ungrib.exe
>
> 产生
>
> SST:2006-09-23_00
>
> SST:2006-09-23_06
>
> SST:2006-09-23_12
>
> SST:2006-09-23_18
>
> SST:2006-09-24_00
* **运行 metgrid**
> * 更改**namelist.wps **中的 **&metgrid** 部分
> &metgrid
>
> fg_name = 'FILE', 'SST'
> * 执行 **metgrid**
> ./metgrid.exe
>
> 产生
>
> met_em.d01.2006-09-23_00:00:00.nc
>
> met_em.d01.2006-09-23_06:00:00.nc
>
> met_em.d01.2006-09-23_12:00:00.nc
>
> met_em.d01.2006-09-23_18:00:00.nc
>
> met_em.d01.2006-09-24_00:00:00.nc
**运行 WRF 主模块**
> 编辑 **namelist.input**,在**&time_control** 和**&physics **部分添加
>
> &time_control
>
> auxinput4_inname = "wrflowinp_d<domain>"
>
> auxinput4_interval = 360,
>
> io_form_auxinput4 = 2,
>
> &physics
>
> sst_update = 1,
>
> 这些参数将在模式模拟中每6小时更新一次sst。
>
> 不要更改auxinput4_inname。
* **运行 real**
> * 链接 **metgrid**输出文件
> > > ln -sf ../../../WPS/met_em* .
> * 执行 **real.exe**
> 运行成功,产生下列文件,(wrflowinp_d01包含SST信息)
>
> wrfbdy_d01
>
> wrfinput_d01
>
> **wrflowinp_d01**
* **运行 WRF**
> * 与一般情况一样,执行 wrf.exe