./real.exe报错，num_metgrid_levels不符，请求赐教

亚当i

之前用FNL一直跑的可以，刚改用了EAR5的欧洲气象数据，WPS部分都运行的很快，需要链接的VTABLE和数据都更换了：ln -sf ungrib/Variable_Tables/Vtable.ECMWF Vtablemet_em文件复制到WRF开始报错。显示num_metgrid_levels不同，不知道为什么WPS处理完的层数只有一层，请教大家大家如何改正？



d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  namelist   ide,jde,num_metgrid_levels=          72          73          34
d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  input file ide,jde,kde               =          72          73           1
d01 2022-06-01_00:00:00 ---- ERROR: Mismatch between namelist and input file dimensions
NOTE:       1 namelist vs input data inconsistencies found.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1248
NOTE:  Please check and reset these options


附上namelist.wps
&share wrf_core = 'ARW',
max_dom = 2,
start_date = '2022-06-01_00:00:00','2022-06-01_00:00:00','2022-05-30_00:00:00',
end_date   = '2022-06-20_00:00:00','2022-06-20_00:00:00','2022-06-10_00:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/

&geogrid
parent_id         =   1,   1,  2
parent_grid_ratio =   1,   3,  3
i_parent_start    =   1,  20,  36
j_parent_start    =   1,  21,  35
e_we              =  72,  94, 49
e_sn              =  73,  61, 73
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the MAXSNOALB and ALBEDO12M
! fields have changed in WPS v4.0. These fields are now interpolated
! from MODIS-based datasets.
!
! To match the output given by the default namelist.wps in WPS v3.9.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','default',
dx = 27000,
dy = 27000,
map_proj = 'lambert',
ref_lat   =  32.627,
ref_lon   = 112.601,
truelat1  =  30.0,
truelat2  =  35.5,
stand_lon = 112.601,
geog_data_path = '/home/test/models/Build_WRF/WPS_GEO'
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/



---------
附上namelist.input：
&time_control
run_days                            = 20,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2022, 2022, 2022,
start_month                         = 06,   06,   05,
start_day                           = 01,   01,   30,
start_hour                          = 00,   00,   00,
end_year                            = 2022, 2022, 2022,
end_month                           = 06,   06,   06,
end_day                             = 20,   20,   10,
end_hour                            = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 60,  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                           = 150,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 72,   94,    49,
e_sn                                = 73,   61,    73,
e_vert                              = 33,    33,    33,
p_top_requested                     = 5000,
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
dx                                  = 27000, 9000,  3000,
dy                                  = 27000, 9000,  3000,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     20,    36,
j_parent_start                      = 1,     21,    35,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
physics_suite                       = 'CONUS'
mp_physics                          =  6,     6,     3,
cu_physics                          =  1,     1,     0,
ra_lw_physics                       =  1,     1,     1,
ra_sw_physics                       =  1,     1,     1,
bl_pbl_physics                      =  1,     1,     2,
sf_sfclay_physics                   =  1,     1,     2,
sf_surface_physics                  =  2,     2,     2,
radt                                = 27,    27,    27,
bldt                                = 0,     0,     0,
cudt                                = 0,     0,     0,
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,
/

--------

白啾啾

试试查看生成的met_em文件？这个命令ncdump -h met_em.d01.2023-04-02_00\:00\:00.nc >& log.ncdump
我在log.ncdump里看到了met_em文件里num_metgrid_levels 的值，好像得设置得和met_em文件里的一样

menseye

我记得欧洲的ERA5资料好像有38层，num_metgrid_levels                  = 34,这儿改一下
还有vtable好像错了，好像是这个：Vtable.ERA-interim.pl  Vtable

亚当i

menseye 发表于 2023-4-4 23:00
我记得欧洲的ERA5资料好像有38层，num_metgrid_levels                  = 34,这儿改一下
还有vtable好像 ...

d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  namelist   ide,jde,num_metgrid_levels=          72          73          38
d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  input file ide,jde,kde               =          72          73           1
d01 2022-06-01_00:00:00 ---- ERROR: Mismatch between namelist and input file dimensions
NOTE:       1 namelist vs input data inconsistencies found.

可是还是报错了，这个wps为啥生成的是1层。

menseye

亚当i 发表于 2023-4-5 13:27
d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  namelist   ide,jde,num_metgrid_levels=       ...

查一查下载的ERA5资料有几层，另外查一下met文件

亚当i

亚当i 发表于 2023-4-5 13:27
d01 2022-06-01_00:00:00  input_wrf.F: SIZE MISMATCH:  namelist   ide,jde,num_metgrid_levels=       ...

又报了个别的错
ERROR in psfc: flag_psfc =  1, flag_soilhgt =  0, flag_slp =  0, sfcp_to_sfcp = T
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1396
not enough info for a p sfc computation
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1
:
system msg for write_line failure : Bad file descriptor

亚当i

menseye 发表于 2023-4-5 17:41
查一查下载的ERA5资料有几层，另外查一下met文件

ERROR in psfc: flag_psfc =  1, flag_soilhgt =  0, flag_slp =  0, sfcp_to_sfcp = T
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1396
not enough info for a p sfc computation
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1
:
system msg for write_line failure : Bad file descriptor
又报了个这个错，麻烦给看一下

menseye

这个错我也不太懂