Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

No average surface temperature for use with inland lakes

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

atique9076

New member
I have run a model on FNL data of 5 days.
I want to apply different PBL schemes i.e YSU,MYJ, Boulac etc.
I am getting error when I selected bl_pbl_scheme=1 and all other scheme =0.
Using sfcprs3 to compute psfc
d01 2017-01-01_00:00:00 No average surface temperature for use with inland lakes
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 3415
mismatch_landmask_ivgtyp

Please look over it and help me.
Thanks

Attached- Namelist.Input
Physics scheme screenshot
Note:- If you can guide me how to apply physics(pbl) scheme in model, I will be grateful because I am new to applying physics scheme.
I read user guide also in this regard but that give me only information not how to do.
 

Attachments

  • namelist.input.png
    namelist.input.png
    154.3 KB · Views: 1,296
  • physics_schemes.png
    physics_schemes.png
    194.5 KB · Views: 1,294
&time_control

run_days = 04,

run_hours = 18,
run_minutes = 0,
run_seconds = 0,

start_year = 2017, 2000, 2000,
start_month = 01, 01, 01,
start_day = 01, 24, 24,
start_hour = 00, 12, 12,
end_year = 2017, 2000, 2000,
end_month = 01, 01, 01,
end_day = 05, 25, 25,
end_hour = 18, 12, 12,

interval_seconds = 10800,
input_from_file = .true.,.true.,.true.,
history_interval = 180, 60, 60,
frames_per_outfile = 1000, 1000, 1000,
restart = .false.,
restart_interval = 7200,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
/

&domains

time_step = 180,

time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,

e_we = 33, 112, 94,

e_sn = 23, 97, 91,
e_vert = 33, 33, 33,
p_top_requested = 5000,

num_metgrid_levels = 32,


num_metgrid_soil_levels = 4,
dx = 15000, 10000, 3333.33,
dy = 15000, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 31, 30,
j_parent_start = 1, 17, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/

&physics
mp_physics = 0, -1, -1,
cu_physics = 0, -1, 0,
ra_lw_physics = 0, -1, -1,
ra_sw_physics = 0, -1, -1,
bl_pbl_physics = 1, -1, -1,
sf_sfclay_physics = 0, -1, -1,
sf_surface_physics = 0, -1, -1,
radt = 30, 30, 30,
bldt = 0, 0, 0,
cudt = 5, 5, 5,
icloud = 1,
num_land_cat = 21,
sf_urban_physics = 0, 0, 0,
/

&fdda
/

&dynamics
hybrid_opt = 2,
w_damping = 0,
diff_opt = 1, 1, 1,
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
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,
gwd_opt = 1,
/

&bdy_control
spec_bdy_width = 5,
specified = .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
 
This was error showing in rsl.error.0000

askid: 0 hostname: Bhishma-Precision-Tower-7910
module_io_quilt_old.F 2931 T
Ntasks in X 1 , ntasks in Y 1
REAL_EM V4.0 PREPROCESSOR
*************************************
Parent domain
ids,ide,jds,jde 1 33 1 23
ims,ime,jms,jme -4 38 -4 28
ips,ipe,jps,jpe 1 33 1 23
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
alloc_space_field: domain 1 , 37820044 bytes allocated
d01 2017-01-01_00:00:00 Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.0
d01 2017-01-01_00:00:00 Input data is acceptable to use: met_em.d01.2017-01-01_00:00:00.nc
metgrid input_wrf.F first_date_input = 2017-01-01_00:00:00
metgrid input_wrf.F first_date_nml = 2017-01-01_00:00:00
d01 2017-01-01_00:00:00 Timing for input 0 s.
d01 2017-01-01_00:00:00 flag_soil_layers read from met_em file is 1
d01 2017-01-01_00:00:00 Turning off use of MAX WIND level data in vertical interpolation
d01 2017-01-01_00:00:00 Turning off use of TROPOPAUSE level data in vertical interpolation
Max map factor in domain 1 = 0.94. Scale the dt in the model accordingly.
Using sfcprs3 to compute psfc
d01 2017-01-01_00:00:00 No average surface temperature for use with inland lakes
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 3415
mismatch_landmask_ivgtyp
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
 
I wonder why you only run with PBl scheme? In this case the physics is not completed and real data case will not run as expected.
Can you turn on all necessary physics and retry the case?
 
Top