一.对研究区域选定以及数据预处理

本步操作主要在/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

转载本文请联系原作者获取授权，同时请注明本文来自张乐乐科学网博客。

收藏

分享

分享到: