登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我想模拟台风的实际个例,没有加bogus,模拟2006年的珍珠,起始时刻设置为2006年5月9日00时,终止时刻为5月19日00时,大概在5月16日左右台风中心会开始向西偏移最终离开domain,想请教一下这个是什么原因造成的啊。 初始时刻设置的domian如图,初始时刻JTWC定义的中心在(133.7E,8.0N)
namelist.input如下:&time_control
run_days = 10,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2006, 2006, 2006, 2006,
start_month = 05, 05, 05, 05,
start_day = 09, 09, 09, 09,
start_hour = 00, 00, 00, 00,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2006, 2006, 2006, 2006,
end_month = 05, 05, 05, 05,
end_day = 19, 19, 19, 19,
end_hour = 00, 00, 00, 00,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 21600,
input_from_file = .true.,.false.,.false.,.false.,
history_interval = 180, 180, 15, 15,
frames_per_outfile = 1, 1, 1, 1,
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2io_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 = 4,
s_we = 1, 1, 1, 1,
e_we = 121, 181, 301, 301,
s_sn = 1, 1, 1, 1,
e_sn = 121, 181, 301, 301,
s_vert = 1, 1, 1, 1,
e_vert = 28, 28, 28, 28,
num_metgrid_levels = 27,
dx =54000, 18000, 6000, 2000,
dy =54000, 18000, 6000, 2000,
grid_id = 1, 2, 3, 4,
parent_id = 0, 1, 2, 3,
i_parent_start = 1, 31, 41, 181,
j_parent_start = 1, 31, 41, 31,
parent_grid_ratio = 1, 3, 3, 3,
parent_time_step_ratio = 1, 3, 3, 3,
feedback = 1,
smooth_option = 1,
vortex_interval = 15, 15, 15, 15,
max_vortex_speed = 40, 40, 40, 40,
corral_dist = 28, 28, 18, 18,
track_level = 85000,
time_to_move = 720,720,180,180,
/
&physics
mp_physics = 2, 2, 2, 2,
ra_lw_physics = 1, 1, 1, 1,
ra_sw_physics = 1, 1, 1, 1,
radt = 54, 18, 6, 2,
sf_sfclay_physics = 2, 2, 2, 2,
sf_surface_physics = 1, 1, 1, 1,
bl_pbl_physics = 2, 2, 2, 2,
bldt = 0, 0, 0, 0,
cu_physics = 1, 1, 1, 0
cudt = 5, 5, 5, 5,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&fdda
grid_fdda = 1
gfdda_inname = "wrffdda_d01"
gfdda_interval_m = 14400
gfdda_end_h = 240
io_form_gfdda = 2
fgdt = 1
if_no_pbl_nudging_uv = 0
if_no_pbl_nudging_t = 1
if_no_pbl_nudging_q = 1
if_zfac_uv = 1
k_zfac_uv = 10
if_zfac_t = 1
k_zfac_t = 10
if_zfac_q = 1
k_zfac_q = 10
guv = 0.0001
gt = 0.0001
gq = 0.0001
if_ramping = 0
dtramp_min =60.0
/
&dynamics
w_damping = 0,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0, 0,0
diff_6th_factor = 0.12,0.12,0.12,0.12
base_temp = 290.
use_baseparam_fr_nml = .t.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,5000.,
dampcoef = 0.2, 0.2, 0.2, 0.2
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
non_hydrostatic = .true., .true., .true.,.true.,
moist_adv_opt = 1, 1, 1, 1,
scalar_adv_opt = 1, 1, 1, 1,
/k_zfac_q = 10
guv = 0.0001
gt = 0.0001
gq = 0.0001
if_ramping = 0
dtramp_min =60.0
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,.false.,.false.,
nested = .false.,.true., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
|