求助：wrf-chem怎么做沙尘气溶胶的辐射反馈？

chickai1109

我用的长波和短波辐射方案都是rrtmg，namelist.input里：aer_ra_feedback                     = 1, aer_op_opt                          = 1,
opt_pars_out                        = 1,
得出的结果和这三个等于0的结果基本上没啥差别，和WRF直接模拟的结果也差不多。。。
请问辐射反馈namelist除了这3项还要设那些？
这是namalist：
&time_control
run_days                            = 6,
run_hours                           = 18,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2019, 2019, 2000,
start_month                         = 07,   07,   01,
start_day                           = 23,   23,   24,
start_hour                          = 00,   00,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2019, 2019, 2000,
end_month                           = 07,   07,   01,
end_day                             = 29,   29,   25,
end_hour                            = 18,   18,   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                  = 1000, 1000, 1000,
restart                             = .false.,
restart_interval                    = 10000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2

auxinput6_inname                    = 'wrfbiochemi_d01',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
auxinput12_inname                   = 'wrf_chem_input',

auxinput5_interval_m                = 21600, 21600, 1440,
auxinput7_interval_m                = 21600, 21600, 1440,
auxinput8_interval_m                = 21600, 21600, 1440,

io_form_auxinput2                   = 2,
io_form_auxinput5                   = 2,
io_form_auxinput6                   = 0,
io_form_auxinput7                   = 2,
io_form_auxinput8                   = 2,
io_form_auxinput12                  = 0,

frames_per_auxinput5  = 24, 24, 24,
frames_per_auxinput6  = 24, 24, 24,
frames_per_auxinput7  = 24, 24, 24,
frames_per_auxinput8  = 24, 24, 24,
frames_per_auxinput12 = 24, 24, 24,

debug_level                         = 0,
auxinput1_inname                    = "met_em.d<domain>.<date>",

/

&dfi_control
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 160,   247,    94,
e_sn                                = 150,   163,    91,
e_vert                              = 48,    48,     30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 27000, 9000,  3333.33,
dy                                  = 27000, 9000,  3333.33,
grid_id                             = 1,     2,     3,
parent_id                           = 1,     1,     2,
i_parent_start                      = 1,     30,    30,
j_parent_start                      = 1,     45,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
zap_close_levels                    = 50,
interp_type                         = 1,
t_extrap_type                       = 2,
force_sfc_in_vinterp                = 0,
use_levels_below_ground             = .true.,
use_surface                         = .true.,
lagrange_order                      = 1,
/
sfcp_to_sfcp                        = .true.,

&physics
num_land_cat                        = 21,
mp_physics                          = 10,    10,    10,
progn                               = 0,     0,     0,
ra_lw_physics                       = 4,     4,     4,
ra_sw_physics                       = 4,     4,     4,
radt                                = 27,    9,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 5,     5,     5,
cu_diag                             = 1,     1,     1,
cudt                                = 0,     0,     0,
ishallow                            = 0,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .true.,
/


&fdda
/

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



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

&chem
kemit                               = 1,
chem_opt                            = 301,      301,      301,
bioemdt                             = 30,       30,        0,
photdt                              = 30,       30,        0,
chemdt                              = 0.1,      0.1,        1,
io_style_emissions                  = 2,
emiss_opt                           = 5,        5,        3,
emiss_opt_vol                       = 0,        0,        0,
emiss_ash_hgt                       = 20000.,
emiss_inpt_opt                      = 1,    1,
chem_in_opt                         = 0,        0,        0,
phot_opt                            = 2,        2,        0,
gas_drydep_opt                      = 0,        0,        0,
aer_drydep_opt                      = 1,        1,        1,
bio_emiss_opt                       = 0,        0,        0,
ne_area                             = 0,
dust_opt                            = 1,

dmsemis_opt                         = 0,
seas_opt                            = 0,
depo_fact                           = 0.25,     0.25,     0.25,
gas_bc_opt                          = 0,        0,        0,
gas_ic_opt                          = 0,        0,        0,
aer_bc_opt                          = 1,        1,        1,
aer_ic_opt                          = 1,        1,        1,
gaschem_onoff                       = 0,        0,        0,
aerchem_onoff                       = 1,        1,        1,
wetscav_onoff                       = 0,        0,        0,
cldchem_onoff                       = 0,        0,        0,
vertmix_onoff                       = 1,        1,        1,
chem_conv_tr                        = 1,        1,        1,
conv_tr_wetscav                     = 0,        0,        0,
conv_tr_aqchem                      = 0,        0,        0,
biomass_burn_opt                    = 0,        0,        0,
plumerisefire_frq                   = 30,       30,       30,
have_bcs_chem                       = .false.,  .false.,  .false.,
aer_ra_feedback                     = 1,
aer_op_opt                          = 1,
opt_pars_out                        = 1,
diagnostic_chem                     = 0,
/

森森小星星

请问你的dustdrydep是零吗？

chickai1109

森森小星星 发表于 2020-4-15 09:57
请问你的dustdrydep是零吗？

是的，这个对辐射反馈有影响么？

森森小星星

chickai1109 发表于 2020-4-15 15:49
是的，这个对辐射反馈有影响么？

不好意思啊 我对辐射不是很清楚 我比较关注沉降
之前我跑出来的dustdrydep都是零 所以在找原因
但是我又用301做了一次发现我的结果不是dustdrydep不是零哦 会不会你还有其他错误

chickai1109

森森小星星 发表于 2020-4-15 21:34
不好意思啊 我对辐射不是很清楚 我比较关注沉降
之前我跑出来的dustdrydep都是零 所以在找原因
但是我 ...

能方便看下你的namelist的chem部分吗？我现在干沉降也出不来

森森小星星

chickai1109 发表于 2020-4-17 21:18
能方便看下你的namelist的chem部分吗？我现在干沉降也出不来

我设置的chem_opt=301
gas_drydep_opt =1
aer_drydep_opt =1

爱橙子不爱橘子

麻烦问一下楼主，做排放源吗？我没做排放源就会出现三个error。

lss19980205

气溶胶辐射反馈是不是只关aer_ra_feedback就行了，aer_op_opt和opt_pars_out不用关吧？