Simulate lightning flashes and NOx
Arshini Saikia
Mary Barth
Dear Arshini,
It looks like the iccg_method should be set differently. Option 2 is specific to the United States. For the Eastern Himalayan region, I recommend using Option 1 and set the IC:CG value via iccg_prescribed_num (numerator) and iccg_prescribed_den (denominator). My recommendation would be:
iccg_method = 1
iccg_prescribed_num = 2
iccg_prescribed_den = 1
Do you have any information on the IC to CG ratio for your region? If so, that would be more reliable. If not, then my suggested value of IC/CG = 2 is fairly conservative. The Boccippio et al. (2001) paper shows a map for the United States with IC/CG values varying from 1 to 8.
Good luck!
Mary Barth
Mary Barth
Hello Arshini,
First, I noticed that you are using one domain (max_dom=1) and
dx=15km with dt=90seconds. Note that WRF will read only the first
column of your namelist settings because max_dom=1.
lightning_option = 11,
iccg_method = 1,
iccg_prescribed_num = 2,
iccg_prescribed_den = 1,
lightning_start_seconds = 600,
cldtop_adjustment = 2,
flashrate_factor = 1.0,
lightning_dt = 90,
cellcount_method is not used by the lightning_option = 11. You should either not include it or set it to 0.
The lightning NOx settings that you have look good. Just to repeat them, they are set to:
lnox_opt = 1, 1,
N_IC = 125., 125.,
N_CG = 125., 125.,
-- Mary
Dear Mam,
As per your suggestion I re-run the model but in wrfoutput i got the IC flashcount and CG flashout value = 0, even lnox_total value is also zero. I attached the edited namelist below. Please give your valuable suggestion.
with regardsArshini
--
You received this message because you are subscribed to the Google Groups "wrf-chem-lightning" group.
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem-lightning+unsub...@ucar.edu.
-- ^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^ Mary Barth Phone: 303-497-8186 Senior Scientist email: bar...@ucar.edu National Center for Atmospheric Research P.O. Box 3000 Boulder, CO 80307 https://staff.ucar.edu/users/barthm DC3 Web Site: http://www2.acom.ucar.edu/dc3/ +--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+
Mary Barth
Arshini,
Have you looked at RAINC to see if it is greater than 0? If RAINC
= 0, then there's no convective rain and therefore there would not
be lightning. If RAINC>0, it may be that ktop (level of neutral
buoyancy) is not reaching heights high enough to produce
lightning. Do you have KTOP in the output file? If not, you can
add it to wrfout (via iofields) and do another run (please make
sure I have the correct name/spelling).
You could also set cldtop_adjustment = 0 to give higher values of flashrate. (see Wong et al., 2013 for what effect cldtop_adjustment has on the parameterization)
-- Mary
Dear Mam,
A very wealthy and healthy New Year to you.
Sorry, still my model run is showing IC_FLASHCOUNT , CG_FLASHCOUNT and lnox_total value zero. Should i try to run my namelist for America, to check weather my namelist is working or not or have some issues. Please share your valuable suggestion.with regardsArshini
To unsubscribe from this group and stop receiving emails from it, send an email to wrf-chem-lightning+unsubscribe@ucar.edu.
-- ^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^ Mary Barth Phone: 303-497-8186 Senior Scientist email: bar...@ucar.edu National Center for Atmospheric Research P.O. Box 3000 Boulder, CO 80307 https://staff.ucar.edu/users/barthm DC3 Web Site: http://www2.acom.ucar.edu/dc3/ +--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+