Idealized case对称性问题

妄想者

我正在尝试在WRF中实现一个新的表面拖拽方法，所以我想检查x和y方向的对称性。

我的出发点是默认的WRF4.1.2，2d理想化的钟形山体案例（WRF/test/em_hill2d_x）。我把默认的namelist.input改成了更高的分辨率，并修改了dyn_em/module_initialize_ideal.f、namelist.input、input_sounding中的相关部分代码，并在y方向上重新跑了一遍模式（之前是x方向）然而结果并不相同。

然后，我又试着测试了一下上面看到的不对称性，我把运行方式从2D改成了3D。山体仍然是钟形，在x和y方向上完全对称。我的两个运行的唯一区别是输入文件，风速保持在10m/s，所有层的potential temperature都是273K，唯一的区别是风向是西风/南风。附上的图片是用ncview绘制的，可以看到两张图的运行情况明显不对称（左图中画的是的U场，右图中画的是V场，colorbar一样）。

附上了图片和namelist.input，是3D的namelist,边界设置为open。

&time_control
run_days                            = 0,
run_hours                           = 0,
run_minutes                         = 5,
run_seconds                         = 0,
start_year                          = 0001,
start_month                         = 01,
start_day                           = 01,
start_hour                          = 00,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 0001,
end_month                           = 01,
end_day                             = 01,
end_hour                            = 10,
end_minute                          = 00,
end_second                          = 00,
history_interval_s                  = 300,
frames_per_outfile                  = 1000,
restart                             = .false.,
restart_interval                    = 1440,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
/

&domains
time_step                           = 0,
time_step_fract_num                 = 7,
time_step_fract_den                 = 100,
max_dom                             = 1,
s_we                                = 1,
e_we                                = 242,
s_sn                                = 1,
e_sn                                = 242,
s_vert                              = 1,
e_vert                              = 251,
dx                                  = 40,
dy                                  = 40,
ztop                                = 5000.,
/

&physics
mp_physics                          = 0,
ra_lw_physics                       = 0,
ra_sw_physics                       = 0,
radt                                = 0,
sf_sfclay_physics                   = 0,
sf_surface_physics                  = 0,
bl_pbl_physics                      = 0,
bldt                                = 0,
cu_physics                          = 0,
cudt                                = 0,
norm_stress                         = 0,     1,     1,
isfflx                              = 1
/

&fdda
/

&dynamics
hybrid_opt                          = 2,
rk_ord                              = 3,
diff_opt                            = 2,      2,      2,
km_opt                              = 1,      1,      1,
damp_opt                            = 2,
zdamp                               = 2000.,
dampcoef                            = 0.2,
khdif                               = 100,
kvdif                               = 100,
smdiv                               = 0.1,
emdiv                               = 0.01,
epssm                               = 0.4,
mix_isotropic                       = 1,      1,      1,
time_step_sound                     = 6,
h_mom_adv_order                     = 5,
v_mom_adv_order                     = 3,
h_sca_adv_order                     = 5,
v_sca_adv_order                     = 3,
non_hydrostatic                     = .true.,

/

&bdy_control
periodic_x                          = .false.,
symmetric_xs                        = .false.,
symmetric_xe                        = .false.,
open_xs                             = .true.,
open_xe                             = .true.,
periodic_y                          = .false.,
symmetric_ys                        = .false.,
symmetric_ye                        = .false.,
open_ys                             = .true.,
open_ye                             = .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&ideal
ideal_case = 1
/