更改地形后wrf中途崩，未发现错误

鹿先森

进行地形敏感性实验，之前已经尝试过改一小块（wrf正常运行）直接改D02抹平（wrf正常运行）
改地形的方法：直接更改geo_em**或者met_em**的地形高度，wrf均可运行（插一句，这俩结果一样）
题来了，这次我尝试只改一个山脉，范围比上一个小区域大，D02小。但是wrf一转就崩。rsl.erro日.0000没问题，real.exe的时候rsl.erro.0000也没问题。
改地形的方法：直接更改geo_em**或者met_em**的地形高度，两者均尝试，都崩。
自己尝试找问题：
1、重装WRF,换原始的geo_em**，直接改D02抹平的geo_em**，WRF运行正常
2、使用改山脉的geo_em**，换不同的参数化方案，WRF崩
3、改山脉geo_em**高度为1500，2000，WRF崩
4、怀疑服务器原因换了一个服务器，WRF崩

我实在是不知道怎么回事了～～～～明明改地形的操作步骤一样～
跪求知道的大神解惑，实在是困惑了好几天～～先谢谢各位前辈了！

ps:
改高度脚本
a=addfile("./geo_em.d01.nc","w")
hgt = a->HGT_M
hgt(0,180:182,196:199) = 1500
a->HGT_M = hgt
namelist.input

&time_control
run_days                            = 0,
run_hours                           = 18,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017,
start_month                         = 05,   05,
start_day                           = 9,   9,
start_hour                          = 00,   00,
start_minute                        = 00,   00,
start_second                        = 00,   00,
end_year                            = 2017, 2017,
end_month                           = 05,   05,
end_day                             = 9,   9,
end_hour                            = 18,   18,
end_minute                          = 00,   00,
end_second                          = 00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 180,  30,
frames_per_outfile                  = 1, 1,
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                           = 18,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 367,    550,
e_sn                                = 367,    550,
e_vert                              = 30,    30,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 28,
num_metgrid_soil_levels             = 4,
dx                                  = 6000, 2000,
dy                                  = 6000, 2000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     92,
j_parent_start                      = 1,     92,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

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

rsl.error.0000
taskid: 0 hostname: c0214
module_io_quilt_old.F        2931 F
Quilting with   1 groups of   0 I/O tasks.
Ntasks in X            8, ntasks in Y            8
--- WARNING: traj_opt is zero, but num_traj is not zero; setting num_traj to zero.
--- NOTE: （这部分帖子，提交说有非法字符）
--- NOTE: RRTMG radiation is not used, setting:  o3input=0 to avoid data pre-processing
--- NOTE: num_soil_layers has been set to      4
WRF V3.8.1 MODEL
Parent domain
ids,ide,jds,jde            1         367           1         367
ims,ime,jms,jme           -4          53          -4          53
ips,ipe,jps,jpe            1          46           1          46
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1,               54174836 bytes allocated
  med_initialdata_input: calling input_input
Timing for processing wrfinput file (stream 0) for domain        1:   13.73938 elapsed seconds
Max map factor in domain 1 =  1.02. Scale the dt in the model accordingly.
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           17 SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS


LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19 CATEGORIES
Nesting domain
ids,ide,jds,jde            1         550           1         550
ims,ime,jms,jme           -4          81          -4          81
ips,ipe,jps,jpe            1          69           1          69
INTERMEDIATE domain
ids,ide,jds,jde           90         278          90         278
ims,ime,jms,jme           85         124          85         124
ips,ipe,jps,jpe           88         114          88         114
*************************************
d01 2017-05-09_00:00:00  alloc_space_field: domain            2,                7699200 bytes allocated
d01 2017-05-09_00:00:00  alloc_space_field: domain            2,              116002564 bytes allocated
d01 2017-05-09_00:00:00 *** Initializing nest domain # 2 from an input file. ***
d01 2017-05-09_00:00:00 med_initialdata_input: calling input_input
Timing for processing wrfinput file (stream 0) for domain        2:   35.09691 elapsed seconds

INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           17 SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19 CATEGORIES
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           17 SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS
LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19 CATEGORIES
Max map factor in domain 1 =  1.02. Scale the dt in the model accordingly.
INPUT LandUse = "MODIFIED_IGBP_MODIS_NOAH"
LANDUSE TYPE = "MODIFIED_IGBP_MODIS_NOAH" FOUND          33 CATEGORIES           2 SEASONS WATER CATEGORY =           17 SNOW CATEGORY =           15
INITIALIZE THREE Noah LSM RELATED TABLES
Skipping over LUTYPE = USGS

LANDUSE TYPE = MODIFIED_IGBP_MODIS_NOAH FOUND          20 CATEGORIES
INPUT SOIL TEXTURE CLASSIFICATION = STAS
SOIL TEXTURE CLASSIFICATION = STAS FOUND          19 CATEGORIES
mediation_integrate.b        1724 DATASET=HISTORY
mediation_integrate.b        1725  grid%id            1  grid%oid            1
Timing for Writing wrfout_d01_2017-05-09_00:00:00 for domain        1:   10.90402 elapsed seconds
Timing for processing lateral boundary for domain        1:    2.38494 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     46 JS      1 JE     46
WRF NUMBER OF TILES =   1
mediation_integrate.b        1724 DATASET=HISTORY
mediation_integrate.b        1725  grid%id            2  grid%oid            2
Timing for Writing wrfout_d02_2017-05-09_00:00:00 for domain        2:   27.23703 elapsed seconds
Tile Strategy is not specified. Assuming 1D-Y
WRF TILE   1 IS      1 IE     69 JS      1 JE     69
WRF NUMBER OF TILES =   1
Timing for main: time 2017-05-09_00:01:00 on domain   2:   34.08164 elapsed seconds
Timing for main: time 2017-05-09_00:02:00 on domain   2:    5.71025 elapsed seconds

tbag

surface_input_source                = 1,
这个改到3试一下
平滑的问题我觉得也很重要

werewolf

地形梯度太大了吧，把改后的wrfinput画个图看看

鹿先森

werewolf 发表于 2020-07-23 23:37
地形梯度太大了吧，把改后的wrfinput画个图看看

改完地形高度后的，geo的地形图

                               
登录/注册后可看大图

来自苹果APP客户端

鹿先森

werewolf 发表于 2020-07-23 23:37
地形梯度太大了吧，把改后的wrfinput画个图看看

wrfout的我明天画了再补上来，谢谢您！

来自苹果APP客户端

werewolf

鹿先森 发表于 2020-7-24 00:04
wrfout的我明天画了再补上来，谢谢您！

不用了，地形是静态数据，不会变的，都一样。这个感觉是地形梯度太大了，计算气压梯度等可能会出错，最好做个平滑

鹿先森

werewolf 发表于 2020-07-24 00:14
不用了，地形是静态数据，不会变的，都一样。这个感觉是地形梯度太大了，计算气压梯度等可能会出错，最好做个平滑

                               
登录/注册后可看大图

                               
登录/注册后可看大图

好奇的问一下，我之前直接把整个D02抹成1500，WRF也可以转结果，那个时候青藏高原东北边的地形落差也很大～～但是WRF也能运行啊～

来自苹果APP客户端

islandowner

鹿先森 发表于 2020-7-24 00:20
好奇的问一下，我之前直接把整个D02抹成1500，WRF也可以转结果，那个时候青藏高原东北边的地形落差也很大 ...

感觉和楼上说的，你这地形落差太大，做个平滑（估计得多平滑几次，有专门平滑地形的函数）。至于你第一次直接抹掉D02地形，落差也大，为什么能运行？个人猜测应该是这样，你这个落差大是D01和D02之间的，对D02边界的影响虽然比较大，但是D01传给D02的边界条件会加一些“海绵区”，使得过渡更平稳。

HANSEN

请问楼主解决了吗？遇到了同样的问题 请问楼主如果解决了的话 是怎么操作的？谢谢！

鹿先森

HANSEN 发表于 2022-7-28 20:27
请问楼主解决了吗？遇到了同样的问题 请问楼主如果解决了的话 是怎么操作的？谢谢！

由于地形落差太大无法运行，后期放弃了，楼上提供的平滑思路你可以尝试一下，加油！