WRF ./real.exe报错 No average surface temperature for use with inland lakes

muyuqingye

刚接触WRF，想试试不同的参数化方案组合，但是./real.exe 报错了，有大神帮忙看看吗？

浏览WRF论坛有人说是静态地理数据缺失，我把官网的所有地理数据都解压缩到WPS_GEOG里重新运行还是这个错误；参考家园里
用ERA资料运行wrf出现错误，请各位大神帮帮忙，谢谢-专业气象研究-气象家园_气象人自己的家园 (06climate.com)
这篇中的指导添加surface_input_source参数也报错。

原本的错误：
rsl.error.0000文件如下：
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         153           1          61
ims,ime,jms,jme           -4         158          -4          66
ips,ipe,jps,jpe            1         153           1          61
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             432096776  bytes allocated
d01 2016-07-30_00:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.0
d01 2016-07-30_00:00:00  Input data is acceptable to use: met_em.d01.2016-07-30_00:00:00.nc
metgrid input_wrf.F first_date_input = 2016-07-30_00:00:00
metgrid input_wrf.F first_date_nml = 2016-07-30_00:00:00
d01 2016-07-30_00:00:00 Timing for input          0 s.
d01 2016-07-30_00:00:00          flag_soil_layers read from met_em file is  1
Max map factor in domain 1 =  1.03. Scale the dt in the model accordingly.

Using sfcprs  to compute psfc
d01 2016-07-30_00:00:00 No average surface temperature for use with inland lakes
Assume non-RUC LSM input
from Noah to RUC - compute Noah bucket
from Noah to RUC - compute RUC bucket
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    3062
grid%tsk unreasonable
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0



添加surface_input_source的报错：
module_io_quilt_old.F        2931 T
  ------ ERROR while reading namelist physics ------
Maybe here?:      sf_surface_physics                  = 3,    3,   3,
Maybe here?:      surface_input_source                = 1,    1,   1,
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:   10851
ERRORS while reading one or more namelists from namelist.input.
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0



namelist.input文件如下：
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2016, 2016,
start_month                         = 07,   07,
start_day                           = 30,   30,
start_hour                          = 00,   00,
end_year                            = 2016, 2016,
end_month                           = 07,   07,
end_day                             = 30,   30,
end_hour                            = 23,   23,
interval_seconds                    = 3600,
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 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                             = 2,
e_we                                = 153,   184,
e_sn                                = 61,    139,
e_vert                              = 50,    50,
p_top_requested                     = 5000,
num_metgrid_levels                  = 38,
num_metgrid_soil_levels             = 4,
dx                                  = 36000, 12000,
dy                                  = 36000, 12000,
grid_id                             = 1,     2,
parent_id                           = 1,     1,
i_parent_start                      = 1,     31,
j_parent_start                      = 1,     15,
parent_grid_ratio                   = 1,     3,
parent_time_step_ratio              = 1,     3,
feedback                            = 1,
smooth_option                       = 0,
/

&physics
mp_physics                          = 2,    2,   2,
cu_physics                          = 1,    1,   0,
ra_lw_physics                       = 1,    1,   1,
ra_sw_physics                       = 1,    1,   1,
bl_pbl_physics                      = 1,    1,   1,
sf_sfclay_physics                   = 1,    1,   1,
sf_surface_physics                  = 3,    3,   3,
surface_input_source              = 1,    1,   1, (这句是按照论坛里的帖子加的，但还是报错)
radt                                = 36,   36,   36,
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,
etac                                = 0.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.,.false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

ayzqs

Maybe here?:      sf_surface_physics                  = 3,    3,   3,
Maybe here?:      surface_input_source                = 1,    1,   1,

是不是逗号的问题？

muyuqingye

ayzqs 发表于 2023-3-16 17:06
Maybe here?:      sf_surface_physics                  = 3,    3,   3,
Maybe here?:      surface_inp ...

更换Noah方案和热交换方案都可以运行，应该不是逗号的问题，应该是RUC方案需要的陆面数据比较多吧，我正在尝试。

muyuqingye

ayzqs 发表于 2023-3-16 17:06
Maybe here?:      sf_surface_physics                  = 3,    3,   3,
Maybe here?:      surface_inp ...

您好，我改成一层嵌套，试了Noah方案和热传输的陆面方案都可以成功，但RUC方案还是不行，这是与RUC方案需要的路面数据多有关吗？

铃铛tsuki

我感觉改成surface_input_source   = 1, （不要写= 1,    1,   1,）