Xochitl H.
New member
Hello,
I am attempting to run a 3 domain simulation for a two week period over the Houston area. My simulation produces just over 14 hours of output before a segmentation fault occurs. I have checked my rsl files for cfl errors and none are present. From looking at my output, there appears to be a single grid point with unrealistic values (i.e. TSK and T2 at this point are significantly lower than all other points and HFX is negative). My thought is that this point is causing the simulation to crash.
Do you know why this point would have such values and, if so, how to fix it?
See below for my namelist.input. I will try to attach pictures of the fields in question as soon as possible (the upload function is currently not working for me). I am using ERA5 forcing data and the default MODIS for my landuse.
Thank you for any and all help. It is greatly appreciated!
UPDATE: While I have still not discovered why my simulation stops working, I have found that my IVGTYP at the point in question is 15 (snow and ice). This is the only point in my domain that is classified this way. Could this be producing the error I am having?
&time_control
run_days = 0,
run_hours = 360,
run_minutes = 0,
run_seconds = 0,
start_year = 2017, 2017, 2017,
start_month = 08, 08, 08,
start_day = 18, 18, 18,
start_hour = 00, 00, 00,
end_year = 2017, 2017, 2017,
end_month = 09, 09, 09,
end_day = 02, 02, 02,
end_hour = 00, 00, 00,
interval_seconds = 21600,
auxinput4_interval = 360, 360, 360,
auxinput4_inname = "wrflowinp_d<domain>",
input_from_file = .true., .true., .true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1, 1, 1,
restart = .false.,
restart_interval = 2880,
history_outname = "./wrf_modis/wrfout_d<domain>_<date>"
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_auxinput4 = 2,
io_form_boundary = 2,
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 265, 626, 298,
e_sn = 220, 511, 298,
e_vert = 48, 48, 48,
sfcp_to_sfcp = .false.
p_top_requested = 5000,
eta_levels = 1.000,0.994,0.988,0.981,0.973,0.964,0.954,0.943,0.931,
0.918,0.904,0.889,0.873,0.856,0.838,0.819,0.799,0.778,
0.756,0.733,0.709,0.684,0.658,0.631,0.603,0.574,0.544,
0.513,0.481,0.448,0.414,0.379,0.343,0.307,0.271,0.236,
0.202,0.170,0.140,0.112,0.086,0.062,0.040,0.025,0.015,
0.008,0.003,0.000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
interp_type = 2,
lagrange_order = 2,
zap_close_levels = 500,
lowest_lev_from_sfc = .false.,
force_sfc_in_vinterp = 1,
dx = 15000, 3000, 1000,
dy = 15000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 68, 273,
j_parent_start = 1, 62, 192,
parent_grid_ratio = 1, 5, 3,
parent_time_step_ratio = 1, 5, 3,
feedback = 0,
smooth_option = 0,
/
&physics
physics_suite = 'CONUS'
mp_physics = 8, 8, 8,
cu_physics = 3, 0, 0,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
bl_pbl_physics = 5, 5, 5,
bl_mynn_edmf = 1, 1, 1,
ishallow = 0,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
radt = 15, 15, 15,
bldt = 0, 0, 0,
cudt = 0, 0, 0,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 0, 0, 1,
sst_update = 1,
do_radar_ref = 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
prec_acc_dt = 60, 60, 60,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 2, 2, 2,
km_opt = 4, 4, 4,
diff_6th_opt = 0,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
smdiv = 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01,
epssm = 0.1, 0.5, 0.5,
non_hydrostatic = .true., .true.,.true.,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am attempting to run a 3 domain simulation for a two week period over the Houston area. My simulation produces just over 14 hours of output before a segmentation fault occurs. I have checked my rsl files for cfl errors and none are present. From looking at my output, there appears to be a single grid point with unrealistic values (i.e. TSK and T2 at this point are significantly lower than all other points and HFX is negative). My thought is that this point is causing the simulation to crash.
Do you know why this point would have such values and, if so, how to fix it?
See below for my namelist.input. I will try to attach pictures of the fields in question as soon as possible (the upload function is currently not working for me). I am using ERA5 forcing data and the default MODIS for my landuse.
Thank you for any and all help. It is greatly appreciated!
UPDATE: While I have still not discovered why my simulation stops working, I have found that my IVGTYP at the point in question is 15 (snow and ice). This is the only point in my domain that is classified this way. Could this be producing the error I am having?
&time_control
run_days = 0,
run_hours = 360,
run_minutes = 0,
run_seconds = 0,
start_year = 2017, 2017, 2017,
start_month = 08, 08, 08,
start_day = 18, 18, 18,
start_hour = 00, 00, 00,
end_year = 2017, 2017, 2017,
end_month = 09, 09, 09,
end_day = 02, 02, 02,
end_hour = 00, 00, 00,
interval_seconds = 21600,
auxinput4_interval = 360, 360, 360,
auxinput4_inname = "wrflowinp_d<domain>",
input_from_file = .true., .true., .true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1, 1, 1,
restart = .false.,
restart_interval = 2880,
history_outname = "./wrf_modis/wrfout_d<domain>_<date>"
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_auxinput4 = 2,
io_form_boundary = 2,
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 265, 626, 298,
e_sn = 220, 511, 298,
e_vert = 48, 48, 48,
sfcp_to_sfcp = .false.
p_top_requested = 5000,
eta_levels = 1.000,0.994,0.988,0.981,0.973,0.964,0.954,0.943,0.931,
0.918,0.904,0.889,0.873,0.856,0.838,0.819,0.799,0.778,
0.756,0.733,0.709,0.684,0.658,0.631,0.603,0.574,0.544,
0.513,0.481,0.448,0.414,0.379,0.343,0.307,0.271,0.236,
0.202,0.170,0.140,0.112,0.086,0.062,0.040,0.025,0.015,
0.008,0.003,0.000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
interp_type = 2,
lagrange_order = 2,
zap_close_levels = 500,
lowest_lev_from_sfc = .false.,
force_sfc_in_vinterp = 1,
dx = 15000, 3000, 1000,
dy = 15000, 3000, 1000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 68, 273,
j_parent_start = 1, 62, 192,
parent_grid_ratio = 1, 5, 3,
parent_time_step_ratio = 1, 5, 3,
feedback = 0,
smooth_option = 0,
/
&physics
physics_suite = 'CONUS'
mp_physics = 8, 8, 8,
cu_physics = 3, 0, 0,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
bl_pbl_physics = 5, 5, 5,
bl_mynn_edmf = 1, 1, 1,
ishallow = 0,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
radt = 15, 15, 15,
bldt = 0, 0, 0,
cudt = 0, 0, 0,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 0, 0, 1,
sst_update = 1,
do_radar_ref = 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
prec_acc_dt = 60, 60, 60,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 2, 2, 2,
km_opt = 4, 4, 4,
diff_6th_opt = 0,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
smdiv = 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01,
epssm = 0.1, 0.5, 0.5,
non_hydrostatic = .true., .true.,.true.,
time_step_sound = 4, 4, 4,
h_mom_adv_order = 5, 5, 5,
v_mom_adv_order = 3, 3, 3,
h_sca_adv_order = 5, 5, 5,
v_sca_adv_order = 3, 3, 3,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
Last edited: