Page 1 of 1

How to calculate lighting and LPI using wrf3.6.1?

PostPosted: Wed Sep 24, 2014 6:04 am
by lynn831018
hi everyone,

Recently I made a numeric model test about a thunderstorm for calculating LPI index (it's a lightning potential index) using WRFV3.6.1. But the result is very poor. The LPI is zero in second and third grid except mother grid. In my opinion, thunderstorm happened in third grid, so LPI should not equal to zero.I think the lightning parameterization options may be wrong in “namelist.input” file . And the namelist.input setting as follows:
&time_control
run_days = 0,
run_hours = 25,
run_minutes = 0,
run_seconds = 0,
start_year = 2008, 2008, 2008, 2008,
start_month = 06, 06, 06, 06,
start_day = 29, 29, 29, 29,
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2008, 2008, 2008, 2008,
end_month = 06, 06, 06, 06,
end_day = 30, 30, 30, 30,
end_hour = 00, 00, 00, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,.true.,
history_interval = 60, 60, 60, 60,
frames_per_outfile = 1000, 1000, 1000, 1000,
restart = .false.,
restart_interval = 60,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/

&domains
time_step = 81,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
s_we = 1, 1, 1, 1,
e_we = 100, 169, 211, 235,
s_sn = 1, 1, 1, 1,
e_sn = 100, 169, 211, 238,
s_vert = 1, 1, 1, 1,
e_vert = 38, 38, 38,38,
num_metgrid_levels = 27,
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, 23, 50, 68,
j_parent_start = 1, 23, 50, 64,
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 1,
smooth_option = 0,
num_metgrid_soil_levels = 4,
force_sfc_in_vinterp = 3,
zap_close_levels = 500
sfcp_to_sfcp = .false.
adjust_heights = .false.
eta_levels = 1.000, 0.995, 0.990, 0.985,
0.980,
0.970, 0.960, 0.950, 0.940,
0.930,
0.920, 0.910, 0.900, 0.880,
0.860,
0.830, 0.800, 0.770, 0.740,
0.710, 0.680,
0.640, 0.600, 0.560, 0.520,
0.480,
0.440, 0.400, 0.360, 0.320,
0.280,
0.240, 0.200, 0.160, 0.120,
0.080,0.040,0


&physics
mp_physics = 8, 8, 8, 8,
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 15, 15, 15, 15,
sf_sfclay_physics = 2, 2, 2, 2,
sf_surface_physics = 2, 2, 2, 2,
bl_pbl_physics = 2, 2, 2, 2,
bldt = 0, 0, 0, 0,
cu_physics = 1, 1, 0, 0,
cudt = 5, 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
sf_urban_physics = 0, 0, 0, 0,
topo_wind = 0, 0, 1, 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,


do_radar_ref = 1,
lightning_option = 3,3,3,3,
lightning_dt = 81,
lightning_start_seconds = 0., 0.,0.,0.,
flashrate_factor = 3.5, 3.5,3.5,3.5,
cellcount_method = 1, 1, 2, 2,
cldtop_adjustment = 2., 2., 2., 2.,
iccg_method = 2, 2, 2, 2,
iccg_prescribed_num = 0.,
iccg_prescribed_den = 1.,
/

&fdda
grid_fdda = 0, 0, 0, 0,
fgdt = 0, 0, 0, 0,
if_no_pbl_nudging_uv = 0, 0, 0, 1,
if_no_pbl_nudging_t = 0, 0, 0, 1,
if_no_pbl_nudging_q = 0, 0, 0, 1,
if_zfac_uv = 0, 0, 0, 1,
k_zfac_uv = 10, 10, 10, 1,
if_zfac_t = 0, 0, 0, 1,
k_zfac_t = 10, 10, 10, 1,
if_zfac_q = 0, 0, 0, 1,
k_zfac_q = 10, 10, 10, 1,
guv = 0.0003, 0.0003, 0.0003, 0.0003,
gt = 0.0003, 0.0003, 0.0003, 0.0003,
gq = 0.0003, 0.0003, 0.0003, 0.0003,
if_ramping = 0,
dtramp_min = 120.0
/

&dynamics
dyn_opt = 2,
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.0, 0.0, 0.0, 0.0,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
smdiv = 0.1, 0.1, 0.1, 0.1,
emdiv = 0.01, 0.01, 0.01, 0.01,
epssm = 0.1, 0.1, 0.1, 0.1,
non_hydrostatic = .true., .true., .true., .true.,
time_step_sound = 4, 4, 4, 4,
h_mom_adv_order = 5, 5, 5, 5,
v_mom_adv_order = 3, 3, 3, 3,
h_sca_adv_order = 3, 5, 5, 5,
v_sca_adv_order = 2, 3, 3, 3,
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.,
periodic_x = .false.,.false.,.false.,.false.,
symmetric_xs = .false.,.false.,.false.,.false.,
symmetric_xe = .false.,.false.,.false.,.false.,
open_xs = .false.,.false.,.false.,.false.,
open_xe = .false.,.false.,.false.,.false.,
periodic_y = .false.,.false.,.false.,.false.,
symmetric_ys = .false.,.false.,.false.,.false.,
symmetric_ye = .false.,.false.,.false.,.false.,
open_ys = .false.,.false.,.false.,.false.,
open_ye = .false.,.false.,.false.,.false.,
nested = .false., .true., .true.,.true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,

Could you help me to find the problem? How can I do to get more reasonable result? Thank you very much!

Re: How to calculate lighting and LPI using wrf3.6.1?

PostPosted: Fri Jul 21, 2017 12:34 pm
by jaymtee
Did you ever get a resolution to this problem?

Re: How to calculate lighting and LPI using wrf3.9.1?

PostPosted: Sun Nov 18, 2018 1:06 am
by darshana
How to get value for LPI ?