Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

No vertical co-ordinate data in WRF-CHEM simulation

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

madhavhmk

New member
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,
/
 
Hi,

What do you mean by this:
"After the run, however, there is no data in vertical levels of the wrfout.
Only surface data is available."

Do you mean that there is no chemical data above the surface (o3, no2, no, etc). Or do you mean that there is no emissions data above the surface (e_no, e_so2, etc)?

How are you providing emissions data to the model? Do you have 1 file per hour, and is the emissions file naming convention:
wrfchemi_d01_2016-11-16_00:00:00
wrfchemi_d01_2016-11-16_01:00:00
wrfchemi_d01_2016-11-16_02:00:00
wrfchemi_d01_2016-11-16_03:00:00
etc?

Or are you organising your emission files differently?

On general namelist issues (which might not be directly related to your problem):
1) you shouldn't set "run_hours" if you have set the start_* and end_* variables. So I recommend changing this to zero.
2) you should make sure that your chemistry time-steps are integer multiples of the WRF time_step, otherwise you'll have issues with model timings. I'd recommend trying time_step = 30 (at least until you've got everything working).
 
Thankyou very much for your reply.......

By no data, I am refering to NO2 in wrfout file, where I am not having any values in vertical levels.
Temporally also, the evolution is dropping to zero and not realistic.

I am providing the emission files as mentioned:
wrfchemi_d01_2016-11-16_00:00:00
wrfchemi_d01_2016-11-16_01:00:00
wrfchemi_d01_2016-11-16_02:00:00
wrfchemi_d01_2016-11-16_03:00:00 etc

But I am finding that this is not hte case with other species like o3 or co. They are having representative values both temporally and vertically.

So I guess my problem is something to do with the way NO2 is mapped from the inventory to the wrfchemi input file.
Any guesses for not having a proper NO2 ?

Many thanks.
 
It does sound like you haven't any NOx emissions.

What tool are you using to create your wrfchemi* emission files? And have you checked the variables within your wrfchemi* files using ncview (or a similar tool), to make sure there are values in the E_NO or E_NO2 arrays?
 
I am using anthro_emis tool for creating wrfchemi. Here nox emissions are read from EDGAR-HTAP. We can set the mapping manually for NO2.
And the wrfchemi files seems to have e_no2 in them properly with all the proper hotspots in place etc.......

I am checking the emission mapping once again and I will see....


Thanks,
 
If you would like us to check your wrfchemi* and wrfout* files for you too (so we can check them against other simulations, to see if they look reasonable) - then you could upload some example files to the nextcloud account following the instructions on the front page of the forum?
 
Top