[Wrf-users] Unrealistic PBL heights when using CLM as compared to NOAH
Jatin Kala
jatin.kala.jk at gmail.com
Wed Nov 29 02:24:54 MST 2017
Hi all,
We ran WRF simulations (3.8.1) which are identical, except for the LSM
choice, i.e., NOAH versus CLM.
While the simulation using NOAH look very reasonable, the one using CLM
looks rather strange.
CLM_PBLH_ncview.png shows PBL heights of about 900 m within the
relaxation zone, but about 100-200 m inside.
While i always expect some discontinuity within the relaxation zone,
this looks spurious.
NOAH_PBLH_ncview.png shows PBL heights at the same output time, and we
do not have this issue.
The same thing issue reflected in the sensible heat flux (HFX) outputs
(also attached).
CLM_HFX_ncview.png shows higher sensible heat flux over inland water
points as compared to the land points.....
NOAH_HFX_ncview.png shows near zero sensible heat flux over the inland
water points, and much higher over land points, which makes a lot more
sense.
When i look at SWDOWN from the CLM outputs (CLM_SWDOWN_ncview.png
attached), it's exactly zero across the domain, except over the inland
water point where it is 600 Wm-2. How is this possible??
SWDOWN for NOAH looks fine (NOAH_SWDOWN_ncview.png attached).
Has anyone else experienced similar issues when using CLM in WRF? The
namelists are attached, and the only difference is sf_surface_physics
and num_soil_layers:
134:temp_CLM Jatin$ diff namelist.input_CLM
../temp_NOAH/namelist.input.NOAH
71c71
< sf_surface_physics = 5, 5, 5, 5,
---
> sf_surface_physics = 2, 2, 2, 2,
80c80
< num_soil_layers = 10,
---
> num_soil_layers = 4,
Both simulations are driven with Era-Interim. real.exe was run for each
simulation separately, using the same met_em files.
tar files with all input and few output files can be downloaded from
here if that's useful:
https://www.dropbox.com/s/nigy40p1pqb1sdz/CLM.tar?dl=0
https://www.dropbox.com/s/b1dfs3wum5ytb8p/NOAH.tar?dl=0
Thanks in advance,
Jatin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CLM_HFX_ncview.png
Type: image/png
Size: 217740 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0006.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CLM_PBLH_ncview.png
Type: image/png
Size: 207141 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0007.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NOAH_HFX_ncview.png
Type: image/png
Size: 150099 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0008.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NOAH_PBLH_ncview.png
Type: image/png
Size: 55492 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0009.png>
-------------- next part --------------
&time_control
run_days = 46,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2005, 2005, 2005 2005,
start_month = 11, 11, 11, 11
start_day = 15, 15, 15, 15
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2005, 2005, 2005, 2005
end_month = 12, 12, 12, 12
end_day = 31, 31, 31, 31
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 = 24, 24, 24, 24,
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
io_form_auxinput4 = 2 ! this is to use sst data
auxinput4_inname = "wrflowinp_d<domain>" ! this is to use sst data
auxinput4_interval = 360,360,360, 360, ! to use sst data
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 4,
e_we = 150,142,202,187
e_sn = 105,112,133,106
e_vert = 45, 45, 45, 45,
eta_levels = 1.00,0.995,0.99,0.985,0.98,0.975,0.97,0.965,0.96,0.955,0.95,0.945,0.94,0.935,0.93,0.925,0.92,0.91,0.90,0.89,0.88,0.87,0.86,0.85,0.84,0.83,0.80,0.77,0.72,0.67,0.62,0.57,0.52,0.47,0.42,0.37,0.32,0.27,0.22,0.17,0.12,0.07,0.04,0.02,0.00
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000, 1000,
dy = 27000, 9000, 3000, 1000,
grid_id = 1, 2, 3, 4,
parent_id = 0, 1, 2, 3,
i_parent_start = 1, 51,30,120
j_parent_start = 1, 31,37,74
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 0,
smooth_option = 0,
/
&physics
mp_physics = 4, 4, 4, 4,
mp_zero_out = 2 ! from J.Evans namelist
mp_zero_out_thresh = 1.e-8 ! from J.Evans namelist
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 10, 10, 10,
cam_abs_freq_s = 10800
levsiz = 59 ! same as above comment
paerlev = 29 ! same as above comment
cam_abs_dim1 = 4 ! same as above comment
cam_abs_dim2 = 28 ! same as above comment
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 = 2, 2, 0, 0,
cudt = 0, 0, 0, 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
sst_update = 1, ! from J.Evans script to use sst data
sst_skin = 1,
tmn_update = 1,
lagday = 150,
usemonalb = .true. ! from J.Evans script, this is to use monthly albedo maps
slope_rad = 1, ! from J.Evans script, this turns on slope effects for ra_sw_physics
sf_urban_physics = 0, 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
grid_fdda = 0, 0, 0, 0, ! apply spectral nudging to outer grid
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360,
fgdt = 0,
xwavenum = 5, ! total grid size in x divide by 1000
ywavenum = 4, ! total grid size in y divide by 1000
if_no_pbl_nudging_ph = 1, ! no nudge in pbl
if_no_pbl_nudging_uv = 1,
if_no_pbl_nudging_t = 1,
if_zfac_ph = 1, ! nudge above k_zfac_ph only
k_zfac_ph = 10,
dk_zfac_ph = 1,
if_zfac_uv = 1,
k_zfac_uv = 10,
dk_zfac_uv = 1,
if_zfac_t = 1,
k_zfac_t = 10,
dk_zfac_t = 1,
gph = 0.0003,
guv = 0.0003,
gt = 0.0003,
io_form_gfdda = 2,
/
&dynamics
rk_ord = 3, ! from J.Evans script
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 = 3,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.05, 0.05, 0.05, 0.05,
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 = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false.,.false.,.false.,
nested = .false., .true., .true.,.true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
-------------- next part --------------
&time_control
run_days = 46,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2005, 2005, 2005 2005,
start_month = 11, 11, 11, 11
start_day = 15, 15, 15, 15
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2005, 2005, 2005, 2005
end_month = 12, 12, 12, 12
end_day = 31, 31, 31, 31
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 = 24, 24, 24, 24,
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
io_form_auxinput4 = 2 ! this is to use sst data
auxinput4_inname = "wrflowinp_d<domain>" ! this is to use sst data
auxinput4_interval = 360,360,360, 360, ! to use sst data
/
&domains
time_step = 60,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 4,
e_we = 150,142,202,187
e_sn = 105,112,133,106
e_vert = 45, 45, 45, 45,
eta_levels = 1.00,0.995,0.99,0.985,0.98,0.975,0.97,0.965,0.96,0.955,0.95,0.945,0.94,0.935,0.93,0.925,0.92,0.91,0.90,0.89,0.88,0.87,0.86,0.85,0.84,0.83,0.80,0.77,0.72,0.67,0.62,0.57,0.52,0.47,0.42,0.37,0.32,0.27,0.22,0.17,0.12,0.07,0.04,0.02,0.00
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 27000, 9000, 3000, 1000,
dy = 27000, 9000, 3000, 1000,
grid_id = 1, 2, 3, 4,
parent_id = 0, 1, 2, 3,
i_parent_start = 1, 51,30,120
j_parent_start = 1, 31,37,74
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 0,
smooth_option = 0,
/
&physics
mp_physics = 4, 4, 4, 4,
mp_zero_out = 2 ! from J.Evans namelist
mp_zero_out_thresh = 1.e-8 ! from J.Evans namelist
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 10, 10, 10,
cam_abs_freq_s = 10800
levsiz = 59 ! same as above comment
paerlev = 29 ! same as above comment
cam_abs_dim1 = 4 ! same as above comment
cam_abs_dim2 = 28 ! same as above comment
sf_sfclay_physics = 1, 1, 1, 1,
sf_surface_physics = 5, 5, 5, 5,
bl_pbl_physics = 1, 1, 1, 1,
bldt = 0, 0, 0, 0,
cu_physics = 2, 2, 0, 0,
cudt = 0, 0, 0, 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 10,
num_land_cat = 24,
sst_update = 1, ! from J.Evans script to use sst data
sst_skin = 1,
tmn_update = 1,
lagday = 150,
usemonalb = .true. ! from J.Evans script, this is to use monthly albedo maps
slope_rad = 1, ! from J.Evans script, this turns on slope effects for ra_sw_physics
sf_urban_physics = 0, 0, 0, 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
grid_fdda = 0, 0, 0, 0, ! apply spectral nudging to outer grid
gfdda_inname = "wrffdda_d<domain>",
gfdda_interval_m = 360,
fgdt = 0,
xwavenum = 5, ! total grid size in x divide by 1000
ywavenum = 4, ! total grid size in y divide by 1000
if_no_pbl_nudging_ph = 1, ! no nudge in pbl
if_no_pbl_nudging_uv = 1,
if_no_pbl_nudging_t = 1,
if_zfac_ph = 1, ! nudge above k_zfac_ph only
k_zfac_ph = 10,
dk_zfac_ph = 1,
if_zfac_uv = 1,
k_zfac_uv = 10,
dk_zfac_uv = 1,
if_zfac_t = 1,
k_zfac_t = 10,
dk_zfac_t = 1,
gph = 0.0003,
guv = 0.0003,
gt = 0.0003,
io_form_gfdda = 2,
/
&dynamics
rk_ord = 3, ! from J.Evans script
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 = 3,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.05, 0.05, 0.05, 0.05,
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 = 10,
spec_zone = 1,
relax_zone = 9,
specified = .true., .false.,.false.,.false.,
nested = .false., .true., .true.,.true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: CLM_SWDOWN_ncview.png
Type: image/png
Size: 143288 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0010.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NOAH_SWDOWN_ncview.png
Type: image/png
Size: 167064 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/wrf-users/attachments/20171129/edfedd28/attachment-0011.png>
More information about the Wrf-users
mailing list