How to calculate lighting and LPI using wrf3.6.1?

Topics related to current and future physics in the WRF as well as any problems you may have.

How to calculate lighting and LPI using wrf3.6.1?

Postby lynn831018 » Wed Sep 24, 2014 6:04 am

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!
lynn831018
 
Posts: 4
Joined: Tue Dec 04, 2012 5:41 am

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

Postby jaymtee » Fri Jul 21, 2017 12:34 pm

Did you ever get a resolution to this problem?
jaymtee
 
Posts: 5
Joined: Tue Aug 23, 2016 11:20 am


Return to Physics

Who is online

Users browsing this forum: No registered users and 8 guests