No average surface temperature for use with inland lakes

[atique9076]

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.

 

[mcurry]

This post was posted under MPAS Science. I have moved it to the WRF section of the Forum.

 

[Ming Chen]

Please upload your namelist.input. The one you attached is not complete and hard to read.

 

[atique9076]

&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,
/

 

[atique9076]

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

 

[Ming Chen]

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?