关于WRF-Chem中沉降通量问题

huashenbing

WRF-Chem中选用RADM2_MADE/SORG_KPP模拟，dry_dep_flux为0 ，查了很久的资料也开了部分参数，不知道哪里出列问题，有人遇到过吗？怎么解决？以下是namelist

&time_control
run_days                            = 0,
run_hours                           = 6,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2019, 2019, 2999,
start_month                         = 08,   03,   06,
start_day                           = 21,   18,   11,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2019, 2019, 2999,
end_month                           = 08,   03,   06,
end_day                             = 21,   21,   12,
end_hour                            = 06,   00,   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                  = 1, 1, 1000,
restart                             = .false.,
restart_interval                    = 5000,
auxinput5_interval_m                = 60, 60, 60
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput4                   = 0
io_form_auxinput5                   = 2
auxinput5_inname = '/home/ncar_fdda/hua/WRF-Chem/WRF-TEST/run/wrfchemi_{00z,12z}_d<domain>',
debug_level                         = 00
iofields_filename                   = "myoutfields_d01.txt","myoutfields_d02.txt"
ignore_iofields_warning             = .false.,
/

&domains
time_step                           = 30,
max_dom                             = 1,
e_we                                = 540,   763,   94,
e_sn                                = 478,   613,    91,
e_vert                              = 30,    30,    20,
eta_levels                          =
          1.000000, 0.99300, 0.985, 0.9750, 0.950,0.932347, 0.918907,0.904075, 0.887721, 0.869715, 0.849928, 0.828211,
          0.804436, 0.778472, 0.750192, 0.719474, 0.686214,
          0.650339, 0.611803, 0.570656, 0.526958, 0.480854,
          0.432582, 0.382474, 0.330973, 0.278674, 0.226390,
          0.175086, 0.092183, 0.00000,
dx                                  = 9000, 3000, 6666.66,
dy                                  = 9000, 3000, 6666.66,
p_top_requested                     = 5000,
num_metgrid_levels                  = 34,
num_metgrid_soil_levels             = 4,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     78,    30,
j_parent_start                      = 1,     91,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

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

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

&chem
kemit                               = 4,
chem_opt                            = 106,     401,
bioemdt                             = 30,       0,
photdt                              = 30,       0,
chemdt                              = 30.,     10.,
io_style_emissions                  = 1,        
emiss_inpt_opt                      = 102,        0,
emiss_opt                           = 3,        0,
emiss_opt_vol                       = 0,        0,  
chem_in_opt                         = 0,        0,
phot_opt                            = 2,        0,
gas_drydep_opt                      = 1,        0,
aer_drydep_opt                      = 1,        1,
bio_emiss_opt                       = 1,        0,
dust_opt                            = 0,        
dust_schme                          = 0,
dmsemis_opt                         = 0,        
seas_opt                            = 1,        
gas_bc_opt                          = 1,        0,
gas_ic_opt                          = 1,        0,
aer_bc_opt                          = 1,        1,
aer_ic_opt                          = 1,        1,
gaschem_onoff                       = 1,        0,
aerchem_onoff                       = 1,        1,
wetscav_onoff                       = 0,        0,
cldchem_onoff                       = 1,        0,
vertmix_onoff                       = 1,        1,
chem_conv_tr                        = 1,        0,
biomass_burn_opt                    = 0,        0,
plumerisefire_frq                   = 30,       30,
aer_ra_feedback                     = 0,        0,
have_bcs_chem                       = .false., .false.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

zhumomo

FYI
http://bbs.06climate.com/forum.php?mod=viewthread&tid=57557

huashenbing

zhumomo 发表于 2019-9-25 11:43
FYI
http://bbs.06climate.com/forum.php?mod=viewthread&tid=57557

按照这个设置了，也是不可以

huasheng

你好请问后面有解决么

chickai1109

同问，解决了么？

xyan88

同问，解决了么？
干沉降DRY_DEP_LEN 也都是0

一大碗年糕

xyan88 发表于 2021-8-27 21:23
同问，解决了么？
干沉降DRY_DEP_LEN 也都是0

同问楼主，解决了么？干沉降DRY_DEP_LEN 也都是0，DRYDEPVEL也是0，是不是必须是特定的方案才行？

一大碗年糕

用问请问楼主解决了吗？