wrf运行参数设置

小怪兽...

各位前辈好，我刚接触wrf，一个人摸索，求帮忙。Ubuntu14.04+WRF3.6.1，1个物理cpu、16核，安装了mpich库及WRFV3(dmpar版)，可并行运算。
实验目的：获得长时间序列的高空间分辨率（3-1km?）气象数据，时间分辨率只要到日就可以。
namelist.wps和namelist.input设置如下，请各位帮忙看下，如何设置比较合适、怎么修改能运行的快点。

namelist.wps:
&share
wrf_core = 'ARW',
max_dom =3,
start_date = '2017-09-01_00:00:00','2017-09-01_00:00:00','2017-09-01_00:00:00',
end_date   = '2019-01-01_00:00:00','2019-01-01_00:00:00','2019-01-01_00:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/

&geogrid
parent_id         =   1,   1,   2,
parent_grid_ratio =   1,   5,   5,
i_parent_start    =   1,   13,  38,   
j_parent_start    =   1,   13,  40,   
e_we              =   51,131,276,
e_sn              =   46,126,226,
geog_data_res     =   '10m','2m',
dx = 30000,
dy = 30000,
map_proj = 'lambert',
ref_lat   =  36.8,
ref_lon   =  109.5,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon =  109.5,
geog_data_path = '/home/linux/data/geog'
/

&ungrib
out_format = 'WPS',
prefix = '/media/linux/18252035589/WPS/2018/FILE'
/

&metgrid
fg_name = '/media/linux/18252035589/WPS/2018/FILE'
io_form_metgrid = 2,
opt_output_from_metgrid_path='/media/linux/18252035589/WPS/2018'
/



namelist.input
&time_control
run_days                            = 487,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017,2017,2017,
start_month                         = 09,   09,   09,
start_day                           = 01,   01,   01,
start_hour                          = 00,   00,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2019, 2019, 2019,
end_month                           = 01,   01,   01,
end_day                             = 01,   01,   01,
end_hour                            = 00,   00,   00,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 1440,1440,1440,
frames_per_outfile                  = 1,1,1,
restart                             = .false.,
restart_interval                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 51,   131,   276,
e_sn                                = 46,   126,   226,
e_vert                              = 30,   30,    30,
p_top_requested                     = 10000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 30000, 6000,  1200,
dy                                  = 30000, 6000,  1200,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     13,    38,
j_parent_start                      = 1,     13,    40,
parent_grid_ratio                   = 1,     5,     5,
parent_time_step_ratio              = 1,     5,     5,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 3,     3,     3,
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 1,     1,     1,
radt                                = 27,    27,    27,
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,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 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,
/

小其其格

按照这个模拟时长和空间分辨率，考虑到你的计算机硬件条件，我估计全部模拟结束需要……嗯，10年？

小怪兽...

小其其格 发表于 2019-9-9 19:42
按照这个模拟时长和空间分辨率，考虑到你的计算机硬件条件，我估计全部模拟结束需要……嗯，10年？

谢谢您的回复，对模拟时长没什么经验。您看如果我改成双层嵌套，dx=27000, 比例1:9，降到3公里，情况怎么样

crazytiy

这个时长有点夸张

小其其格

小怪兽... 发表于 2019-9-10 09:12
谢谢您的回复，对模拟时长没什么经验。您看如果我改成双层嵌套，dx=27000, 比例1:9，降到3公里，情况怎么 ...

建议借个超算或者大型机来跑吧

小怪兽...

小其其格 发表于 2019-9-10 22:02
建议借个超算或者大型机来跑吧

好的，谢谢您

管麦斯

时间太长了，越到后面积分误差会越大。做什么要模拟这么久啊

一大碗年糕

长时间模拟可以用区域气候模式RegCM之类的，WRF的话需要调试好稳定
e. Regional climate case at 10 – 30 km grid sizes (e.g. used in NCAR’s regional climate runs):
mp_physics                          = 6,
ra_lw_physics                       = 3,
ra_sw_physics                       = 3,
radt                                = 30,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 2,
bl_pbl_physics                      = 1,
bldt                                = 0,
cu_physics                          = 1,
cudt                                = 5,
sst_update                          = 1,
tmn_update                          = 1,
sst_skin                            = 1,
bucket_mm                           = 100.0,
bucket_J                            = 1.e9,
ptop_requested                      = 1000,
e_vert                              = 51,
spec_bdy_width                      = 10,
spec_zone                           = 1,
relax_zone                          = 9,
spec_exp                            = 0.33,