[Wrf-users] WRF/chem error with MEGAN bio-emissions and RACM_VBS mechanism

Chris Tessum tess0050 at umn.edu
Fri Feb 8 09:02:58 MST 2013


Hello,

I am trying to run WRF/chem v3.4 with the RACM_VBS chemical mechanism and
MEGAN biogenic emissions. Everything works fine in the beginning, but after
one hour of simulation time the program crashes when it tries to write out
the results, giving this error when I set the debug level to 1000:

 d03 2004-12-27_01:00:00 module_io.F (md_calls.m4) : in
wrf_put_dom_ti_integer_s
 ca
 output_wrf.b writing 2d real sebio_oli
  date 2004-12-27_01:00:00
  ds    903596857   903539941   904037630
  de    903966924   903372531   903691536
  ps    904183828   904434164   903202198
  pe    903706240   903644900   903560598
  ms    903548162          -4   903853002
  me    903802871   903809105   903731388
 d03 2004-12-27_01:00:00 module_io.F: in wrf_write_field
wrf.exe: Rank 0:1: MPI_Gatherv: Internal MPI error

If I set bio_emiss_opt=1 (switching from the MEGAN bio-emissions to the
simple bio-emissions), the problem goes away. If I change chem_opt to 2
(the default RADM2 mechanism), that also causes the problem to go away.
Does anyone have any ideas about what could be causing this problem or how
to fix it? I have included the namelist below:

&time_control
run_days                      = 1,
run_hours                     = 0,
run_minutes                   = 0,
run_seconds                   = 0,
start_year                    = 2004, 2004, 2004,
start_month                   = 12, 12, 12,
start_day                     = 27, 27, 27,
start_hour                    = 00, 00, 00,
start_minute                  = 00, 00, 00,
start_second                  = 00, 00, 00,
end_year                      = 2005, 2005, 2005,
end_month                     = 02, 02, 02,
end_day                       = 16, 16, 16,
end_hour                      = 00, 00, 00,
end_minute                    = 00, 00, 00,
end_second                    = 00, 00, 00,
interval_seconds              = 21600,
input_from_file               = .true., .true., .true.,
history_interval_h            = 24, 24, 24,
auxhist1_interval             = 60, 60, 60,
frames_per_outfile            = 1, 1, 1,
frames_per_auxhist1           = 24, 24, 24,
restart                       = .false.,
restart_interval_h            = 24,
io_form_history               = 2,
io_form_auxhist1              = 2,
io_form_restart               = 2,
io_form_input                 = 2,
io_form_boundary              = 11,
io_form_auxinput5             = 0,
io_form_auxinput6             = 0,
io_form_auxinput7             = 11,
io_form_auxinput12            = 0,
frames_per_auxinput5          = 24, 24, 24,
auxinput5_interval_m          = 60, 60, 60,
frames_per_auxinput7          = 1, 1, 1,
auxinput7_interval_m          = 60, 60, 60,
auxinput7_inname              = wrffirechemi_d<domain>_<date>,
auxinput12_inname             = wrf_chem_input_d<domain>,
nocolons                      = .true.,
debug_level                   = 10001,
/

&domains
time_step                     = 225,
use_adaptive_time_step        = .false.,
time_step_fract_num           = 0,
time_step_fract_den           = 1,
max_dom                       = 3,
e_we                          = 149, 55, 73,
e_sn                          = 113, 55, 91,
e_vert                        = 28, 28, 28,
p_top_requested               = 5000,
num_metgrid_levels            = 27,
num_metgrid_soil_levels       = 4,
dx                            = 36000, 12000, 2000,
dy                            = 36000, 12000, 2000,
grid_id                       = 1, 2, 3,
parent_id                     = 1, 1, 2,
i_parent_start                = 1, 77, 21,
j_parent_start                = 1, 65, 19,
parent_grid_ratio             = 1, 3, 6,
parent_time_step_ratio        = 1, 3, 6,
feedback                      = 1,
smooth_option                 = 0,
/

