davegill
commented
4 years ago @bmalecic Given the suggestions that I initially had below, please move this over to the WRF Forum: forum.mmm.ucar.edu. Add your attached namelist file, and include difference plots of a couple of fields that you mentioned. I'll leave this issue open until we resolve what is happening.
The forum takes a day to manually process for first-time contributors. Afterwards the back-and-forth responses can be fairly quick.
I used WRF4.1.5 to simulate days with observed convection activity. One of the diagnostics I am interested in analysing is LPI (lightning_option = 3). I've noticed large differences in simulated fields (U, V, W, T, ...) with lightning_option = 3 on and off, even though LPI is a diagnostic variable. Has anyone seen this kind of problem before? If so, any recommendation on how to fix this problem?
Attached is my namelist.input
&time_control run_days = 1, run_hours = 12, run_minutes = 0, run_seconds = 0, start_year = 2017, 2017, 2017, 2017, start_month = 06, 06, 06, 06, start_day = 24, 24, 24, 24, start_hour = 12, 12, 12, 12, start_minute = 00, 00, 00, 00, start_second = 00, 00, 00, 00, end_year = 2017, 2017, 2017, 2017, end_month = 06, 06, 06, 06, end_day = 26, 26, 26, 26, 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, 10, 1, frames_per_outfile = 1000, 24, 6, 60, iofields_filename = "hail.txt","hail.txt","hail.txt","hail.txt", auxhist7_outname = "haild",
auxhist7_interval = 60, 60, 10, 1,
frames_per_auxhist7 = 1000, 1000, 6, 1,
io_form_auxhist7 = 2
restart = .false.,
restart_interval = 360,
write_hist_at_0h_rst = .true.
override_restart_timers = .true.
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
/
&domains time_step = 40, time_step_fract_num = 0, time_step_fract_den = 1, max_ts_locs = 40, max_dom = 3, e_we = 204, 304, 424, 265, 226! dimenzija u W-E smjeru e_sn = 180, 232, 289, 310, 226 e_vert = 66, 66, 66, 66, use_levels_below_ground = .true. use_surface = .true. smooth_cg_topo = .true. p_top_requested = 5000, num_metgrid_levels = 26, num_metgrid_soil_levels = 4, dx = 9000, 3000, 1000, 333.3333, dy = 9000, 3000, 1000, 333.3333, grid_id = 1, 2, 3, 4, parent_id = 0, 1, 2, 3, parent_grid_ratio = 1, 3, 3, 3, i_parent_start = 1, 48, 110, 64, ! ovo je u x-smjeru, slijeva-nadesno j_parent_start = 1, 57, 86, 48, ! ovo je u y-smjeru, odozdo-nagore parent_time_step_ratio = 1, 3, 3, 3, feedback = 0, smooth_option = 2, eta_levels = 1.00000, 0.99726, 0.99452, 0.99113, 0.98873, 0.98594, 0.98305, 0.98026, 0.97737, 0.97458, 0.97178, 0.96899, 0.96620, 0.96341, 0.95783, 0.95095, 0.93923, 0.92752, 0.91581, 0.90409, 0.87836, 0.85334, 0.82901, 0.80538, 0.77677, 0.74915, 0.72253, 0.69681, 0.67208, 0.64816, 0.62054, 0.59412, 0.56052, 0.52124, 0.48455, 0.45045, 0.41874, 0.38913, 0.36171, 0.31216, 0.29003, 0.26929, 0.25005, 0.23210, 0.21545, 0.19990, 0.18544, 0.16560, 0.14776, 0.13170, 0.11735, 0.10439, 0.09272, 0.07906, 0.06720, 0.05693, 0.04806, 0.04028, 0.03061, 0.02283, 0.01645, 0.01137, 0.00728, 0.00399, 0.00130, 0.000, /
&physics mp_physics = 10, 10, 10, 10, lightning_option = 0, 3, 3, 3, lightning_dt = 600, 600, 600, 600, lightning_start_seconds = 3600., 3600., 3600., 3600., flashrate_factor = 1, 1, 1, 1, cldtop_adjustment = 2., 2., 2., 2., iccg_method = 2, 2, 2, 2, do_radar_ref = 1, ra_lw_physics = 1, 1, 1, 1, ra_sw_physics = 1, 1, 1, 1, topo_shading = 1, radt = 9, 9, 9, 9, sf_sfclay_physics = 1, 1, 1, 1, sf_surface_physics = 4, 4, 4, 4, bl_pbl_physics = 8, 8, 8, 8, bldt = 0, 0, 0, 0, cu_physics = 1, 0, 0, 0, hailcast_opt = 0, 1, 1, 1, ! haildt = 0, 0, 0, 0, num_land_cat = 20, cudt = 5, 5, 5, 5, cugd_avedx = 3, isfflx = 1, ifsnow = 0, icloud = 1, surface_input_source = 1, num_soil_layers = 4, sf_urban_physics = 0, maxiens = 1, maxens = 3, maxens2 = 3, maxens3 = 16, ensdim = 144, /
&fdda /
&dynamics rk_ord = 3, w_damping = 1, diff_opt = 1, 1, 1, 1, km_opt = 4, 4, 4, 4, diff_6th_opt = 0, 0, 0, 0, diff_6th_factor = 0.12, 0.12, 0.12, 0.12, base_temp = 290. damp_opt = 3, zdamp = 5000., 5000., 5000., 5000., dampcoef = 0.2, 0.2, 0.2, 0.2, khdif = 0, 0, 0, 0, kvdif = 0, 0, 0, 0, epssm = 0.1, 0.1, 0.2, 0.2, 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, /