## WRF实例运行1
****单层嵌套模型****
step1:定义模型模拟域
- 下载静态地理数据并解压(以本人计算机实践为例)
下载地址WRFV4:https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.html#specifi
下载地址WRFV3:
https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog_V3.html
解压后存储地址:$HOME/Public/Build_WRF/WPS_GEOG
可以不用全部下载,运行geogrid时,程序告诉你缺什么就可以下载什么。
- 修改namelist.wps内容
cd ~/Public/Build_WRF/WRFV4/WPS
vi namelist.wps
修改namelist.wps中的&geogrid和&share中与时间无关的字段,均可按照自己想设置的域范围进行修改。
&share
wrf_core = 'ARW',
max_dom = 1, #单层粗域处理&geogrid #定义模拟域投影及所有模型网格尺寸和位置。具体每项的含义可以参考前文“the WRF preprocessing system(WPS)介绍”
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 31,
j_parent_start = 1, 17,
e_we = 74, 112,
e_sn = 61, 97,geog_data_res = 'default','default',
dx = 30000,
dy = 30000,
map_proj = 'lambert',
ref_lat = 34.83,
ref_lon = -81.03,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = -98.0,
geog_data_path = '/home/hsy/Public/Build_WRF/WPS_GEOG/'
最主要的是一定要输入WPS_GEOG的绝对路径,注意每行后面都有','作为结尾。
- 运行geogrid
./geogrid.exe
运行成功后。
nc文件可以使用ncview查看。
ncview geogrid_output/20000124/geo_em.d01.nc
step2:用ungrib从网格数据中提取气象字段
- 连接气象数据,使用./link_grib.csh脚本连接到自己解压数据后的存储地址。(本文使用的是20000124-20000125的gfs数据,可自行下载,也可以去博主资源主页https://download.csdn.net/download/A18040554844/18937434下载。)
./link_grib.csh ../../../grid_data/JAN00_GFS_FNL/JAN00/
- 连接Vtable。(本文使用的是Vtable.GFS,可根据下载的数据选择不同的Vtable,更可以根据自己下载的数据源创建新的Vtable。)
每个Vtable包含7个或11个字段,具体取决于Vtable是用于GRIB Edition 1数据源还是用于GRIB Edition 2数据源。Vtable的字段可以分为三类:描述如何在GRIB文件中标识数据的字段,描述ungrib和metgrid程序如何标识数据的字段,以及特定于GRIB Edition 2的字段。具体如何创建Vtable可以自行根据官网教程设置。
ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
- 修改namelist.wps中与时间相关的share字段以及ungrib字段。
&share#修改与模拟时间相关的部分
start_date = '2000-01-24_12:00:00',
end_date = '2000-01-25_12:00:00',
interval_seconds = 21600,
&ungrib
out_format = 'WPS',
prefix = 'ungrib_output/20000124/FILE',#输出的临时文件位置为ungrib_output/20000124/,FILE为文件前缀。
- 运行ungrib。
./ungrib.exe
运行成功后。
step3:用metgrid水平插值气象数据
- 修改namelist.wps中&metgrid部分,这里可以再次检查一下&share中时间设置是否恰当。
&metgrid
fg_name = 'ungrib_output/20000124/FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = 'metgrid_output/'
- 运行metgrid。
./metgrid.exe
运行成功后。
附单层粗域namelist.wps设置以供参考(WRF4.0)。
step4:运行real.exe。
- 连接metgrid产生的nc数据。
cd ../WRF/test/em_real/
ls -l ../../../WPS/metgrid_output/met_em.d01.2000-01-2*
ln -sf ../../../WPS/metgrid_output/met_em.d01.2000-01-2* .
连接成功后。
- 修改namelist.input内容
主要修改&time_control和&domains部分。
&time_control
run_days = 0,
run_hours = 24,
run_minutes = 0,
run_seconds = 0,
start_year = 2000,
start_month = 01,
start_day = 24,
start_hour = 12,
end_year = 2000,
end_month = 01,
end_day = 25,
end_hour = 12,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 180, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_auxinput2 = 2,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
/&domains
time_step = 180,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 74, , ,
e_sn = 61, , ,
e_vert = 30, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 30000, 10000, 3333.33,
dy = 30000, 10000, 3333.33,
grid_id = 1,
parent_id = 0,
i_parent_start = 1, , ,
j_parent_start = 1, , ,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
- 运行real
./real.exe
#运行完以后,立马检查real输出。
tail rsl.error.0000
#or
tail rsl.out.0000
#or
vi rsl.error.0000
real运行报错,证明namelist.input设置与气象数据输入不同,请耐心根据报错内容修改namelist.input中&domain字段,直至出现SUCCESS COMPLETE REAL_EM INIT以及
wrfbdy_d01和 wrfinput_d01两个文件。
step5:运行wrf模型。
- 运行wrf
./wrf.exe
运行成功后,会出现SUCCESS COMPLETE WRF提示,并且出现wrfout_d0*文件。