&physics
mp_physics                    = 3, 3, 3,
ra_lw_physics                 = 3, 3, 3,
ra_sw_physics                 = 3, 3, 3,
radt                          = 36, 12, 2,
sf_sfclay_physics             = 1, 1, 1,
sf_surface_physics            = 2, 2, 2,
bl_pbl_physics                = 1, 1, 1,
bldt                          = 0, 0, 0,
cu_physics                    = 5, 5, 5,
cu_rad_feedback               = .true., .true., .true.,
cudt                          = 0, 0, 0,
cu_diag                       = 1, 1, 1,
isfflx                        = 1,
ifsnow                        = 1,
icloud                        = 1,
surface_input_source          = 1,
num_soil_layers               = 4,
sf_urban_physics              = 0, 0, 0,
/

&fdda
grid_fdda                     = 1, 1, 1,
gfdda_inname                  = wrffdda_d<domain>,
gfdda_interval_m              = 360, 360, 360,
gfdda_end_h                   = 8760, 8760, 8760,
io_form_gfdda                 = 11,
/

&dynamics
w_damping                     = 1,
diff_opt                      = 1,
km_opt                        = 4,
diff_6th_opt                  = 0, 0, 0,
diff_6th_factor               = 0.12, 0.12, 0.12,
base_temp                     = 290.0,
damp_opt                      = 0,
zdamp                         = 5000.0, 5000.0, 5000.0,
dampcoef                      = 0.2, 0.2, 0.2,
khdif                         = 0, 0, 0,
kvdif                         = 0, 0, 0,
non_hydrostatic               = .true., .true., .true.,
chem_adv_opt                  = 2, 2, 2,
moist_adv_opt                 = 2, 2, 2,
scalar_adv_opt                = 2, 2, 2,
tke_adv_opt                   = 2, 2, 2,
/

&bdy_control
spec_bdy_width                = 5,
spec_zone                     = 1,
relax_zone                    = 4,
specified                     = .true., .false., .false.,
nested                        = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group           = 0,
nio_groups                    = 1,
/

&chem
kemit                         = 23,
chem_opt                      = 108, 108, 108,
bioemdt                       = 30, 30, 30,
photdt                        = 30, 30, 30,
chemdt                        = 7.5, 2.5, 0.83,
emiss_inpt_opt                = 1, 1, 1,
emiss_opt                     = 3, 3, 3,
io_style_emissions            = 2,
chem_in_opt                   = 0, 0, 0,
phot_opt                      = 2, 2, 2,
gas_drydep_opt                = 1, 1, 1,
aer_drydep_opt                = 1, 1, 1,
bio_emiss_opt                 = 3, 3, 3,
dust_opt                      = 0,
dmsemis_opt                   = 0,
seas_opt                      = 0,
gas_bc_opt                    = 1, 1, 1,
gas_ic_opt                    = 1, 1, 1,
aer_bc_opt                    = 1, 1, 1,
aer_ic_opt                    = 1, 1, 1,
gaschem_onoff                 = 1, 1, 1,
aerchem_onoff                 = 1, 1, 1,
wetscav_onoff                 = 0, 0, 0,
cldchem_onoff                 = 0, 0, 0,
vertmix_onoff                 = 1, 1, 1,
chem_conv_tr                  = 1, 1, 1,
biomass_burn_opt              = 1, 1, 1,
plumerisefire_frq             = 30, 30, 30,
have_bcs_chem                 = .true., .true., .true.,
aer_ra_feedback               = 0, 0, 0,
/

Thank you for your time,

-- 
Chris Tessum

Doctoral Candidate
Environmental Engineering
University of Minnesota
+1 (612) 568 2238
tess0050 at umn.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20130208/0f94dd45/attachment-0001.html 


More information about the Wrf-users mailing list