Segmentation Fault when enabling both WRF-Chem (chem_opt=10) and Urban Canopy Model (sf_urban_physics=1) in WRF v4.5

Maxrr

New member
Hi everyone,

I am encountering a persistent segmentation fault when attempting to run a coupled simulation with both the chemistry module and the Urban Canopy Model (UCM) enabled in WRF-Chem v4.5.

Here are the details of my configuration and the testing I've done:

  1. The Issue: The model crashes with a "segmentation fault" after approximately 6 minutes of simulation time.
  2. Case A (Success): Chemistry OFF (chem_opt = 0), Urban Physics ON (sf_urban_physics = 1). The model runs successfully.
  3. Case B (Success): Chemistry ON (chem_opt = 10), Urban Physics OFF (sf_urban_physics = 0). The model runs successfully.
  4. Case C (Failure): Chemistry ON (chem_opt = 10) AND Urban Physics ON (sf_urban_physics = 1). This leads to the crash.

Debugging Steps Taken:
  1. I have checked the rsl.error.* files, but they only report a generic "Segmentation Fault" without pointing to a specific module.
  2. I have increased the debug_level but haven't identified a clear trigger.
  3. My domain settings and static data (WPS) seem fine since both modules work independently.
Has anyone encountered compatibility issues between chem_opt=10 and sf_urban_physics=1 in version 4.5? Could this be related to memory limits (stack size), or is there a known conflict between the urban tile processing and the chemistry species allocation in the boundary layer?

Any suggestions on how to further isolate the error or specific namelist.input adjustments would be greatly appreciated.

Thank you in advance for your help!
 
namelist.input:
&physics
physics_suite = 'CONUS'
mp_physics = 2, 2,
cu_physics = 3, 3,
ra_lw_physics = 4, 4,
ra_sw_physics = 4, 4,
bl_pbl_physics = 1, 1,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
radt = 27, 27,
bldt = 0, 0,
cudt = 5, 5,
cu_diag = 1, 1,
progn = 1, 1,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 1, 1,
!fractional_seaice = 1,
/

&fdda
/

&dynamics
hybrid_opt = 2,
w_damping = 0,
diff_opt = 2, 2,
km_opt = 4, 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
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,
gwd_opt = 1, 0,
/

&bdy_control
spec_bdy_width = 5,
specified = .true.
/

&grib2
/

&chem
kemit = 11,
chem_opt = 10, 10,
chem_in_opt = 0, 0,
bioemdt = 30, 30,
photdt = 27, 27,
chemdt = 5., 5.,
io_style_emissions = 1,
emiss_opt = 4, 4,
emiss_inpt_opt = 101, 101,
bio_emiss_opt = 3, 3,
ne_area = 300,
phot_opt = 2, 2,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 1,
!biomass_burn_opt = 3, 3,
scale_fire_emiss = .false., .false.,
plumerisefire_frq = 60, 60,
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 2,
gas_bc_opt = 1, 1,
gas_ic_opt = 1, 1,
aer_bc_opt = 1, 1,
aer_ic_opt = 1, 1,
gaschem_onoff = 1, 1,
aerchem_onoff = 1, 1,
wetscav_onoff = 1, 1,
cldchem_onoff = 1, 1,
vertmix_onoff = 1, 1,
chem_conv_tr = 1, 1,
aer_ra_feedback = 1, 1,
have_bcs_chem = .true., .true.,
/
 
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2023-07-02_01:01:12
d01 2023-07-02_01:01:12 Error trying to read metadata
d01 2023-07-02_01:01:12 Input data is acceptable to use:
**WARNING** Time in input file not being checked **WARNING**
[a3103n07:6772 :0:6772] Caught signal 11 (Segmentation fault: address not mapped to object at address 0xfffffffe11623040)
==== backtrace (tid: 6772) ====
0 0x00000000042e7b11 module_sf_sfclayrev_mp_psim_stable_() ???:0
1 0x00000000042e2f14 module_sf_sfclayrev_mp_sfclayrev1d_() ???:0
2 0x00000000042e0ae6 module_sf_sfclayrev_mp_sfclayrev_() ???:0
3 0x000000000393cfb6 module_surface_driver_mp_surface_driver_() ???:0
4 0x000000000312c992 module_first_rk_step_part1_mp_first_rk_step_part1_() ???:0
5 0x00000000021be0af solve_em_() ???:0
6 0x0000000001f3b890 solve_interface_() ???:0
7 0x000000000061ba9f module_integrate_mp_integrate_() ???:0
8 0x0000000000416e11 module_wrf_top_mp_wrf_run_() ???:0
9 0x0000000000416dcf MAIN__() ???:0
10 0x0000000000416d62 main() ???:0
11 0x00000000000223d5 __libc_start_main() ???:0
12 0x0000000000416c69 _start() ???:0
=================================
 
Back
Top