wrfchem笔记

如果要用上一次化学结果驱动本次模拟,需要将上一次模拟结果命名成 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

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容