Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

WRF stops in conjunction with "radt" in nests

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

rameshv

New member
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,
/
 
Top