Re: Fwd: Problem in lightning parameterization
Mary Barth
Hello Madhuparna,
With a 1-km domain, you are no longer using the convective
parameterization, which is good. When you turn off the convective
parameterization, you also need to change lightning options.
lightning_option = 1 or 2 work with cloud-resolving domains. I
would recommend trying the option that uses the cloud top height.
Cellcount_method should also be set to a value (likely the one that counts grid cells in the processor, option 1).
I suggest you use iccg_method = 4 and then set iccg_prescribed_num = 3 and iccg_prescribed_den = 1. This ratio can be varied, it depends on the region. (for example, look at Boccippio et al. (2001))
The flashrate_factor is to be used to adjust the number of flashes that are predicted. Lightning parameterizations still do not predict well compared to observations, so we rely on adjusting the parameterization to get the model representing observations.
-- Mary
Dear Madam,
I wrote you about my problem regarding lightning parameterization one year before. In between I was involved in other objectives and could not follow up in this matter.At present, my post doctoral research aims to predict the lightning events in India. I tried lightning parameterization freshly using WRF3.7.1 with lightning option 11, and 1,2.. My model has 4 nested domains (27,9,3,1 km). I am getting model simulated flash in 27, 9 and 3 km but in innermost domain fails to produce any flash.Can you please suggest about the solution of the error. Model is running successfully and producing normal wrf outputs in each domain. Only problem is in innermost domain and lightning parameterization. We are eager to get model predicted flash in innermost domain. My Namelist is as follows:&time_control
run_days = 0,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2017, 2017, 2017, 2017,
start_month = 05, 05, 05, 05,
start_day = 07, 07, 07, 07,
start_hour = 00,00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2017, 2017, 2017, 2017,
end_month = 05, 05, 05, 05,
end_day = 07, 07, 07, 07,
end_hour = 12,12, 12, 12,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 10800,
input_from_file = .true.,.true.,.true.,.true.,.true.,
history_interval = 30, 30, 30, 30, 02,
frames_per_outfile = 1000, 1000, 1000, 1000, 1000,
restart = .false.,
restart_interval = 50000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/
&domains
time_step =60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 4,
e_we = 172, 397, 802, 1171, 970,
e_sn = 148, 310, 622, 1030, 970,
e_vert = 45, 45, 45, 45, 52,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000, 1000, 333.3333,
dy = 27000, 9000, 3000, 1000, 333.3333,
grid_id = 1, 2, 3, 4, 5,
parent_id = 0, 1, 2, 3, 4,
i_parent_start = 1, 20, 70, 195, 300,
j_parent_start = 1, 21, 50, 170, 165,
parent_grid_ratio = 1, 3, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3, 3,
feedback = 1,
smooth_option = 0
/
/
&physics
mp_physics = 8, 8, 8, 8, 10,
ra_lw_physics = 1, 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1, 1,
radt = 27, 27, 27, 27, 27,
sf_sfclay_physics = 1, 1, 1, 1, 1,
sf_surface_physics = 2, 2, 2, 2, 2,
bl_pbl_physics = 1, 1, 1, 1, 1,
bldt = 0, 0, 0, 0, 0,
cu_physics = 93, 93, 93, 93, 0,
cudt = 5, 5, 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 20,
sf_urban_physics = 0, 0, 0, 0, 0,
do_radar_ref = 1,
lightning_option = 11, 11, 11, 11, 11,
lightning_start_seconds = 600,600,600,600,
lightning_dt = 60, 60, 60, 60,
flashrate_factor =1.0,0.95,0.9,0.8,
cellcount_method =0,0,0,0,
cldtop_adjustment = 2, 2, 2, 2,
iccg_method = 2, 2, 2, 2,
iccg_prescribed_num = 0.,
iccg_prescribed_den = 1.,
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1, 1, 1, 1, 1,
km_opt = 4, 4, 4, 4, 4,
diff_6th_opt = 0, 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.2, 0.2,
khdif = 0, 0, 0, 0, 0,
kvdif = 0, 0, 0, 0, 0,
non_hydrostatic = .true., .true., .true., .true., .true.,
moist_adv_opt = 1, 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,.false., .false.,
nested = .false., .true., .true., .true., .true.,
/
&grib2
/
&namelist_quiltnio_tasks_per_group = 0,
nio_groups = 1,
/
I have changed lightning_dt=60,20,20,20 and flash rate factors to 1.0 everywhere and other combinations also. But nothing is working for 1km domain. What may be the cause of this type of error madam?Please help in this regard.Thanking You,With Regards,Madhuparna
On Fri, Feb 24, 2017 at 5:25 AM, Mary Barth <bar...@ucar.edu> wrote:
Hello Madhuparna,
Is there a place where you can put an example output so that I can see what is happening?
Mary
On 2/23/17 11:56 AM, Madhuparna Halder wrote:
Dear Madam,Thank You for your help. But I am still unable to produce any flash with lightning parameterization (PR92). WRF3.6.1 could not resolve the problem of lightning_dt with lightning_dt = 60, 20, 20, 20.It was giving the same error. Then I installed WRF 3.8.1. With lightning_dt= 60,20,20,20 , I ran the model and there was no error as such but unfortunately it could not simulate any flash in any time step though it was experimented for a lightning event.What could be the problem and how to rectify the same ? Please suggest in this regard.Thanking You,With Regards,Madhuparna
On Thu, Feb 16, 2017 at 9:51 PM, Mary Barth <bar...@ucar.edu> wrote:
Hello Madhuparna,
V3.6.1 has a couple of bugs that we fixed for V3.7 and later. One of them is the lightning_dt. However, since you give the lightning_dt in the namelist, I do not think it was the bug but may be because the "model time step" calculated from the dt at 27 km domain is different than the lightning_dt set in the namelist. For example, the 27-km domain has dt=60s, resulting in the 9-km domain having dt=20s, and the 3-km domain having dt=6.666667s. And the lightning_dt for the 9-km domain is 6.66s which may not give a value of 0 when the mod(lightning_dt, time_step) is checked.
I would suggest trying: lightning_dt = 60, 20, 20, 20
Regards,
Mary
On 2/16/17 2:08 AM, Madhuparna Halder wrote:
Dear Madam,,I am working on estimation and forecast of lightning over India using WRF. I am facing problem in running WRF 3.6.1 with lightning option. I have read your paper in Atmospheric Research.
I am getting error "FATAL CALLED FROM FILE: <stdin> LINE: 87lightning_init: lightning_dt needs to be a multiple of model time step dt".Do you have faced similar kind of error ?
Following is the namelist.input:&time_control
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2015, 2015, 2015, 2015,
start_month = 03, 03, 03, 03,
start_day = 14, 14, 14, 14,
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2015, 2015, 2015, 2015,
end_month = 03, 03, 03, 03,
end_day = 15, 15, 15, 15,
end_hour = 00, 00, 00, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
history_interval = 720, 720, 720, 30,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,.true.,
frames_per_outfile = 1000, 1000, 1000, 1000,
restart = .false.,
restart_interval = 180,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput4 = 2
debug_level = 0
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 1440
/
&domains
time_step = 60
time_step_fract_num = 00
max_dom = 4,
e_we = 91, 148, 262, 391,
e_sn = 85, 160, 262, 391,
e_vert = 31, 31, 31,31,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000, 1000,
dy = 27000, 9000, 3000, 1000,
grid_id = 1, 2, 3, 4,
parent_id = 1, 1, 2, 3,
i_parent_start = 1, 22, 30, 67,
j_parent_start = 1, 19, 30, 67
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 10, 10, 10, 10,
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 30, 30, 30, 30,
sf_sfclay_physics = 1, 1, 1, 1,
sf_surface_physics = 2, 2, 2, 2,
bl_pbl_physics = 1, 1, 1, 1,
bldt = 0, 0, 0, 0,
cu_physics = 93, 93, 0, 0,
cudt = 5, 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
do_radar_ref = 1,
lightning_option = 11,1,1,1,
lightning_start_seconds =600,600,600,600,
lightning_dt = 60, 20, 6.66, 2.22,
flashrate_factor =1.0,1.0,1.0,1.0,
cellcount_method =0,0,0,0,
cldtop_adjustment = 2, 2, 2, 2,
iccg_method = 2, 2, 2, 2,
iccg_prescribed_num = 0.,
iccg_prescribed_den = 1.,
/
&fdda
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.2,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
non_hydrostatic = .true., .true., .true., .true,
moist_adv_opt = 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,.false.,
nested = .false., .true., .true.,.true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
What changes are required in lightning option ? How to modify lightning_dt ?
Please suggest regarding solving the problem. I will be grateful if you kindly help me as I am stuck in middle of the research problem.
Thanking You,
With Regards,
Dr. Madhuparna HalderProject ScientistIndian Institute of Tropical MeteorologyDr. Homi Bhabha Road,PashanPune-411008
--
==========================================With Regards,Dr. Madhuparna HalderPost Doctoral FellowIndian Institute of Tropical MeteorologyDr. Homi Bhabha Road,PashanPune-411008, Maharashtra, India
--
==========================================With Regards,Dr. Madhuparna HalderPost Doctoral FellowIndian Institute of Tropical MeteorologyDr. Homi Bhabha Road,PashanPune-411008
--
==========================================With Regards,Dr. Madhuparna HalderPost Doctoral FellowIndian Institute of Tropical MeteorologyDr. Homi Bhabha Road,PashanPune-411008
-- ^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^ 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/ +--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+
--
==========================================With Regards,Dr. Madhuparna HalderPost Doctoral FellowAllahabad UniversityAllahabad-211002
-- ^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^--^ 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/ +--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+