Dear WRFhelp,
I am intending to run a 3-domain WRF run (12,4,1.333 km) with the following setup.
Unfortunately, I have some issues when I add the third domain into the simulation.
It exactly stops in correspondence with radt which is set to 12 (when I change it to
a different value for radt, the same happens at the minute specified).
A few other things also I noticed with from WRF-3.9 onwards for the same situation
as mentioned above:
(1) when you set the number of grids nx and ny differently (instead of a square domain,
if I choose a rectangular domain) in the nests.
(2) when you set the grid spacing below 4 km in the nest (4, 1.333 only in the simulation).
I don't understand the reasoning behind this. Could you help me with it. I am attaching
the namelist file below.
Note: If I am running only one domain run, the above issue is not found at all.
Timing for main: time 2014-05-30_12:11:53 on domain 2: 3.77984 elapsed seconds
Timing for main: time 2014-05-30_12:11:54 on domain 3: 0.80333 elapsed seconds
Timing for main: time 2014-05-30_12:11:55 on domain 3: 0.80033 elapsed seconds
Timing for main: time 2014-05-30_12:11:56 on domain 3: 0.79460 elapsed seconds
Timing for main: time 2014-05-30_12:11:56 on domain 2: 3.77581 elapsed seconds
Timing for main: time 2014-05-30_12:11:57 on domain 3: 0.80447 elapsed seconds
Timing for main: time 2014-05-30_12:11:58 on domain 3: 0.79110 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 3: 0.79499 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 2: 3.76707 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 1: 12.68133 elapsed seconds
-----------------------------------
&time_control
run_days = 0,
run_hours = 9,
run_minutes = 0,
run_seconds = 0,
start_year = 2014, 2014, 2014, 2017,
start_month = 05, 05, 05, 08,
start_day = 30, 30, 30, 10,
start_hour = 12, 12, 12, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2014, 2014, 2014, 2017,
end_month = 05, 05, 05, 08,
end_day = 30, 30, 30, 18,
end_hour = 21, 21, 21, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 3600
input_from_file = .true.,.true.,.true.,.true.,
fine_input_stream = 0, 0, 0, 0,
history_interval = 60, 10, 10, 180,
frames_per_outfile = 1, 6, 6, 1,
restart = .false.,
restart_interval = 720,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput2 = 2
io_form_auxinput10 = 2
debug_level = 0
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 60, 60, 60, 60,
io_form_auxinput4 = 2
/
&domains
time_step = 10,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 501, 454, 355, 355,
e_sn = 501, 454, 355, 355,
e_vert = 49, 49, 49, 49,
num_metgrid_levels = 38
num_metgrid_soil_levels = 4
p_top_requested = 10000
lagrange_order = 2
dx = 12000, 4000, 1333.3333, 444.4444,
dy = 12000, 4000, 1333.3333, 444.4444,
grid_id = 1, 2, 3, 4,
parent_id = 0, 1, 2, 3,
i_parent_start = 1, 240, 190, 160,
j_parent_start = 1, 225, 210, 120,
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 0
smooth_option = 0
use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.2, 1.2, 1.2, 1.2,
target_hcfl = .84, .84, .84, .84,
max_step_increase_pct = 5, 31, 31, 31,
starting_time_step = 12, 4, 2, 1,
max_time_step = 36, 12, 1, 1,
min_time_step = 6, 4, 1, 1,
min_time_step_den = 0, 9, 9, 9,
adaptation_domain = 2,
/
&dfi_control
/
&physics
mp_physics = 8, 8, 8, 8,
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 12, 12, 12, 12,
sf_sfclay_physics = 2, 2, 2, 2,
sf_surface_physics = 2, 2, 2, 2,
bl_pbl_physics = 2, 2, 2, 2,
bldt = 0, 0, 0, 0,
cu_physics = 14, 0, 0, 0,
cudt = 0, 0, 0, 0,
kfeta_trigger = 2,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
tmn_update = 0,
lagday = 140,
sst_update = 1,
mp_zero_out = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
slope_rad = 1,
topo_shading = 0,
co2tf = 1
levsiz = 59
paerlev = 29
cam_abs_dim1 = 4
cam_abs_dim2 = 49
cam_abs_freq_s = 3600
prec_acc_dt = 1440, 1440, 1440, 1440,
/
&tc
/
&fdda
grid_fdda = 0, 0, 0, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_end_h = 780, 780, 24, 24,
gfdda_interval_m = 360, 360, 360, 360,
fgdt = 0, 0, 0, 0,
if_no_pbl_nudging_uv = 0, 0, 0, 0,
if_no_pbl_nudging_t = 1, 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1, 1,
if_zfac_uv = 0, 0, 0, 0,
k_zfac_uv = 10, 10, 10, 10,
if_zfac_t = 0, 0, 0, 0,
k_zfac_t = 10, 10, 10, 10,
if_zfac_q = 0, 0, 0, 0,
k_zfac_q = 10, 10, 10, 10,
guv = 0.005, 0.005, 0.0003, 0.0003,
gt = 0.005, 0.005, 0.0003, 0.0003,
gq = 0.005, 0.005, 0.0003, 0.0003,
if_ramping = 1,
dtramp_min = 60.0,
io_form_gfdda = 2,
/
&scm
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
damp_opt = 3,
base_temp = 290.
use_baseparam_fr_nml = .t.
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.025,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
smdiv = 0.1, 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01, 0.01,
time_step_sound = 24, 48, 72, 36,
h_mom_adv_order = 5, 5, 5, 5,
v_mom_adv_order = 3, 3, 3, 3,
h_sca_adv_order = 5, 5, 5, 5,
v_sca_adv_order = 3, 3, 3, 3,
non_hydrostatic = .true.,.true.,.true.,.true.,
moist_adv_opt = 2, 2, 2, 2,
scalar_adv_opt = 2, 2, 2, 2,
tke_adv_opt = 2, 2, 2, 2,
/
&fire
/
&bdy_control
spec_bdy_width = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false.,.false.,.false.,
nested = .false.,.true.,.true.,.true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am intending to run a 3-domain WRF run (12,4,1.333 km) with the following setup.
Unfortunately, I have some issues when I add the third domain into the simulation.
It exactly stops in correspondence with radt which is set to 12 (when I change it to
a different value for radt, the same happens at the minute specified).
A few other things also I noticed with from WRF-3.9 onwards for the same situation
as mentioned above:
(1) when you set the number of grids nx and ny differently (instead of a square domain,
if I choose a rectangular domain) in the nests.
(2) when you set the grid spacing below 4 km in the nest (4, 1.333 only in the simulation).
I don't understand the reasoning behind this. Could you help me with it. I am attaching
the namelist file below.
Note: If I am running only one domain run, the above issue is not found at all.
Timing for main: time 2014-05-30_12:11:53 on domain 2: 3.77984 elapsed seconds
Timing for main: time 2014-05-30_12:11:54 on domain 3: 0.80333 elapsed seconds
Timing for main: time 2014-05-30_12:11:55 on domain 3: 0.80033 elapsed seconds
Timing for main: time 2014-05-30_12:11:56 on domain 3: 0.79460 elapsed seconds
Timing for main: time 2014-05-30_12:11:56 on domain 2: 3.77581 elapsed seconds
Timing for main: time 2014-05-30_12:11:57 on domain 3: 0.80447 elapsed seconds
Timing for main: time 2014-05-30_12:11:58 on domain 3: 0.79110 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 3: 0.79499 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 2: 3.76707 elapsed seconds
Timing for main: time 2014-05-30_12:12:00 on domain 1: 12.68133 elapsed seconds
-----------------------------------
&time_control
run_days = 0,
run_hours = 9,
run_minutes = 0,
run_seconds = 0,
start_year = 2014, 2014, 2014, 2017,
start_month = 05, 05, 05, 08,
start_day = 30, 30, 30, 10,
start_hour = 12, 12, 12, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2014, 2014, 2014, 2017,
end_month = 05, 05, 05, 08,
end_day = 30, 30, 30, 18,
end_hour = 21, 21, 21, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 3600
input_from_file = .true.,.true.,.true.,.true.,
fine_input_stream = 0, 0, 0, 0,
history_interval = 60, 10, 10, 180,
frames_per_outfile = 1, 6, 6, 1,
restart = .false.,
restart_interval = 720,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput2 = 2
io_form_auxinput10 = 2
debug_level = 0
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 60, 60, 60, 60,
io_form_auxinput4 = 2
/
&domains
time_step = 10,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 501, 454, 355, 355,
e_sn = 501, 454, 355, 355,
e_vert = 49, 49, 49, 49,
num_metgrid_levels = 38
num_metgrid_soil_levels = 4
p_top_requested = 10000
lagrange_order = 2
dx = 12000, 4000, 1333.3333, 444.4444,
dy = 12000, 4000, 1333.3333, 444.4444,
grid_id = 1, 2, 3, 4,
parent_id = 0, 1, 2, 3,
i_parent_start = 1, 240, 190, 160,
j_parent_start = 1, 225, 210, 120,
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 0
smooth_option = 0
use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.2, 1.2, 1.2, 1.2,
target_hcfl = .84, .84, .84, .84,
max_step_increase_pct = 5, 31, 31, 31,
starting_time_step = 12, 4, 2, 1,
max_time_step = 36, 12, 1, 1,
min_time_step = 6, 4, 1, 1,
min_time_step_den = 0, 9, 9, 9,
adaptation_domain = 2,
/
&dfi_control
/
&physics
mp_physics = 8, 8, 8, 8,
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 12, 12, 12, 12,
sf_sfclay_physics = 2, 2, 2, 2,
sf_surface_physics = 2, 2, 2, 2,
bl_pbl_physics = 2, 2, 2, 2,
bldt = 0, 0, 0, 0,
cu_physics = 14, 0, 0, 0,
cudt = 0, 0, 0, 0,
kfeta_trigger = 2,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
tmn_update = 0,
lagday = 140,
sst_update = 1,
mp_zero_out = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
slope_rad = 1,
topo_shading = 0,
co2tf = 1
levsiz = 59
paerlev = 29
cam_abs_dim1 = 4
cam_abs_dim2 = 49
cam_abs_freq_s = 3600
prec_acc_dt = 1440, 1440, 1440, 1440,
/
&tc
/
&fdda
grid_fdda = 0, 0, 0, 0,
gfdda_inname = "wrffdda_d<domain>",
gfdda_end_h = 780, 780, 24, 24,
gfdda_interval_m = 360, 360, 360, 360,
fgdt = 0, 0, 0, 0,
if_no_pbl_nudging_uv = 0, 0, 0, 0,
if_no_pbl_nudging_t = 1, 1, 1, 1,
if_no_pbl_nudging_q = 1, 1, 1, 1,
if_zfac_uv = 0, 0, 0, 0,
k_zfac_uv = 10, 10, 10, 10,
if_zfac_t = 0, 0, 0, 0,
k_zfac_t = 10, 10, 10, 10,
if_zfac_q = 0, 0, 0, 0,
k_zfac_q = 10, 10, 10, 10,
guv = 0.005, 0.005, 0.0003, 0.0003,
gt = 0.005, 0.005, 0.0003, 0.0003,
gq = 0.005, 0.005, 0.0003, 0.0003,
if_ramping = 1,
dtramp_min = 60.0,
io_form_gfdda = 2,
/
&scm
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
damp_opt = 3,
base_temp = 290.
use_baseparam_fr_nml = .t.
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.025,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
smdiv = 0.1, 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01, 0.01,
time_step_sound = 24, 48, 72, 36,
h_mom_adv_order = 5, 5, 5, 5,
v_mom_adv_order = 3, 3, 3, 3,
h_sca_adv_order = 5, 5, 5, 5,
v_sca_adv_order = 3, 3, 3, 3,
non_hydrostatic = .true.,.true.,.true.,.true.,
moist_adv_opt = 2, 2, 2, 2,
scalar_adv_opt = 2, 2, 2, 2,
tke_adv_opt = 2, 2, 2, 2,
/
&fire
/
&bdy_control
spec_bdy_width = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false.,.false.,.false.,
nested = .false.,.true.,.true.,.true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/