<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Hello,</span><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">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:</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><div><font face="courier new, monospace"> d03 2004-12-27_01:00:00 module_io.F (md_calls.m4) : in wrf_put_dom_ti_integer_s</font></div>
<div><font face="courier new, monospace"> ca</font></div><div><font face="courier new, monospace"> output_wrf.b writing 2d real sebio_oli</font></div><div><font face="courier new, monospace"> date 2004-12-27_01:00:00</font></div>
<div><font face="courier new, monospace"> ds 903596857 903539941 904037630</font></div><div><font face="courier new, monospace"> de 903966924 903372531 903691536</font></div><div><font face="courier new, monospace"> ps 904183828 904434164 903202198</font></div>
<div><font face="courier new, monospace"> pe 903706240 903644900 903560598</font></div><div><font face="courier new, monospace"> ms 903548162 -4 903853002</font></div><div><font face="courier new, monospace"> me 903802871 903809105 903731388</font></div>
<div><font face="courier new, monospace"> d03 2004-12-27_01:00:00 module_io.F: in wrf_write_field</font></div><div><font face="courier new, monospace">wrf.exe: Rank 0:1: MPI_Gatherv: Internal MPI error</font></div><div><br>
</div><div>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:</div>
<div><br></div><div><div><font face="courier new, monospace">&time_control</font></div><div><font face="courier new, monospace">run_days = 1,</font></div><div><font face="courier new, monospace">run_hours = 0,</font></div>
<div><font face="courier new, monospace">run_minutes = 0,</font></div><div><font face="courier new, monospace">run_seconds = 0,</font></div><div><font face="courier new, monospace">start_year = 2004, 2004, 2004, </font></div>
<div><font face="courier new, monospace">start_month = 12, 12, 12, </font></div><div><font face="courier new, monospace">start_day = 27, 27, 27, </font></div><div><font face="courier new, monospace">start_hour = 00, 00, 00, </font></div>
<div><font face="courier new, monospace">start_minute = 00, 00, 00, </font></div><div><font face="courier new, monospace">start_second = 00, 00, 00, </font></div><div><font face="courier new, monospace">end_year = 2005, 2005, 2005, </font></div>
<div><font face="courier new, monospace">end_month = 02, 02, 02, </font></div><div><font face="courier new, monospace">end_day = 16, 16, 16, </font></div><div><font face="courier new, monospace">end_hour = 00, 00, 00, </font></div>
<div><font face="courier new, monospace">end_minute = 00, 00, 00, </font></div><div><font face="courier new, monospace">end_second = 00, 00, 00, </font></div><div><font face="courier new, monospace">interval_seconds = 21600,</font></div>
<div><font face="courier new, monospace">input_from_file = .true., .true., .true., </font></div><div><font face="courier new, monospace">history_interval_h = 24, 24, 24, </font></div><div><font face="courier new, monospace">auxhist1_interval = 60, 60, 60, </font></div>
<div><font face="courier new, monospace">frames_per_outfile = 1, 1, 1, </font></div><div><font face="courier new, monospace">frames_per_auxhist1 = 24, 24, 24, </font></div><div><font face="courier new, monospace">restart = .false.,</font></div>
<div><font face="courier new, monospace">restart_interval_h = 24,</font></div><div><font face="courier new, monospace">io_form_history = 2,</font></div><div><font face="courier new, monospace">io_form_auxhist1 = 2,</font></div>
<div><font face="courier new, monospace">io_form_restart = 2,</font></div><div><font face="courier new, monospace">io_form_input = 2,</font></div><div><font face="courier new, monospace">io_form_boundary = 11,</font></div>
<div><font face="courier new, monospace">io_form_auxinput5 = 0,</font></div><div><font face="courier new, monospace">io_form_auxinput6 = 0,</font></div><div><font face="courier new, monospace">io_form_auxinput7 = 11,</font></div>
<div><font face="courier new, monospace">io_form_auxinput12 = 0,</font></div><div><font face="courier new, monospace">frames_per_auxinput5 = 24, 24, 24, </font></div><div><font face="courier new, monospace">auxinput5_interval_m = 60, 60, 60, </font></div>
<div><font face="courier new, monospace">frames_per_auxinput7 = 1, 1, 1, </font></div><div><font face="courier new, monospace">auxinput7_interval_m = 60, 60, 60, </font></div><div><font face="courier new, monospace">auxinput7_inname = wrffirechemi_d<domain>_<date>,</font></div>
<div><font face="courier new, monospace">auxinput12_inname = wrf_chem_input_d<domain>,</font></div><div><font face="courier new, monospace">nocolons = .true.,</font></div><div><font face="courier new, monospace">debug_level = 10001,</font></div>
<div><font face="courier new, monospace">/</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&domains</font></div><div><font face="courier new, monospace">time_step = 225,</font></div>
<div><font face="courier new, monospace">use_adaptive_time_step = .false.,</font></div><div><font face="courier new, monospace">time_step_fract_num = 0,</font></div><div><font face="courier new, monospace">time_step_fract_den = 1,</font></div>
<div><font face="courier new, monospace">max_dom = 3,</font></div><div><font face="courier new, monospace">e_we = 149, 55, 73, </font></div><div><font face="courier new, monospace">e_sn = 113, 55, 91, </font></div>
<div><font face="courier new, monospace">e_vert = 28, 28, 28, </font></div><div><font face="courier new, monospace">p_top_requested = 5000,</font></div><div><font face="courier new, monospace">num_metgrid_levels = 27,</font></div>
<div><font face="courier new, monospace">num_metgrid_soil_levels = 4,</font></div><div><font face="courier new, monospace">dx = 36000, 12000, 2000, </font></div><div><font face="courier new, monospace">dy = 36000, 12000, 2000, </font></div>
<div><font face="courier new, monospace">grid_id = 1, 2, 3, </font></div><div><font face="courier new, monospace">parent_id = 1, 1, 2, </font></div><div><font face="courier new, monospace">i_parent_start = 1, 77, 21, </font></div>
<div><font face="courier new, monospace">j_parent_start = 1, 65, 19, </font></div><div><font face="courier new, monospace">parent_grid_ratio = 1, 3, 6, </font></div><div><font face="courier new, monospace">parent_time_step_ratio = 1, 3, 6, </font></div>
<div><font face="courier new, monospace">feedback = 1,</font></div><div><font face="courier new, monospace">smooth_option = 0,</font></div><div><font face="courier new, monospace">/</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&physics</font></div><div><font face="courier new, monospace">mp_physics = 3, 3, 3, </font></div><div>
<font face="courier new, monospace">ra_lw_physics = 3, 3, 3, </font></div><div><font face="courier new, monospace">ra_sw_physics = 3, 3, 3, </font></div><div><font face="courier new, monospace">radt = 36, 12, 2, </font></div>
<div><font face="courier new, monospace">sf_sfclay_physics = 1, 1, 1, </font></div><div><font face="courier new, monospace">sf_surface_physics = 2, 2, 2, </font></div><div><font face="courier new, monospace">bl_pbl_physics = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">bldt = 0, 0, 0, </font></div><div><font face="courier new, monospace">cu_physics = 5, 5, 5, </font></div><div><font face="courier new, monospace">cu_rad_feedback = .true., .true., .true., </font></div>
<div><font face="courier new, monospace">cudt = 0, 0, 0, </font></div><div><font face="courier new, monospace">cu_diag = 1, 1, 1, </font></div><div><font face="courier new, monospace">isfflx = 1,</font></div>
<div><font face="courier new, monospace">ifsnow = 1,</font></div><div><font face="courier new, monospace">icloud = 1,</font></div><div><font face="courier new, monospace">surface_input_source = 1,</font></div>
<div><font face="courier new, monospace">num_soil_layers = 4,</font></div><div><font face="courier new, monospace">sf_urban_physics = 0, 0, 0, </font></div><div><font face="courier new, monospace">/</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&fdda</font></div><div><font face="courier new, monospace">grid_fdda = 1, 1, 1, </font></div><div>
<font face="courier new, monospace">gfdda_inname = wrffdda_d<domain>,</font></div><div><font face="courier new, monospace">gfdda_interval_m = 360, 360, 360, </font></div><div><font face="courier new, monospace">gfdda_end_h = 8760, 8760, 8760, </font></div>
<div><font face="courier new, monospace">io_form_gfdda = 11,</font></div><div><font face="courier new, monospace">/</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&dynamics</font></div>
<div><font face="courier new, monospace">w_damping = 1,</font></div><div><font face="courier new, monospace">diff_opt = 1,</font></div><div><font face="courier new, monospace">km_opt = 4,</font></div>
<div><font face="courier new, monospace">diff_6th_opt = 0, 0, 0, </font></div><div><font face="courier new, monospace">diff_6th_factor = 0.12, 0.12, 0.12, </font></div><div><font face="courier new, monospace">base_temp = 290.0,</font></div>
<div><font face="courier new, monospace">damp_opt = 0,</font></div><div><font face="courier new, monospace">zdamp = 5000.0, 5000.0, 5000.0, </font></div><div><font face="courier new, monospace">dampcoef = 0.2, 0.2, 0.2, </font></div>
<div><font face="courier new, monospace">khdif = 0, 0, 0, </font></div><div><font face="courier new, monospace">kvdif = 0, 0, 0, </font></div><div><font face="courier new, monospace">non_hydrostatic = .true., .true., .true., </font></div>
<div><font face="courier new, monospace">chem_adv_opt = 2, 2, 2, </font></div><div><font face="courier new, monospace">moist_adv_opt = 2, 2, 2, </font></div><div><font face="courier new, monospace">scalar_adv_opt = 2, 2, 2, </font></div>
<div><font face="courier new, monospace">tke_adv_opt = 2, 2, 2, </font></div><div><font face="courier new, monospace">/</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&bdy_control</font></div>
<div><font face="courier new, monospace">spec_bdy_width = 5,</font></div><div><font face="courier new, monospace">spec_zone = 1,</font></div><div><font face="courier new, monospace">relax_zone = 4,</font></div>
<div><font face="courier new, monospace">specified = .true., .false., .false., </font></div><div><font face="courier new, monospace">nested = .false., .true., .true., </font></div>
<div><font face="courier new, monospace">/</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&grib2</font></div><div><font face="courier new, monospace">/</font></div>
<div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&namelist_quilt</font></div><div><font face="courier new, monospace">nio_tasks_per_group = 0,</font></div><div>
<font face="courier new, monospace">nio_groups = 1,</font></div><div><font face="courier new, monospace">/</font></div><div><font face="courier new, monospace"><br></font></div><div><font face="courier new, monospace">&chem</font></div>
<div><font face="courier new, monospace">kemit = 23,</font></div><div><font face="courier new, monospace">chem_opt = 108, 108, 108, </font></div><div><font face="courier new, monospace">bioemdt = 30, 30, 30, </font></div>
<div><font face="courier new, monospace">photdt = 30, 30, 30, </font></div><div><font face="courier new, monospace">chemdt = 7.5, 2.5, 0.83, </font></div><div><font face="courier new, monospace">emiss_inpt_opt = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">emiss_opt = 3, 3, 3, </font></div><div><font face="courier new, monospace">io_style_emissions = 2,</font></div><div><font face="courier new, monospace">chem_in_opt = 0, 0, 0, </font></div>
<div><font face="courier new, monospace">phot_opt = 2, 2, 2, </font></div><div><font face="courier new, monospace">gas_drydep_opt = 1, 1, 1, </font></div><div><font face="courier new, monospace">aer_drydep_opt = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">bio_emiss_opt = 3, 3, 3, </font></div><div><font face="courier new, monospace">dust_opt = 0,</font></div><div><font face="courier new, monospace">dmsemis_opt = 0,</font></div>
<div><font face="courier new, monospace">seas_opt = 0,</font></div><div><font face="courier new, monospace">gas_bc_opt = 1, 1, 1, </font></div><div><font face="courier new, monospace">gas_ic_opt = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">aer_bc_opt = 1, 1, 1, </font></div><div><font face="courier new, monospace">aer_ic_opt = 1, 1, 1, </font></div><div><font face="courier new, monospace">gaschem_onoff = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">aerchem_onoff = 1, 1, 1, </font></div><div><font face="courier new, monospace">wetscav_onoff = 0, 0, 0, </font></div><div><font face="courier new, monospace">cldchem_onoff = 0, 0, 0, </font></div>
<div><font face="courier new, monospace">vertmix_onoff = 1, 1, 1, </font></div><div><font face="courier new, monospace">chem_conv_tr = 1, 1, 1, </font></div><div><font face="courier new, monospace">biomass_burn_opt = 1, 1, 1, </font></div>
<div><font face="courier new, monospace">plumerisefire_frq = 30, 30, 30, </font></div><div><font face="courier new, monospace">have_bcs_chem = .true., .true., .true., </font></div><div><font face="courier new, monospace">aer_ra_feedback = 0, 0, 0, </font></div>
<div><font face="courier new, monospace">/</font></div><div><br></div><div>Thank you for your time,</div></div></div><div><br></div>-- <br>Chris Tessum<br><br>Doctoral Candidate<br>Environmental Engineering<br>University of Minnesota<br>
+1 (612) 568 2238<br><a href="mailto:tess0050@umn.edu" target="_blank">tess0050@umn.edu</a>
</div>