WRF模式最底层如何设低一点？

zhaolinanshiwo

求问各位大佬：模式最顶层可以设置namelist里的p_top_requested,那我想把模式最底层设置低一点该如何设置呢？

我知道eta_levels可以具体设置垂直分层的每一层，但是最底层eta始终等于1，肯定不是在这里设置的吧

还请知道的大佬们不吝赐教~

tbag

第一直觉是在不改变P top的情况下增加层数，但如此同时肯定会增加计算时间，我记得可以手动去设定eta_level,你就把第二个eta_level弄的接近于1一点不就行了么，比如原来是1，0.995变成1，0.0998

zhaolinanshiwo

tbag 发表于 2017-12-7 03:10
第一直觉是在不改变P top的情况下增加层数，但如此同时肯定会增加计算时间，我记得可以手动去设定eta_level ...

谢谢解答！
可是我比较疑惑的是最底层eta值不是始终为1？这样设置eta_levels只能加密低层吧？

weinihou

建议lz 去看下WRF的模式说明，不是那个用户手册

你设置的eta是w层，staggered level，代码里是叫znw，你现在想要要的是scalar level，znu，是根据znw插值出来的
所以“始终为一”的其实是边界条件

zhaolinanshiwo

weinihou 发表于 2017-12-7 11:26
建议lz 去看下WRF的模式说明，不是那个用户手册

你设置的eta是w层，staggered level，代码里是叫znw，你 ...

答主你说的模式说明是？

weinihou

zhaolinanshiwo 发表于 2017-12-7 20:55
答主你说的模式说明是？

技术说明，官网 Doc/Pub 里能找到
http://202.119.32.195/cache/7/03 ... 0f15f8ab/arw_v3.pdf

zhaolinanshiwo

tbag 发表于 2017-12-7 03:10
第一直觉是在不改变P top的情况下增加层数，但如此同时肯定会增加计算时间，我记得可以手动去设定eta_level ...

你好！我试着手动设置了eta_levels（直接在namelist.input里面添加了这一参数），运行出现了问题，但却没有报错！
您知道手动设置的时候还需要修改其他什么参数么？
namelist如下：

&time_control

run_days                            = 0,

run_hours                           = 30,

run_minutes                         = 0,

run_seconds                         = 0,

start_year                          = 2014, 2014, 2014

start_month                         = 07, 07, 07

start_day                           = 07, 07, 07

start_hour                          = 12, 12, 12

start_minute                        = 00, 00, 00

start_second                        = 00, 00, 00

end_year                            = 2014, 2014, 2014

end_month                           = 07, 07, 07

end_day                             = 08, 08, 08

end_hour                            = 18, 18, 18

end_minute                          = 00, 00, 00

end_second                          = 00, 00, 00

interval_seconds                    = 21600,

input_from_file                     = .true., .true., .true.

history_interval                    = 90, 30, 30

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                           = 135,

time_step_fract_num                 = 0,

time_step_fract_den                 = 1,

max_dom                             = 3,

e_we                                = 65, 91, 124

e_sn                                = 49, 61, 97

e_vert                              = 40, 40, 40

eta_levels                          = 1.000, 0.99999,0.9998, 0.998, 0.997, 0.9954, 0.993, 0.983, 0.970, 0.954, 0.934, 0.909, 0.880,0.845, 0.807, 0.765, 0.719, 0.672, 0.622, 0.571, 0.520, 0.468, 0.420, 0.376,0.335, 0.298, 0.263, 0.231, 0.202, 0.175, 0.150, 0.127, 0.106, 0.088, 0.070,0.055, 0.040, 0.026, 0.013, 0.000

p_top_requested                     = 1000,

num_metgrid_levels                  = 27,

num_metgrid_soil_levels             = 4,

dx                                  = 27000,9000, 3000

dy                                  = 27000,9000, 3000

grid_id                            = 1,     2,    3,

parent_id                           = 1, 1, 2,

i_parent_start                      = 1, 18, 25,

j_parent_start                      = 1, 15, 14,

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                       = 2,     2,    2,

radt                                = 30,    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                              = 1,

icloud                              = 1,

surface_input_source                = 1,

num_soil_layers                     = 4,

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,

/

如下图，wrf就停在这一步也不报错

                               
登录/注册后可看大图

                               
登录/注册后可看大图

小小小鱼

zhaolinanshiwo 发表于 2017-12-10 23:11
你好！我试着手动设置了eta_levels（直接在namelist.input里面添加了这一参数），运行出现了问题，但却没 ...

你好，请问这个问题解决了吗

zhaolinanshiwo

小小小鱼 发表于 2018-11-11 21:11
你好，请问这个问题解决了吗

没有解决，之前被分配学习其他的东西了
现在开始继续学习WRF，互相学习~

chongzika

zhaolinanshiwo 发表于 2018-11-12 10:53
没有解决，之前被分配学习其他的东西了
现在开始继续学习WRF，互相学习~

没看出来有什么问题