WPS在LANDSEA上的异常

yangshujian

在跑WPS之后出来的met_em文件里有LANDSEA这个变量，海洋部分与陆地部分有明显的垂直分界线，在海洋部分是0，陆地部分是1.这样的后果就是后续做WRF之后的结果也会有明显的分界。具体请看图。一并贴上我的namelist.wps和namelist.input。各位同仁有遇到这种情况的吗，请问有可能是哪里出了问题。


namelist.input

&time_control
run_days                            = 1,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2018, 2000, 2000,
start_month                         = 11,   01,   01,
start_day                           = 29,   24,   24,
start_hour                          = 00,   12,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2018, 2000, 2000,
end_month                           = 11,   01,   01,
end_day                             = 30,   25,   25,
end_hour                            = 00,   12,   12,
end_minute                          = 00,   00,   00,
end_second                          = 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                    = 5000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 120,    112,   94,
e_sn                                = 120,    97,    91,
e_vert                              = 31,    30,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 9000, 10000,  3333.33,
dy                                  = 9000, 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                          = 8,     3,     3,
ra_lw_physics                       = 4,     1,     1,
ra_sw_physics                       = 4,     1,     1,
radt                                = 10,    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,     1,     0,
cudt                                = 5,     5,     5,
isfflx                              = 1,
ifsnow                              = 0,
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                            = 0,      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,     
/

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


namelist.wps

&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2018-11-29_00:00:00','2006-08-16_12:00:00',
end_date   = '2018-11-30_00:00:00','2006-08-16_12:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/

&geogrid
parent_id         =   1,   1,
parent_grid_ratio =   1,   3,
i_parent_start    =   1,  31,
j_parent_start    =   1,  17,
e_we              =  120, 112,
e_sn              =  120,  97,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! 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','default',
dx = 9000,
dy = 9000,
map_proj = 'lambert',
ref_lat   =  38.6,
ref_lon   =  113.2,
truelat1  =  30.0,
truelat2  =  60.0,
stand_lon =  113.2,
geog_data_path = '/home/user/geog/'
/

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

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

yangshujian

问题已解决，3.6.1的WPS中分辨率设置不同和静态地理资料的差异都有可能导致这种情况

安静的天分

这种情况不会影响WRF的结果吧

yangshujian

安静的天分 发表于 2021-3-9 10:43
这种情况不会影响WRF的结果吧

我做的其实是有影响的，WRF结果也会显示一样的形状，不过也受模拟区域的影响，模拟区域大一些就不会