GFS资料转WRF ./real.exe显示mismatch_landmask_ivgtyp错误

jljyyh

尝试用GFS（格式：gfs.0p25.YYYYMMDDHH.f000.grib2）资料转WRF，前面的步骤都没有出现问题，但是在运行./real.exe查看rsl.erroe.00文件发现显示mismatch_landmask_ivgtyp问题，使用的WRF版本是V3.9.1
查看网上的解决办法，是将namelist.input中surface_input_source 改成1，但我的那么namelist.input默认就为1，并且后来查到wrf解决问题的官网上说的是初始场用欧洲中心的数据有出现过这个问题
http://bbs.06climate.com/forum.php?mod=viewthread&tid=46183

有考虑是不是GFS资料的问题，但查找了WRF上GFS资料的使用（http://www2.mmm.ucar.edu/wrf/OnLineTutorial/DATA/GFS/index.php），没有发现什么异常（∵之前有考虑是不是Vtable.GFS变了），然后更换官网上的第一个例子（Jan 2000）的例子，同样出现了这个问题（由此排除了初始场、资料、分辨率、积分步长的问题），所以想麻烦大家帮我看一下问题出现在哪里？感谢大家~~~

附上最后调试的namelist.wps和namelist.input
namelist.wps
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2000-01-24_12:00:00','2017-08-24_00:00:00','2017-08-25_00:00:00'
end_date   = '2000-01-25_12:00:00','2017-08-25_00:00:00','2017-08-25_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, 138,  95
j_parent_start    =   1, 128, 107
e_we              = 74, 334, 421
e_sn              = 61, 343, 397
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the HGT_M, GREENFRAC,
! and LU_INDEX/LANDUSEF fields have changed in WPS v3.8. The HGT_M field
! is now interpolated from 30-arc-second USGS GMTED2010, the GREENFRAC
! field is interpolated from MODIS FPAR, and the LU_INDEX/LANDUSEF fields
! are interpolated from 21-class MODIS.
!
! To match the output given by the default namelist.wps in WPS v3.7.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'gtopo_10m+usgs_10m+nesdis_greenfrac+10m','gtopo_2m+usgs_2m+nesdis_greenfrac+2m',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','5m','2m'
dx = 30000,
dy = 30000,
map_proj = 'lambert',
ref_lat   =  34.83,
ref_lon   = -81.03,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon = -98.0,
geog_data_path = '/home/WRF/geog/'
/

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

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,


namelist.input
&time_control
run_days                            = 0,
run_hours                           = 24,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2000, 2017, 2017,
start_month                         = 01,   08,   08,
start_day                           = 24,   24,   24,
start_hour                          = 12,   00,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2000, 2017, 2017,
end_month                           = 01,   08,   08,
end_day                             = 25,   25,   25,
end_hour                            = 12,   00,   00,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,   60,   60,
frames_per_outfile                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 74,   334,   421,
e_sn                                = 61,   343,   397,
e_vert                              = 33,    30,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 2,
dx                                  = 30000, 3000, 1000,
dy                                  = 30000, 3000, 1000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,   138,    95,
j_parent_start                      = 1,   128,   107,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 4,      4,     4,
ra_lw_physics                       = 1,      1,     1,
ra_sw_physics                       = 1,      1,     1,
&physics_suite                       = 'CONUS'
radt                                = 30,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 1,     0,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
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                            = 0,
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,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

werewolf

阿羊出没请注意 发表于 2021-11-26 09:28
您好，我用的就是WRF4版本，仍然出现问题请问该如何解决呀？网上搜的办法，我试了依然报错。。

用的GSI吗，出现了什么问题

werewolf

用WRF4版本

jljyyh

werewolf 发表于 2019-7-30 08:42
用WRF4版本

好的，我尝试用WRFV4再试一次

jljyyh

werewolf 发表于 2019-7-30 08:42
用WRF4版本

十分感谢~~

冽xin

请问你0.25的资料在哪下的呢

阿羊出没请注意

werewolf 发表于 2019-7-30 08:42
用WRF4版本

您好，我用的就是WRF4版本，仍然出现问题请问该如何解决呀？网上搜的办法，我试了依然报错。。

阿羊出没请注意

werewolf 发表于 2021-11-26 15:11
用的GSI吗，出现了什么问题

我用的ERA-5的逐小时pressure和single level的数据，在./real.exe时，出现报错信息：


the grid%landmask and soil/veg cats do not match
i,j=         106           8
grid%landmask=   0.00000000
grid%ivgtyp=          17
grid%isltyp=          12
iswater=          17
grid%tslb=   0.00000000       0.00000000       0.00000000       0.00000000       0.00000000
grid%sst=   0.00000000
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    3388
mismatch_landmask_ivgtyp

是不是WRF4.0的存在bug？我数据里不包含SST时就可以运行，加了SST信息就出现这种报错

myj0219

感谢分享！！！！

sillynana11

阿羊出没请注意 发表于 2021-11-26 18:34
我用的ERA-5的逐小时pressure和single level的数据，在./real.exe时，出现报错信息：

请问解决问题了吗 如何解决的呢