WRF-Chem V3.8.1 lightning no emissions problem

178 views
Skip to first unread message

dominique fonteyn

unread,
Aug 8, 2017, 5:59:38 AM8/8/17
to wrf-chem-lightning
Hello,

I'm trying to run WRF_Chem (MOZCART) for the August-September 2013 period over the US (seac4rs campaign), but i seem to have a bit of a problem with the NO lightning emissions, there were none.

At first I thought the namelsit options were problematic (see below)

But trying to understandmore I had a look at the code, and had look at solve_em.F (as mentioned in the user guide) .

The use of the module_lightning_driver has been commented out and there is no call to the lightning driver present in the V3.8.1 version of solve_em.F

I presume this has a major impact on the lightning emissions.

Any comments/suggestions on how to solve the problem.

Thanks very much

 &physics
 mp_physics                          = 10,
 progn                               = 1,
 ra_lw_physics                       = 4,
 ra_sw_physics                       = 4,
 radt                                = 30,
 sf_sfclay_physics                   = 1,
 sf_surface_physics                  = 2,
 bl_pbl_physics                      = 1,
 bldt                                = 0,
 cu_physics                          = 5,
 cu_diag                             = 1,
 cudt                                = 0,
 ishallow                            = 0,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 1,
 mp_zero_out                         = 2,
 mp_zero_out_thresh                  = 1.e-8
 maxiens                             = 1,
 maxens                              = 3,
 maxens2                             = 3,
 maxens3                             = 16,
 ensdim                              = 144,
 cu_rad_feedback                     = .true.,
 num_land_cat                        = 21,
 do_radar_ref                        = 1
 lightning_option                    = 11,
 iccg_method                         = 2,
 lightning_start_seconds             = 600,
 cldtop_adjustment                   = 2,
 flashrate_factor                    = 20,

&chem
 kemit                               = 1,
 chem_opt                            = 112,
 bioemdt                             = 30.0,
 photdt                              = 30.0,
 chemdt                              = 10,
 io_style_emissions                  = 2,
 emiss_inpt_opt                      = 111,
 emiss_opt                           = 8,    1,
 emiss_opt_vol                       = 0,
 chem_in_opt                         = 1,
 phot_opt                            = 3,
 gas_drydep_opt                      = 1,
 aer_drydep_opt                      = 1,
 bio_emiss_opt                       = 0,
 ne_area                             = 354,
 dust_opt                            = 1,
 dmsemis_opt                         = 0,
 seas_opt                            = 1,
 gas_bc_opt                          = 1,
 gas_ic_opt                          = 1,
 aer_bc_opt                          = 1,
 aer_ic_opt                          = 1,
 gaschem_onoff                       = 1,
 aerchem_onoff                       = 1,
 wetscav_onoff                       = 1,
 cldchem_onoff                       = 0,
 vertmix_onoff                       = 1,
 chem_conv_tr                        = 1,
 conv_tr_wetscav                     = 1,
 conv_tr_aqchem                      = 1,
 biomass_burn_opt                    = 2,
 aer_op_opt                          = 1,
 aer_ra_feedback                     = 1,
 opt_pars_out                        = 0,
 diagnostic_chem                     = 0,
 scale_fire_emiss                    = .true.,
 have_bcs_chem                        = .true.,
 have_bcs_upper                      = .true.,
 fixed_ubc_inname                    = ubvals_b40.20th.track1_1996-2005.nc,
 lightning_opt                       = 101,
 lightning_time_step                 = 1,
 

Stacy Walters

unread,
Aug 8, 2017, 8:29:43 AM8/8/17
to dominique fonteyn, wrf-chem-lightning
Dominic,

Please use the discussion forum at:

https://www2.acom.ucar.edu/wrf-chem/discussion-forum

Most likely you will need to register before you use a specific forum.

Stacy

bar...@ucar.edu

unread,
Aug 8, 2017, 11:01:59 AM8/8/17
to wrf-chem-lightning
Hello Dominique,

The lightning flash rate routine is now called from phys/module_diagnostics_driver.F which is called from dyn_em/module_after_all_rk_steps.F.

It seems to me that the namelist is set up okay for getting lightning flashrates. Can you verify that? There should be a 2D array called flashcount or flashrate in the wrfout file. The only issue I have had in the physics part is making sure the cu_physics is set to one of the Grell schemes, which is what you have. As a note, if you change it to the Grell-Freitas scheme, the program will stop because we haven't yet altered the lightning_init subroutine to allow for the GF scheme (sorry!), but you can easily alter that code.

The lightning-NOx namelist options do not look right. It should be lnox_option = 1 or 2 and there's no lightning_time_step here. Check the User's Guide again. To help, I have these as my namelist options in the &chem section:
 lnox_opt                            = 1,
 N_IC                                = 125.,
 N_CG                                = 125.,

Let us know how it goes.
Mary

dominique fonteyn

unread,
Aug 9, 2017, 5:39:39 AM8/9/17
to wrf-chem-lightning
Hello Mary,

Now It seems to be working fine.
Thanks very much, it is all very sensitive to namelist input

again thanks for your kind support

Op dinsdag 8 augustus 2017 17:01:59 UTC+2 schreef bar...@ucar.edu:
Reply all
Reply to author
Forward
0 new messages