Hello
I am running wrf 4.2 for a coastal domain (city and port facilities)
I get fairly low mixing heights (20-40 m above ground level) overnight at port site stations, even though I would expect higher PBLs because of the marine boundary layer influence. I do not understand what could be going on wrong.
The simulation runs on may (austral fall season)
Below I post the namefile physics section
Regards
Hector
&physics
sst_update = 1
rdlai2d = .true.
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 12, 12, 12,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
ysu_topdown_pblmix = 1,
topo_wind = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 2, 0, 0,
cu_diag = 0, 0, 0,
cudt = 0, 0, 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 20,
sf_urban_physics = 0, 0, 0,
slope_rad = 1, 1, 1,
topo_shading = 1, 1, 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false., .false., .false.,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-8,
/
I am running wrf 4.2 for a coastal domain (city and port facilities)
I get fairly low mixing heights (20-40 m above ground level) overnight at port site stations, even though I would expect higher PBLs because of the marine boundary layer influence. I do not understand what could be going on wrong.
The simulation runs on may (austral fall season)
Below I post the namefile physics section
Regards
Hector
&physics
sst_update = 1
rdlai2d = .true.
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 12, 12, 12,
sf_sfclay_physics = 1, 1, 1,
sf_surface_physics = 2, 2, 2,
bl_pbl_physics = 1, 1, 1,
ysu_topdown_pblmix = 1,
topo_wind = 1, 1, 1,
bldt = 0, 0, 0,
cu_physics = 2, 0, 0,
cu_diag = 0, 0, 0,
cudt = 0, 0, 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 20,
sf_urban_physics = 0, 0, 0,
slope_rad = 1, 1, 1,
topo_shading = 1, 1, 1,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false., .false., .false.,
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-8,
/