Hi:
My data for initial and boundary conditions are FNL 0.25deg x 0.25deg. WRF 4.2 is aborting shortly after startup (please see below). I am also including my namelist.
I had a look at the variable GHT in my met_em.* files. Over complex terrain areas there is an inconsistency in the lowest two levels: the second level stays in several areas below the first level. I was told that with ERA5 forcing data I should only use 4-point interpolators for GHT in the metgrid table and it worked fine. However, this procedure failed with NCEP FNL analysis. A priori real.exe seems to fix this problem as in wrfinput the pressure/geopotential height against level number show a decreasing/increasing behaviour for increasing lower levels. However, wrf.exe makes a few initial steps and then crashes. I wonder if I should verify for integrity some additional variables in wrfinput or wrfbdy to see if the problem in the lowest levels of met_em.* becomes somehow transferred to the initial and boundary conditions. It is remarkable that SFCLAY is the last subroutine before the fault, as it deals with the lowest levels. Any suggestion about how to solve this metgrid issue ?
My data for initial and boundary conditions are FNL 0.25deg x 0.25deg. WRF 4.2 is aborting shortly after startup (please see below). I am also including my namelist.
Code:
d02 2019-09-11_00:00:40 DEBUG wrf_timetoa(): returning with str = [2019-09-11_00:00:40]
d02 2019-09-11_00:00:40 call radiation_driver
d02 2019-09-11_00:00:40 Top of Radiation Driver
d02 2019-09-11_00:00:40 calling inc/HALO_PWP_inline.inc
d02 2019-09-11_00:00:40 call surface_driver
d02 2019-09-11_00:00:40 in SFCLAY
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error:
#0 0x2aaaabc5333f in ???
#1 0x271f2dd in ???
#2 0x2722a46 in ???
#3 0x2727c52 in ???
#4 0x1f0c1ba in ???
#5 0x17bab87 in ???
#6 0x12f9901 in ???
#7 0x11b7ae4 in ???
#8 0x4735fa in ???
#9 0x473bda in ???
#10 0x406213 in ???
#11 0x405bcc in ???
#12 0x2aaaabc3f494 in ???
#13 0x405c03 in ???
#14 0xffffffffffffffff in ???
Code:
&time_control
run_days = 0,
run_hours = 36,
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 = 12, 12,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 21600,
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 = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 364, 697,
e_sn = 382, 574,
e_vert = 100, 100,
p_top_requested = 100,
num_metgrid_levels = 34,
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,
max_dz = 600.,
auto_levels_opt = 2,
zap_close_levels = 0.1,
lagrange_order = 1,
smooth_cg_topo = .true.,
/
&physics
mp_physics = 3, 3,
ra_lw_physics = 1, 1,
ra_sw_physics = 1, 1,
radt = 30, 30,
sf_sfclay_physics = 1, 1,
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.,
iso_temp = 200.,
base_pres_strat = 1300,
base_lapse_strat = -40.,
damp_opt = 3,
zdamp = 10000., 10000.,
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,
/