WRF实例运行（2）

## WRF实例运行2

参考自https://www2.mmm.ucar.edu/wrf/OnLineTutorial/CASES/SingleDomain/index.php
****单层嵌套模型****

本案例研究适用于 2016-10-06_00 至 2016-10-08_00 期间。

我们要设置的域如下所示。

step1：解压输入的 GRIB 数据(ungrib.exe)

我们计划将此案例研究用于许多不同的模型运行，包括将 SST 数据添加到模型运行的示例。所以在我们开始之前，让我们先解压数据。

1. 下载静态地理数据并储存（以作者计算机实践为例）

下载地址：https://www2.mmm.ucar.edu/wrf/TUTORIAL_DATA/matthew_1deg.tar.gz

解压后的储存地址：~/Public/data_test/single/matthew

GFS（Global Forecast System）是NCEP的示范产品，数据类型为grib2数据。

      2. 连接GFS Vtable以及grib数据

分辨率：1°的全球数据；输出频率为6小时；27层大气压（1000-10hpa；不包括地表）

ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
./link_grib.csh ../../data_test/single/matthew/fnl

连接建立后输出内容：

GRIBFILE.AAA -> ../DATA/matthew/fnl_20161006_00_00.grib2
GRIBFILE.AAB -> ../DATA/matthew/fnl_20161006_06_00.grib2
GRIBFILE.AAC -> ../DATA/matthew/fnl_20161006_12_00.grib2
GRIBFILE.AAD -> ../DATA/matthew/fnl_20161006_18_00.grib2
GRIBFILE.AAE -> ../DATA/matthew/fnl_20161007_00_00.grib2
GRIBFILE.AAF -> ../DATA/matthew/fnl_20161007_06_00.grib2
GRIBFILE.AAG -> ../DATA/matthew/fnl_20161007_12_00.grib2
GRIBFILE.AAH -> ../DATA/matthew/fnl_20161007_18_00.grib2
GRIBFILE.AAI -> ../DATA/matthew/fnl_20161008_00_00.grib2

        3. 修改namelist.wps

max_dom = 1
start_date = '2016-10-06_00:00:00',
end_date = '2016-10-08_00:00:00',
interval_seconds = 21600,
prefix = 'FILE',

4.运行ungrib建立中间文件

字段作为二维水平数据块写入中间文件。每块水平板只包含一个单一变量(即500mb 相对湿度、表面 t 等)。这些文件被写入为未格式化的 FORTRAN 记录。对于每个水平数据块，写入了许多记录。记录数取决于中间文件版本。

./ungrib.exe

FILE:2016-10-06_00
FILE:2016-10-06_06
FILE:2016-10-06_12
FILE:2016-10-06_18
FILE:2016-10-07_00
FILE:2016-10-07_06
FILE:2016-10-07_12
FILE:2016-10-07_18
FILE:2016-10-08_00

熟悉中间文件。

./util/rd_intermediate.exe FILE:2016-10-06_00

 step2：设置模型域(geogrid.exe)

1、确保在WPS文件目录下，修改namelist.wps

max_dom = 1 
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 91,
e_sn = 100,
geog_data_res = 'default',
dx = 27000,
dy = 27000,
map_proj = 'mercator',
ref_lat = 28.00,
ref_lon = -75.00,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -75.0,
geog_data_path = 'Your WPS_GEOG data location'最好是绝对路径

 2、查看域的位置。

ncl util/plotgrids_new.ncl

 

3、运行geogrid建立该域的静态数据

注意：现在运行 geogrid.exe 为这个域创建静态数据： 注意：这将覆盖您以前运行时可能拥有的任何 geo_em.d01.nc 文件， 因此，如果您想保留旧文件，请将其保存在其他地方。

./geogrid.exe

4、检查静态数据geo_em.d01.nc是否创建成功，并检查其文件内容。

ncdump -h geo_em.d01.nc

 

step3：将输入数据插入我们的模型域 (metgrid.exe)

1、修改namelist.wps内容

无需对namelist.wps文件进行任何其他更改。只需确保为 Matthew 案例正确设置开始和结束日期，并且将 fg_name 设置为 FILE。

2、运行metgrid.exe 用于在我们的模型域中插入输入数据：

./metgrid.exe

 met_em.d01.2016-10-06_00:00:00.nc
met_em.d01.2016-10-06_06:00:00.nc
met_em.d01.2016-10-06_12:00:00.nc
met_em.d01.2016-10-06_18:00:00.nc
met_em.d01.2016-10-07_00:00:00.nc
met_em.d01.2016-10-07_06:00:00.nc
met_em.d01.2016-10-07_12:00:00.nc
met_em.d01.2016-10-07_18:00:00.nc
met_em.d01.2016-10-08_00:00:00.nc

3、检查该文件内容。

ncdump -h met_em.d01.2016-10-06_00:00:00.nc

step4：运行模型 (real.exe & wrf.exe)

1、在em_real文件夹内连接metgrid.exe创建的met_em文件。

cd ../WRF/test/em_real/
ln -sf ../../../WPS/met_em.d01.2016-10* .

2、编辑namelist.wps内容。

 

run_days = 0,
run_hours = 48,
run_minutes = 0,
run_seconds = 0,
start_year = 2016,
start_month = 10,
start_day = 06,
start_hour = 00,
end_year = 2016,
end_month = 10,
end_day = 08,
end_hour = 00,
interval_seconds = 21600
input_from_file = .true.,
history_interval = 180,
frames_per_outfile = 1000,
restart = .false.,
restart_interval = 1440,
time_step = 150,
max_dom = 1,
e_we = 91,
e_sn = 100,
s_vert = 1,
e_vert = 35,
num_metgrid_levels = 32
dx = 27000,
dy = 27000,

注意：上面的 restart_interval 设置 - 我们计划稍后重新启动，因此我们需要设置此运行以预期重新启动。这将创建额外的重新启动输出文件，稍后我们可以从中重新启动模型。

3、运行real.exe。

wrfinput_d01

wrfbdy_d01

 

4、运行wrf.exe。

输出：

wrfout_d01_2016-10-06_00:00:00
wrfrst_d01_2016-10-07_00:00:00
wrfrst_d01_2016-10-08_00:00:00

检查文件内容：

 ncdump -h wrfout_d01_2016-10-06_00:00:00
 ncdump -v Times wrfout_d01_2016-10-06_00:00:00