小白跑WRF第七天,运行WRF案例
官方提供了一个运行案例及数据,但是作者并没有使用,因为跑官方的案例失败了,但是自己去下了一些fnl的数据,跑通了。
这个帖子不介绍参数设置,目的是怕初学者搞晕。
所设计的数据大家可以留下邮箱,作者提供。或者自行准备相同格式的数据。
准备一、创建DATA文件夹
在Build_WRF下创建一个DATA文件夹用于存放要用的数据。当然,你也可以自己搞个名字,自己找个路径,只要后面合并的时候能搞清楚Linux的路径命令即可。
这里就不演试了。
准备二、替换namelist.input内容
该文件路径:Build_WRF\WRF\test\em_real
方法一样,复制下面的,然后全部替换。
&time_control
run_days = 0,
run_hours = 30,
run_minutes = 0,
run_seconds = 0,
start_year = 2021, 2000, 2000,
start_month = 10, 01, 01,
start_day = 15, 24, 24,
start_hour = 00, 12, 12,
end_year = 2021, 2000, 2000,
end_month = 10, 01, 01,
end_day = 17, 25, 25,
end_hour = 18, 12, 12,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 7200,
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, 112, 94,
e_sn = 61, 97, 91,
e_vert = 33, 33, 33,
p_top_requested = 5000,
num_metgrid_levels = 34,
num_metgrid_soil_levels = 4,
dx = 30000, 10000, 3333.33,
dy = 30000, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 31, 30,
j_parent_start = 1, 17, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
physics_suite = 'CONUS'
mp_physics = -1, -1, -1,
cu_physics = -1, -1, 0,
ra_lw_physics = -1, -1, -1,
ra_sw_physics = -1, -1, -1,
bl_pbl_physics = -1, -1, -1,
sf_sfclay_physics = -1, -1, -1,
sf_surface_physics = -1, -1, -1,
radt = 30, 30, 30,
bldt = 0, 0, 0,
cudt = 5, 5, 5,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
/
&fdda
/
&dynamics
hybrid_opt = 2,
w_damping = 0,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
gwd_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
specified = .true.
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
接下来开始,才是日常运行WRF该有的步骤
第一步、运行静态地理数据
一、进入WPS文件夹
cd ~/Build_WRF/WPS
二、修改namelist.wps
这里建议直接通过记事本去修改,具体的内容后面会出一期专门的帖子来介绍里面各参数的含义,这里先带大家跑一遍。
该文件位置:\Build_WRF\WPS\namelist.wps
注意:改动的前提是大家按照之前帖子中的namelist.wps已做了覆盖
改动如下图
简单介绍一下:
其中max_dom = 1,代表嵌套层数为1。
后面的start_date和end_date代表了起始和结束时间,代表了2层,如果3层,后面继续加。
照理说现在设置的是1,只需要改第一列的时间,但是本人是本着改了2个更保险。
那为什么
start_date = ‘2020-07-27_00:00:00’
end_date = ‘2020-07-27_18:00:00’
因为一会儿我要跑的数据是从 '2020-07-27_00:00:00’到’2020-07-27_18:00:00’的。
修改完成后保存退出。
三、生成地理数据
这些操作都在WPS下进行。
输入命令:
./geogrid.exe
当出现下面文字代表成功
如果出现下面的内容,请去下载对应的地理信息包,方法在之前的帖子以及交代。
成功后,输入命令,检查地理信息文件
ls -lah geo_em.d01.nc
第二步、处理气象数据
一、将数据放入DATA文件夹中
这一步,可以在DATA文件夹下再创建子文件夹(像作者一样),也可以不创建,直接DATA。但是路径要自己搞清楚。
作者放的是 ‘2020-07-27_00:00:00’到’2020-07-27_18:00:00’,逐6小时的FNL资料,grib2格式。
提示:大家务必要知道自己数据的格式,时间等信息,这对设置各类namelist文件有用,具体后面帖子再介绍
二、链接数据
这些操作都在WPS下进行。
输入命令:
./link_grib.csh ~/Build_WRF/DATA/20200727/fnl
ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable
二、启动ungrib.exe
这些操作都在WPS下进行。
输入命令:
./ungrib.exe
成功则出现:
不成功就根据报错找问题。
检查文件信息 输入:
ls -lah FILE*
得到处理好的气象数据
这些文件再WPS根目录下。
第三步、融合气象和地理数据
这些操作都在WPS下进行。
这个就是启动metgrid.exe程序
./metgrid.exe
到这里,WPS部分,WRF的预处理部分完成了。
第四步、修改WRF中的namelist.input
跳转至em_real文件夹(后面的操作都是WRF\test\em_real下进行的)
cd ~/Build_WRF/WRF/test/em_real/
下面修改namelist.input,文件位置:\Build_WRF\WRF\test\em_real\namelist.input
修改的是起止时间以及需要运行模拟几个小时的数据。一层所以只需要改第一列。
第五步、链接WPS到WRF
后面的操作都是WRF\test\em_real下进行的
ln -sf ~/Build_WRF/WPS/met_em* .
mpirun -np 1 ./real.exe
查看文件
ls -alh wrfbdy_d01 wrfinput_d01
第六步、运行WRF(时间较久,耐心等待)
后面的操作都是WRF\test\em_real下进行的
mpirun -np 2 ./wrf.exe
命令中的2,代表了2核并行运行,看自己电脑的能力,也可以8核并行。
成功后会出现
这样,WRF算是跑完了,结果也成功了,用python可以直接读取。
第七步、ARWpost处理数据文件
进入Build_WRF/ARWpost/
cd ~/Build_WRF/ARWpost/
修改namelist.arwpost的起止时间
运行ARWpost.exe
./ARWpost.exe
最后得到
这样,最后在设定好的文件夹里找到对应的ctl和dat文件就可以用grads处理啦
大家如果需要数据可以邮件我,也可以自己找数据尝试。