一.对研究区域选定以及数据预处理
本步操作主要在/home/zll/WRF/WPS目录下完成.
(1) 编辑namelist.wps文件(红字是需要修改的部分)
&share
wrf_core = 'ARW',
max_dom = 2, % 代表嵌套层数
start_date = '2014-01-01_00:00:00','2014-01-01_00:00:00', % 模拟开始时间
end_date = '2014-03-31_18:00:00','2014-03-31_18:00:00',%模拟结束时间
interval_seconds = 21600% 驱动数据时间间隔(s)
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3, % 这里是3,代表内层分辨率是外层的3倍
i_parent_start = 1, 31,
j_parent_start = 1, 17,
e_we = 74, 112, % 外层和内层东西方向的格点数
e_sn = 61, 97, % 外层和内层南北方向的格点数
geog_data_res = '10m','2m',
dx = 30000, % 外层x方向的一个个点代表的距离(m)
dy = 30000, % 外层x方向的一个个点代表的距离(m)
map_proj = 'lambert',
ref_lat = 47.5, % 中心点纬度
ref_lon = 90, % 中心点经度
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = 90, % 中心点经度
geog_data_path ='/home/zll/WRF/geog' % geog文件所在的位置
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
% 设置完namelist.wps文件之后保存一下,然后运行nclplotgrids.ncl,可以生成研究区域的范围。如果生成的研究区偏大,可以通过设置e_we和e_sn来调整。
(2)运行命令:./geogrid.exe。运行成功后将会生成几个nc文件。
(3)链接Vtable.GFS,并连接下载的驱动文件
./link_grib.csh/home/zll/WRF/Data/fnl* ./
ln -sf /home/zll/WRF/WPS/ungrib/Variable_Tables/Vtable.GFS Vtable
(4)解压气象数据
./ungrib.exe
(5)把气象数据水平插入到模式领域内
./metgrid.exe
二、WRF主程序的运行
这部分工作主要在/home/zll/WRF/WRFV3/run文件夹下完成
(1)修改namelist.input文件,红色是需要修改的部分,参数化方案在这个文件里修改
&time_control
run_days = 89, % 模型运行天数,1~3月共90天
run_hours = 18, % 每天0:00~18:00
run_minutes = 0, % 分钟
run_seconds = 0, % 秒
start_year = 2014, 2014, 2000,
start_month = 01, 01, 01,
start_day = 01, 01, 24,
start_hour = 00, 00, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2014, 2014, 2000,
end_month = 03, 03, 01,
end_day = 31, 31, 25,
end_hour = 18, 18, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60, % 模式输出时间间隔(分钟)
frames_per_outfile = 1000,1000, 1000, %运行1000小时生成一个文件,防止服务器中断。
restart = .false.,
restart_interval = 1000,%运行1000小时生成一个文件,防止服务器中断。
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/
&domains % 这部分参数参照namelist.wps文件设置
time_step = 180, % 这个设置为空间分辨率(30km)的6倍
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 = 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, 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
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 1, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
/
&fdda
/
&dynamics
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 = 0,
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,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,.false.,.false.,
nested = .false.,.true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
(2)链接 WPS 过程中产生的met_em.d01.2009-12-25_00:00:00.nc 等 nc
文件
ln -sf/home/zll/WRF/WPS/met_em.d* ./
(3)运行 real.exe
./real.exe
(4)运行 wrf.exe,用于生成最终数据。
./wrf.exe
备注:如果在超算上运行,需提交作业:
bsub -n 16 -qNormal ./wrf.exe
其中16代表使用的节点数。
查询目前运行的作业:bjobs -u all
转载本文请联系原作者获取授权,同时请注明本文来自张乐乐科学网博客。
收藏
分享
分享到: