No anthropogenic emissions in the wrfout file

Any issues with running WRF Chemistry

No anthropogenic emissions in the wrfout file

Postby onkeymodec » Mon Apr 20, 2015 1:58 am

I compiled WRF-CHEM V3.6.1 and performed exercises in the tutorial.
These tutorials are provided on WRF-CHEM's website.
Things went on well until the Exercise 2, I checked the wrfinput file but there is no items of anthropogenic emissions (for instance, the E_CO item). But I CHECKED all the files required for namelist.input and confirm that auxiliary input files are correct using NCVIEW.
As a result, in the wrfout file, there is no items like PM2_5_DRY.Here is the namelist.input file:
Code: Select all
 &time_control
 run_days                            = 2,
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2010, 2010, 2000,
 start_month                         = 07,   07,   01,
 start_day                           = 14,   14,   24,
 start_hour                          = 00,   00,   12,
 start_minute                        = 00,   00,   00,
 start_second                        = 00,   00,   00,
 end_year                            = 2010, 2010, 2000,
 end_month                           = 07,   07,   01,
 end_day                             = 16,   16,   25,
 end_hour                            = 00,   00,   12,
 end_minute                          = 00,   00,   00,
 end_second                          = 00,   00,   00,
 interval_seconds                    = 21600,
 input_from_file                     = .true.,.true.,.false.,
 history_interval                    = 60,  60,   60,
 frames_per_outfile                  = 72, 1000, 1000,
 restart                             = .false.,
 restart_interval                    =    0,
 io_form_history                     = 2,
 io_form_restart                     = 2,
 io_form_input                       = 2,
 io_form_boundary                    = 2,
 auxinput5_inname                    = 'wrfchemi_00z_d<domain>'
 auxinput6_inname                    = 'wrfbiochemi_d01',
 auxinput7_inname                    = 'wrffirechemi_d<domain>',
 auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
 auxinput12_inname                   = 'wrf_chem_input',
 auxinput13_inname                   = 'wrfchemv_d<domain>',
 auxinput5_interval_m                = 86400, 86400, 60,
 auxinput7_interval_m                = 86400, 86400, 60,
 auxinput8_interval_m                = 86400, 86400, 60,
 auxinput13_interval_m               = 86400, 86400, 60,
 io_form_auxinput2                   = 2,
 io_form_auxinput5                   = 2,
 io_form_auxinput6                   = 0,
 io_form_auxinput7                   = 2,
 io_form_auxinput8                   = 2,
 io_form_auxinput12                  = 0,
 io_form_auxinput13                  = 0,
 debug_level                         = 0,
 auxinput1_inname                    = "met_em.d<domain>.<date>",
/

 &dfi_control
/

 &domains
 time_step                           = 600,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 s_we                                =   1,    1,     1,
 e_we                                = 41,    41,    94,
 s_sn                                =   1,    1,     1,
 e_sn                                = 41,    41,    91,
 e_vert                              = 31,    31,    28,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 100000, 4000,  3333.33,
 dy                                  = 100000, 4000,  3333.33,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     6,    30,
 j_parent_start                      = 1,     65,    30,
 parent_grid_ratio                   = 1,     5,     3,
 parent_time_step_ratio              = 1,     5,     3,
 p_top_requested                     = 5000,
 feedback                            = 1,
 smooth_option                       = 0
 p_top_requested                     = 5000
 zap_close_levels                    = 50
 interp_type                         = 1
 t_extrap_type                       = 2
 force_sfc_in_vinterp                = 0
 use_levels_below_ground             = .true.
 use_surface                         = .true.
 lagrange_order                      = 1
 /
 sfcp_to_sfcp                        = .true.

 &physics
 mp_physics                          = 4,     2,     2,
 progn                               = 0,
 ra_lw_physics                       = 1,     1,     1,
 ra_sw_physics                       = 2,     1,     1,
 radt                                = 30,    30,    30,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 2,     2,     1,
 bl_pbl_physics                      = 1,     1,     1,
 bldt                                = 0,     0,     0,
 cu_physics                          = 5,     5,     0,
 cu_diag                             = 1,     0,     0,
 cudt                                = 0,     0,     0,
 ishallow                            = 0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,     0,     0,
 mp_zero_out                         = 2,
 mp_zero_out_thresh                  = 1.e-12
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 cu_rad_feedback                     = .true.,
 /

 &fdda
 /

 &dynamics
 rk_ord                              = 3,
 w_damping                           = 1,
 diff_opt                            = 1,
 km_opt                              = 4,
 diff_6th_opt                        = 0,
 diff_6th_factor                     = 0.12,
 base_temp                           = 290.
 damp_opt                            = 0,
 zdamp                               = 5000.,  5000.,  5000.,
 dampcoef                            = 0.01,   0.01,   0.01
 khdif                               = 0,      0,      0,
 kvdif                               = 0,      0,      0,
 non_hydrostatic                     = .true., .true., .true.,
 moist_adv_opt                       = 2,      0,      0,
 scalar_adv_opt                      = 2,      0,      0,
 chem_adv_opt                        = 2,      0,      0,
 tke_adv_opt                         = 2,      0,      0,
 time_step_sound                     = 4,      4,      4,
 h_mom_adv_order                     = 5,      5,      5,
 v_mom_adv_order                     = 3,      3,      3,
 h_sca_adv_order                     = 5,      5,      5,
 v_sca_adv_order                     = 3,      3,      3,
 /

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

 &chem
 kemit                               = 1,
 chem_opt                            = 300,
 emiss_inpt_opt                      = 1,
 bioemdt                             = 0,
 photdt                              = 0,
 chemdt                              = 60,
 io_style_emissions                  = 1, 
 emiss_opt                           = 5,          0,
 emiss_opt_vol                       = 0,          0,
 emiss_ash_hgt                       = 20000., 
 chem_in_opt                         = 0,          0,
 phot_opt                            = 0,          0,
 gas_drydep_opt                      = 1,          1,
 aer_drydep_opt                      = 1,          1,
 bio_emiss_opt                       = 1,          0,
 dust_opt                            = 1,
 dmsemis_opt                         = 1,
 seas_opt                            = 1,
 depo_fact                           = 0.25
 gas_bc_opt                          = 1,          0,
 gas_ic_opt                          = 1,          0,
 aer_bc_opt                          = 1,          0,
 aer_ic_opt                          = 1,          0,
 gaschem_onoff                       = 1,          0,
 aerchem_onoff                       = 1,          0,
 wetscav_onoff                       = 0,          0,
 cldchem_onoff                       = 0,          0,
 vertmix_onoff                       = 1,          0,
 chem_conv_tr                        = 1,          0,
 conv_tr_wetscav                     = 0,          0,
 conv_tr_aqchem                      = 0,          0,
 biomass_burn_opt                    = 1,          0,
 plumerisefire_frq                   = 30,         0,
 have_bcs_chem                       = .false., .false., .false.,
 aer_ra_feedback                     = 1,
 aer_op_opt                          = 1,
 opt_pars_out                        = 0,
 diagnostic_chem                     = 0,
 /



I suspect there are incorrect configurations in 'aux' part.
Thanks in advance if you have any idea solving the problem!
onkeymodec
 
Posts: 2
Joined: Thu Apr 16, 2015 3:13 am

Re: No anthropogenic emissions in the wrfout file

Postby ahwrf » Mon Apr 20, 2015 2:31 am

I do not see any special mistake in your namelist.

How did you generate the emission files? Have you used the program "prep_chem_sources"?

If yes, you should have made some separate emission files which should be linked in WRFV3/run or WRV3/test/em_real directories beside wrfinput, ... to be used in the run process.
ahwrf
 
Posts: 13
Joined: Thu Dec 15, 2011 3:36 am
Location: Persia

Re: No anthropogenic emissions in the wrfout file

Postby villamen » Mon Apr 25, 2016 6:06 pm

Hi, this is a late reply so you might have already found out by now that you are not supposed to getting any anthropogenic emissions in the exercise 2. Just dust, DMS, sea salt, and biomass burning emissions. So you probably did everything fine.
villamen
 
Posts: 8
Joined: Mon Apr 25, 2016 5:52 pm

Re: No anthropogenic emissions in the wrfout file

Postby villamen » Mon Apr 25, 2016 10:23 pm

Sorry, just realized I was looking at a different exercise. I tried to do this exercise 2 as well, and run into the same problem: I do not see anthropogenic emissions in my wrfinput file after running real.exe. If you were able to solve it, could you tell me how? Thanks!
villamen
 
Posts: 8
Joined: Mon Apr 25, 2016 5:52 pm

Re: No anthropogenic emissions in the wrfout file

Postby ahwrf » Tue May 17, 2016 7:29 am

Have you used the program "prep_chem_sources.exe" and global emission inventory? If so, you mustn't have any problem, unless because of namelist.input issues.
ahwrf
 
Posts: 13
Joined: Thu Dec 15, 2011 3:36 am
Location: Persia

Re: No anthropogenic emissions in the wrfout file

Postby hmacintyre » Fri Oct 21, 2016 11:25 am

Hi, this might be a bit late, but I had a similar issue. I had to set io_form_auxinput options to 2 when running real.exe (to get the emissions into tthe wrfinput file), but then turned them off (set io_form_auxinput = 0) to run wrf.exe. I didn't find that it made any difference when I left it set to 2 when running wrf.exe though, so I assume it just reads it in twice and overwrites. Hope this helps.
hmacintyre
 
Posts: 12
Joined: Thu Dec 04, 2014 4:31 pm


Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 2 guests