ali.waseem607
New member
Hello,
I am simulating a mountainous region where at certain times it is expected that the top of the mount enters the free troposphere and is above the PBLH. I was plotting the values from wrf and the vlaues look odd with discontinuities present. I am using the YSU scheme for the boundary layer physics. A graph of PBLH as well as settings in the namelist input file are provided below. I wanted to know if
1. it is possible to have negative values for PBLH?
2. if it isn't possible is there an alternative way to determine if a region is above the boundary layer?
3. If the values obtained from my WRF simulation are expected or if something looks incorrect. The two lines are from different locations in the wrf file.
Thank you for any help that you can provide!
Namelist inputs:
&time_control
run_days = 0,
run_hours = 66,
run_minutes = 0,
run_seconds = 0,
start_year = 2024, 2024, 2024,
start_month = 10, 10, 10,
start_day = 08, 08, 08,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2024, 2024, 2024,
end_month = 10, 10, 10,
end_day = 10, 10, 10,
end_hour = 18, 18, 18,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 30, 30, 5,
frames_per_outfile = 1000, 1000, 1000,
history_outname = './wrfout_d<domain>_<date>.nc',
restart = .false.,
restart_interval = 720,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 1000,
/
&domains
time_step = 36,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 148, 241, 304,
e_sn = 148, 241, 304,
e_vert = 97, 97, 97,
!97 levels
eta_levels = 1, 0.9938147, 0.9859506, 0.9760143, 0.9635575, 0.9483238, 0.9333184, 0.9185379, 0.903979, 0.8896384,
0.8755127, 0.8615989, 0.8478936, 0.8343937, 0.8210963, 0.8079981, 0.7950964, 0.782388, 0.7698702,
0.75754, 0.7453946, 0.7334313, 0.7216474, 0.7100402, 0.6986068, 0.687345, 0.6762519, 0.6653252, 0.6545622,
0.6439605, 0.6335179, 0.6232318, 0.6130998, 0.6031198, 0.5932894, 0.5836063, 0.5740684, 0.5646735, 0.5554194,
0.546304, 0.5373253, 0.5284812, 0.5197697, 0.5111888, 0.5027365, 0.4944109, 0.4862102, 0.4781324, 0.4701756,
0.4623381, 0.4546182, 0.447014, 0.4395238, 0.4321458, 0.4248784, 0.4177201, 0.410669, 0.4037236, 0.3968824,
0.3901436, 0.383506, 0.3769678, 0.3705276, 0.364184, 0.3579355, 0.3517806, 0.3457181, 0.3397464, 0.3338642,
0.3280702, 0.3223631, 0.3167415, 0.3112042, 0.3055, 0.299, 0.2915, 0.28120, 0.2668, 0.249, 0.230, 0.21, 0.1880719,
0.1648721, 0.1439084, 0.1249652, 0.1078478, 0.09238029, 0.07840355, 0.06577397, 0.05436168, 0.0440493, 0.0347309,
0.02631061, 0.01870192, 0.01182657, 0.005613877, 0,
p_top_requested = 5000,
smooth_cg_topo = .true.,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 12000., 3000., 1000.,
dy = 12000., 3000., 1000.,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 45, 68,
j_parent_start = 1, 45, 64,
parent_grid_ratio = 1, 4, 3,
parent_time_step_ratio = 1, 4, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 10, 10, 10,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 9, 3, 1,
sf_sfclay_physics = 1, 1, 1, !2,2,2 5,5,5
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1, !YSU=1,1,1 !MYJ:2,2,2, 5,5,5
bldt = 0, 0, 0,
cu_physics = 1, 0, 0,
cudt = 5, 0, 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
fractional_seaice = 0,
SST_UPDATE = 0,
SEAICE_ALBEDO_OPT = 0,
SEAICE_ALBEDO_DEFAULT = 0.65
/
&fdda
grid_fdda = 1,1,1,
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360,360,360,
gfdda_end_h = 6,6,6,
io_form_gfdda = 2,
fgdt = 0,0,0,
if_no_pbl_nudging_uv = 1,1,1,
if_no_pbl_nudging_t = 1,1,1,
if_no_pbl_nudging_q = 1,1,1,
if_zfac_uv = 0,1,1,
k_zfac_uv = 27,
if_zfac_t = 1,1,1,
k_zfac_t = 10,
if_zfac_q = 1,
k_zfac_q = 10,
guv = 0.0000309,0.0000231,0.0000231,
gt = 0.0000309,0.0,0.0,
gq = 0.0000309,0.0,0.0,
if_ramping = 0,
dtramp_min = 60.0,
/
&dynamics
w_damping = 1,
diff_opt = 2, 2, 2,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
epssm = 0.3, 0.3, 0.3, !0.5
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
tke_adv_opt = 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&diags
!p_lev_diags = 1
!num_press_levels = 4,
!press_levels = 85000, 50000, 20000
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am simulating a mountainous region where at certain times it is expected that the top of the mount enters the free troposphere and is above the PBLH. I was plotting the values from wrf and the vlaues look odd with discontinuities present. I am using the YSU scheme for the boundary layer physics. A graph of PBLH as well as settings in the namelist input file are provided below. I wanted to know if
1. it is possible to have negative values for PBLH?
2. if it isn't possible is there an alternative way to determine if a region is above the boundary layer?
3. If the values obtained from my WRF simulation are expected or if something looks incorrect. The two lines are from different locations in the wrf file.
Thank you for any help that you can provide!
Namelist inputs:
&time_control
run_days = 0,
run_hours = 66,
run_minutes = 0,
run_seconds = 0,
start_year = 2024, 2024, 2024,
start_month = 10, 10, 10,
start_day = 08, 08, 08,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2024, 2024, 2024,
end_month = 10, 10, 10,
end_day = 10, 10, 10,
end_hour = 18, 18, 18,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.true.,.true.,
history_interval = 30, 30, 5,
frames_per_outfile = 1000, 1000, 1000,
history_outname = './wrfout_d<domain>_<date>.nc',
restart = .false.,
restart_interval = 720,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 1000,
/
&domains
time_step = 36,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 3,
e_we = 148, 241, 304,
e_sn = 148, 241, 304,
e_vert = 97, 97, 97,
!97 levels
eta_levels = 1, 0.9938147, 0.9859506, 0.9760143, 0.9635575, 0.9483238, 0.9333184, 0.9185379, 0.903979, 0.8896384,
0.8755127, 0.8615989, 0.8478936, 0.8343937, 0.8210963, 0.8079981, 0.7950964, 0.782388, 0.7698702,
0.75754, 0.7453946, 0.7334313, 0.7216474, 0.7100402, 0.6986068, 0.687345, 0.6762519, 0.6653252, 0.6545622,
0.6439605, 0.6335179, 0.6232318, 0.6130998, 0.6031198, 0.5932894, 0.5836063, 0.5740684, 0.5646735, 0.5554194,
0.546304, 0.5373253, 0.5284812, 0.5197697, 0.5111888, 0.5027365, 0.4944109, 0.4862102, 0.4781324, 0.4701756,
0.4623381, 0.4546182, 0.447014, 0.4395238, 0.4321458, 0.4248784, 0.4177201, 0.410669, 0.4037236, 0.3968824,
0.3901436, 0.383506, 0.3769678, 0.3705276, 0.364184, 0.3579355, 0.3517806, 0.3457181, 0.3397464, 0.3338642,
0.3280702, 0.3223631, 0.3167415, 0.3112042, 0.3055, 0.299, 0.2915, 0.28120, 0.2668, 0.249, 0.230, 0.21, 0.1880719,
0.1648721, 0.1439084, 0.1249652, 0.1078478, 0.09238029, 0.07840355, 0.06577397, 0.05436168, 0.0440493, 0.0347309,
0.02631061, 0.01870192, 0.01182657, 0.005613877, 0,
p_top_requested = 5000,
smooth_cg_topo = .true.,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 12000., 3000., 1000.,
dy = 12000., 3000., 1000.,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
i_parent_start = 1, 45, 68,
j_parent_start = 1, 45, 64,
parent_grid_ratio = 1, 4, 3,
parent_time_step_ratio = 1, 4, 3,
feedback = 1,
smooth_option = 0
/
&physics
mp_physics = 10, 10, 10,
ra_lw_physics = 4, 4, 4,
ra_sw_physics = 4, 4, 4,
radt = 9, 3, 1,
sf_sfclay_physics = 1, 1, 1, !2,2,2 5,5,5
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1, !YSU=1,1,1 !MYJ:2,2,2, 5,5,5
bldt = 0, 0, 0,
cu_physics = 1, 0, 0,
cudt = 5, 0, 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
fractional_seaice = 0,
SST_UPDATE = 0,
SEAICE_ALBEDO_OPT = 0,
SEAICE_ALBEDO_DEFAULT = 0.65
/
&fdda
grid_fdda = 1,1,1,
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360,360,360,
gfdda_end_h = 6,6,6,
io_form_gfdda = 2,
fgdt = 0,0,0,
if_no_pbl_nudging_uv = 1,1,1,
if_no_pbl_nudging_t = 1,1,1,
if_no_pbl_nudging_q = 1,1,1,
if_zfac_uv = 0,1,1,
k_zfac_uv = 27,
if_zfac_t = 1,1,1,
k_zfac_t = 10,
if_zfac_q = 1,
k_zfac_q = 10,
guv = 0.0000309,0.0000231,0.0000231,
gt = 0.0000309,0.0,0.0,
gq = 0.0000309,0.0,0.0,
if_ramping = 0,
dtramp_min = 60.0,
/
&dynamics
w_damping = 1,
diff_opt = 2, 2, 2,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2,
epssm = 0.3, 0.3, 0.3, !0.5
khdif = 0, 0, 0,
kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
moist_adv_opt = 1, 1, 1,
scalar_adv_opt = 1, 1, 1,
tke_adv_opt = 1, 1, 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,
nested = .false., .true., .true.,
/
&grib2
/
&diags
!p_lev_diags = 1
!num_press_levels = 4,
!press_levels = 85000, 50000, 20000
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/