~~关于TKE的输出~~~~

zhanglang86

各位大神，不知道哪里有问题 TKE输出全部一直是0，请大家帮我看看。问题出在哪里呢？
我的namelist：
&time_control            
run_days                            = 0,
run_hours                           = 12,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2016,
start_month                         = 03,
start_day                           = 19,
start_hour                          = 12,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2016,
end_month                           = 03,
end_day                             = 20,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 10800
input_from_file                     = .true.,.true.,.true.,.true,
history_interval                    = 15,  15,   15, 15,
frames_per_outfile                  = 1000, 1, 1, 1,
restart                             = .false.,
restart_interval         = 1800000,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
io_form_auxinput2                   = 2,
io_form_auxinput4                   = 2,
io_form_auxinput5                   = 0,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 0,
auxinput4_interval_m                = 360,  360,  360,
auxinput5_interval_m                = 60,   60,   60,
auxinput7_interval_m                = 1440, 1440, 1440,
auxinput1_inname                    = "./met_em.d<domain>.<date>",
auxinput6_inname                    = "wrfbiochemi_d<domain>",
auxinput7_inname                    = "wrffirechemi_d<domain>",
auxinput4_inname                    = "wrflowinp_d<domain>",
frames_per_auxinput7                = 38,   38,   38,
debug_level                         = 99,
/

&domains                 
time_step                           = 300,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 86,142,127,124,
e_sn                                = 89,127,115,109,
e_vert                              = 39,    39,    39, 39,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 54000,18000,6000,2000,
dy                                  = 54000,18000,6000,2000,
grid_id                             = 1,2,3,4,
parent_id                           = 1,1,2,3,
i_parent_start                      = 1,15,53,45,
j_parent_start                      = 1,22,55,38,
parent_grid_ratio                   = 1,3,3,3,
parent_time_step_ratio              = 1,3,3,3,
feedback                            = 1,
smooth_option                       = 0
num_traj                 = 0,
/

&physics
mp_physics                          = 2,        2,        2,        2,
progn                               = 1,        1,        1,        1,  
naer                                = 1e9,      1e9,      1e9,      1e9,
ra_lw_physics                       = 5,        5,        5,        5,
ra_sw_physics                       = 5,        5,        5,        5,
radt                                = 30,       30,       30,       30,
sf_sfclay_physics                   = 2,        1,        1,        1,
sf_surface_physics                  = 2,        1,        1,        3,
bl_pbl_physics                      = 2,        2,        2,        2,
mfshconv                            = 1,        1,        0,        0,
bldt                                = 0,        0,        0,        0,
mp_zero_out                         = 2,                           
mp_zero_out_thresh                  = 1.e-8,                        
cu_physics                          = 5,        5,        5,        5,
cudt                                = 1,        1,        1,        1,
cugd_avedx                          = 1,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 6,   
sf_urban_physics                    = 2,        1,        1,        1,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
sst_update                          = 1,
cu_rad_feedback                     = .true.,.true.,  .true.,  .true.,
cu_diag                             = 1, 1, 1, 1,
bl_mynn_tkebudget                   = 1, 1,
/


&fdda
/

&dynamics
rk_ord                              = 3,
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 1,            0,            0,            0,
diff_6th_factor                     = 0.12,         0.12,         0.12,         0.12,
base_temp                           = 290.,                                    
damp_opt                            = 0,                                       
zdamp                               = 5000.,        5000.,        5000.,        5000.,
dampcoef                            = 0.01,         0.01,         0.01,         0.01,
khdif                               = 0,            0,            0,            0,
kvdif                               = 0,            0,            0,            0,
non_hydrostatic                     = .true.,       .true.,       .true.,       .true.,
moist_adv_opt                       = 2,            2,             2,            2,
scalar_adv_opt                      = 2,            2,             2,            2,
chem_adv_opt                        = 2,            2,             2,            2,
tke_adv_opt                         = 2,            2,             2,            2,
/

&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,
/

ppap

楼主问题解决了吗？

ppap

楼主，为什么我运行的wrf的结果里没有tke这个参数？请问您的参数化方案是什么样的？

dophineliu

ppap 发表于 2017-9-9 19:41
楼主，为什么我运行的wrf的结果里没有tke这个参数？请问您的参数化方案是什么样的？

TKE 默认不输出， 要在Registry 里修改为输出， 结果才会输出

dophineliu

楼主问题解决了吗，我的tke输出也都是0

irides

dophineliu 发表于 2018-3-16 08:47
楼主问题解决了吗，我的tke输出也都是0

请问您解决了吗？

green_tea789

你的namelist:
sf_sfclay_physics&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160; = 2,&#160; &#160;&#160; &#160;&#160;&#160;1,&#160; &#160;&#160; &#160;&#160;&#160;1,&#160; &#160;&#160; &#160;&#160;&#160;1,
bl_pbl_physics&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160;&#160; &#160; = 2,&#160; &#160;&#160; &#160;&#160;&#160;2,&#160; &#160;&#160; &#160;&#160;&#160;2,&#160; &#160;&#160; &#160;&#160;&#160;2,
手册里写：
Mellor-Yamada-Janjic (Eta) TKE scheme; use sf_sfclay_physics=2     
可能是因为sf_sfclay_physics用的是1

green_tea789

变量输出为0，要去看namelist.input里参数化方案用的是什么，这个参数化方案计算输出这个变量吗？ 如果输出还是为0，就要去查手册，是否要配套其他的参数化方案。