I am using WRF 4.2 to perform a high resolution simulation with adaptive stepsize. My run always stops around 30 min after the initial time in coincidence with the first radiation calculations by WRF (radt=30), when it is calling cldfra1. Boundary and initial conditions have been prepared from ERA5 model levels up to 0.01 hPa. I checked a large variety of combinations of short and long wave radiation and microphysics options in the namelist and the simulation always fails at the same point. I also tested radt=6 following the suggestion of selecting a value close to dx in km and as expected the program crashed at minute 6 instead of 30. I have seen a similar problematic case with cldfra1 with no solution (wrfforum.com/viewtopic.php?f=21&t=11692&sid=b67aa07976ae9bf0c452061ad852ec78). I was not able to find any clear description on how cldfra1 works. Below I am including my namelist and the last lines of one of the rsl.error files (debug_level=1000) in case that someone can give me a suggestion. Thank you in advance
Code:
&time_control
run_days = 0,
run_hours = 31,
run_minutes = 0,
run_seconds = 0,
start_year = 2019, 2019,
start_month = 9, 9,
start_day = 11, 11,
start_hour = 0, 0,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2019, 2019,
end_month = 9, 9,
end_day = 12, 12,
end_hour = 7, 7,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 3600,
input_from_file = .true., .true.,
history_interval = 180, 6,
history_outname = "/scratch/peter/wrfout_d<domain>_<date>"
frames_per_outfile = 1000, 110,
restart = .false.,
restart_interval = 5000,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 1000,
/
&domains
time_step = 9,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 364, 697,
e_sn = 382, 574,
e_vert = 150, 150,
p_top_requested = 1.0002,
num_metgrid_levels = 138,
num_metgrid_soil_levels = 4,
dx = 9000, 3000,
dy = 9000, 3000,
grid_id = 1, 2,
parent_id = 1, 1,
i_parent_start = 1, 46,
j_parent_start = 1, 82,
parent_grid_ratio = 1, 3,
parent_time_step_ratio = 1, 3,
feedback = 1,
smooth_option = 0,
use_adaptive_time_step = .true.,
step_to_output_time = .true.,
target_cfl = 1.2, 1.2,
max_step_increase_pct = 25, 50,
starting_time_step = 3, 1,
starting_time_step_den = 10, 10,
max_time_step = 45, 15,
min_time_step = 3, 1,
min_time_step_den = 100, 100,
adaptation_domain = 1,
/
&physics
mp_physics = 5, 5,
ra_lw_physics = 99, 99,
ra_sw_physics = 99, 99,
radt = 30, 30,
co2tf = 1,
sf_sfclay_physics = 91, 91,
sf_surface_physics = 2, 2,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 0, 0,
cudt = 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.,
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.2, 0.2,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 1, 1,
scalar_adv_opt = 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Timing for wrf_patch_to_global_generic: 0.00000 elapsed seconds
<stdin> writing 2d real sst_input Status = 0
Timing for wrf_ext_write_field: 0.24428 elapsed seconds
d02 2019-09-11_00:30:10+07/** output_wrf: calling wrf_iosync
d02 2019-09-11_00:30:10+07/** module_io.F: in wrf_iosync
d02 2019-09-11_00:30:10+07/** output_wrf: back from wrf_iosync
d02 2019-09-11_00:30:10+07/** output_wrf: end, fid = 3
d02 2019-09-11_00:30:10+07/** med_hist_out: opened /scratch/peter/wrfout_d02_2019-09-11_00:30:10 as DATASET=HISTORY
d02 2019-09-11_00:30:10+07/** in med_latbound_in
d02 2019-09-11_00:30:10+07/** module_integrate: calling solve interface
d02 2019-09-11_00:30:10+07/** grid spacing, dt, time_step_sound= 3000.00000 15.0000000 4
d02 2019-09-11_00:30:10+07/** calling inc/HALO_EM_MOIST_OLD_E_7_inline.inc
d02 2019-09-11_00:30:10+07/** calling inc/PERIOD_BDY_EM_MOIST_OLD_inline.inc
d02 2019-09-11_00:30:10+07/** call rk_step_prep
d02 2019-09-11_00:30:10+07/** calling inc/HALO_EM_A_inline.inc
d02 2019-09-11_00:30:10+07/** calling inc/PERIOD_BDY_EM_A_inline.inc
d02 2019-09-11_00:30:10+07/** call rk_phys_bc_dry_1
d02 2019-09-11_00:30:10+07/** call init_zero_tendency
d02 2019-09-11_00:30:10+07/** calling inc/HALO_EM_PHYS_A_inline.inc
d02 2019-09-11_00:30:10+07/** call phy_prep
d02 2019-09-11_00:30:10+07/** DEBUG wrf_timetoa(): returning with str = [2019-09-11_00:30:10]
d02 2019-09-11_00:30:10+07/** call radiation_driver
d02 2019-09-11_00:30:10+07/** Top of Radiation Driver
d02 2019-09-11_00:30:10+07/** CALL cldfra1
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
#0 0x2af4596bd33f in ???
#1 0x16d9d23 in ???
#2 0x17b279f in ???
#3 0x12f9901 in ???
#4 0x11b7ae4 in ???
#5 0x4735fa in ???
#6 0x473bda in ???
#7 0x406213 in ???
#8 0x405bcc in ???
#9 0x2af4596a9494 in ???
#10 0x405c03 in ???
#11 0xffffffffffffffff in ???