Hello all,
I am able to successfully run the WRF-CHEM with MOZART data (Chem_opt=112, emiss_opt=8).
After the run, however, there is no data in vertical levels of the wrfout.
Only surface data is available.
I have set the kemit and vertical levels in the namelist.
Also, the temporal evolution of the concentrations are not captured.
Only the initial first housr data is good. The rest is falling to zero, even though I set the run for six hours.
I am using 3.8.1 version of the model. Any help regarding this is highly appreciated..!!!
----------------------------Namelist.input--------------------------------------------------------------
&time_control
run_days = 0,
run_hours = 6,
run_minutes = 0,
run_seconds = 0,
start_year = 2016, 2017, 2017,
start_month = 11, 08, 08,
start_day = 16, 15, 15,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2016, 2017, 2017,
end_month = 11, 08, 08,
end_day = 16, 15, 15,
end_hour = 06, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
frames_per_auxinput5 = 1
io_form_input = 2
io_form_boundary = 2
io_form_auxinput5 = 2
io_form_auxinput6 = 0
io_form_auxinput7 = 0
io_form_auxinput8 = 0
debug_level = 00
auxinput5_interval_m = 60, 3600, 3600
/
&domains
time_step = 36,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 85, 112, 94,
e_sn = 35, 97, 91,
e_vert = 30, 30, 30,
dx = 5000, 200, 66.66,
dy = 5000, 200, 66.66,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 30, 30,
j_parent_start = 1, 20, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 2, 2, 2,
progn = 0, 0, 0,
naer = 1e9
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 2, 2, 2,
radt = 30, 10, 10,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 1, 0, 0,
cu_physics = 5, 5, 0,
cudt = 1, 1, 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&chem
kemit = 10,
chem_opt = 112, 0,
chem_in_opt = 0,
bioemdt = 30, 30,
photdt = 30, 30,
chemdt = 2., 2.,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1,
emiss_opt = 8, 2,
phot_opt = 0, 1,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 1,
bio_emiss_opt = 0, 0,
ne_area = 100
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
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 = 0, 0,
cldchem_onoff = 0, 0,
vertmix_onoff = 1, 1,
chem_conv_tr = 0, 0,
biomass_burn_opt = 2, 0,
plumerisefire_frq = 30, 30,
aer_ra_feedback = 0, 0,
have_bcs_chem = .false., .false.,
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am able to successfully run the WRF-CHEM with MOZART data (Chem_opt=112, emiss_opt=8).
After the run, however, there is no data in vertical levels of the wrfout.
Only surface data is available.
I have set the kemit and vertical levels in the namelist.
Also, the temporal evolution of the concentrations are not captured.
Only the initial first housr data is good. The rest is falling to zero, even though I set the run for six hours.
I am using 3.8.1 version of the model. Any help regarding this is highly appreciated..!!!
----------------------------Namelist.input--------------------------------------------------------------
&time_control
run_days = 0,
run_hours = 6,
run_minutes = 0,
run_seconds = 0,
start_year = 2016, 2017, 2017,
start_month = 11, 08, 08,
start_day = 16, 15, 15,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2016, 2017, 2017,
end_month = 11, 08, 08,
end_day = 16, 15, 15,
end_hour = 06, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
frames_per_auxinput5 = 1
io_form_input = 2
io_form_boundary = 2
io_form_auxinput5 = 2
io_form_auxinput6 = 0
io_form_auxinput7 = 0
io_form_auxinput8 = 0
debug_level = 00
auxinput5_interval_m = 60, 3600, 3600
/
&domains
time_step = 36,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 85, 112, 94,
e_sn = 35, 97, 91,
e_vert = 30, 30, 30,
dx = 5000, 200, 66.66,
dy = 5000, 200, 66.66,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 30, 30,
j_parent_start = 1, 20, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 2, 2, 2,
progn = 0, 0, 0,
naer = 1e9
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 2, 2, 2,
radt = 30, 10, 10,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
bldt = 1, 0, 0,
cu_physics = 5, 5, 0,
cudt = 1, 1, 1,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&chem
kemit = 10,
chem_opt = 112, 0,
chem_in_opt = 0,
bioemdt = 30, 30,
photdt = 30, 30,
chemdt = 2., 2.,
io_style_emissions = 2,
emiss_inpt_opt = 1, 1,
emiss_opt = 8, 2,
phot_opt = 0, 1,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 1,
bio_emiss_opt = 0, 0,
ne_area = 100
dust_opt = 0,
dmsemis_opt = 0,
seas_opt = 0,
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 = 0, 0,
cldchem_onoff = 0, 0,
vertmix_onoff = 1, 1,
chem_conv_tr = 0, 0,
biomass_burn_opt = 2, 0,
plumerisefire_frq = 30, 30,
aer_ra_feedback = 0, 0,
have_bcs_chem = .false., .false.,
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/