登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我需要用WRF模拟台风个例,编译的时候选择了自动移动网格(vortex following),但是在编辑namelist的时候没有看到控制移动网格的参数……然后我直接在没有控制移动网格的参数的情况下提交,运行的时候,rsl.out.0000里面每计算一次风暴中心就会有moad can not move字样,应该网格是移动失败的……于是我在网上查到与移动网格有关的参数有这样的:
****移动网格控制参数说明:******** 控制移动嵌套网格的方式有两种:1. 用户指定移动网格,此方式 是指嵌套网格的每一次移动都是由用户通过参数来指定;2. 自动移动网格,此方式是 指移动嵌套网格在启动后,根据模式的计算状态,自动判断下一次的移动参数。 用户指定移动:编译时需要在 ARCHFLAGS 选项中添加“-D MOVE_NESTS”来激活。允许的 最大移动套网格移动次数为 50,不过也可以在源程序 frame/module_driver_constants.F 进行修改。) ) num_moves = 4 移动嵌套网格总移动次数。 move_id = 2,2,2,2, 每一次移动嵌套网格区域编号列表。 move_interval = 60,120,150,180, 每一次移动的启动时间列表,单位为分钟,自模式积分起始时刻算起。 move_cd_x = 1,1,0,-1, 在 i 方向(即东西方向)每一次相对于父网格移动格点数。 move_cd_y = 1,0,-1,1, 在 j 方向(即南北方向)每一次相对于父网格移动格点数。 正整数表示顺着 i/j 值增大的方向,负值表示顺着 i/j 值减小的方向。0 表示不移动。 目前移动距离限制只能为一个网格单元。 自动移动: 编译时需要在 ARCHFLAGS 选项中添加 “-D MOVE_NESTS”和“-DVORTEX_CENTER”来激活。 目前,这些参数是应用中等涡旋追随法(mid-level vortex following algorithm)来确定嵌套网格的移动,还在测试阶段。 vortex_interval = 15 经过多长时间计算一次涡旋的位置,单位为分钟 max_vortex_speed = 40 涡旋的最大移动速度,用于计算新涡旋位置的搜索半径 corral_dist =8 移动嵌套网格靠近粗网格边界允许的最大网格单元数, 此参数也就是规定了移动网格靠 近粗网格允许的最大距离。 然后我查看了configure.wrf的ARCHFLAGS的选项里面有“-D MOVE_NESTS”和“-DVORTEX_CENTER”,那我就试着在namelist.input的&domain里面添加了vortex_interval max_vortex_speed corral_dist这些参数,结果./real.exe的时候说读取domain信息失败…… 请大神教教我怎么解决这个问题!~~ PS:大神们有什么需要看的与问题有关的数据尽管告知!~我把namelist.input贴上来:&time_control
run_days = 10,
run_hours = 18,
run_minutes = 0,
run_seconds = 0,
start_year = 2013, 2013, 2000,
start_month = 07, 07, 01,
start_day = 05, 05, 24,
start_hour = 00, 00, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2013, 2013, 2000,
end_month = 07, 07, 01,
end_day = 15, 15, 25,
end_hour = 18, 18, 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 = 90,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 184, 388, 94,
e_sn = 73, 112, 91,
e_vert = 30, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 30000, 10000, 3333.33,
dy = 30000, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 19, 30,
j_parent_start = 1, 19, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/ &physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 10, 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,
sf_urban_physics = 0, 0, 0,
/ &fdda
/ &dynamics
w_damping = 1,
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, 以及configure.wrf里面的ARCHFLAGSARCHFLAGS = $(COREDEFS) -DIWORDSIZE=$(IWORDSIZE) -DDWORDSIZE=$(DWORDSIZE) -DRWORDSIZE=$(RWORDSIZE) -DLWORDSIZE=$(LWORDSIZE) \
$(ARCH_LOCAL) \
$(DA_ARCHFLAGS) \
-DDM_PARALLEL \
\
-DNETCDF \
\
\
\
\
\
\
\
-DMOVE_NESTS -DVORTEX_CENTER \
-DUSE_ALLOCATABLES \
-DGRIB1 \
-DINTIO \
-DLIMIT_ARGS \
-DCONFIG_BUF_LEN=$(CONFIG_BUF_LEN) \
-DMAX_DOMAINS_F=$(MAX_DOMAINS) \
-DMAX_HISTORY=$(MAX_HISTORY) \
|