Problem with calculating direct aerosol-effect to radiation

Topics related to current and future physics in the WRF as well as any problems you may have.

Problem with calculating direct aerosol-effect to radiation

Postby twis91 » Sun Sep 28, 2014 6:22 am

Hey everyone,

my working group is trying to simulate the direct aerosol effect 'cause of dust. We use the "new goddard radiation scheme" and an emission scheme which we developed at our university.

In normal cases, turning on/off the variable aer_ra_feedback should do exact the thing we want? In fact we thought so and tried it. But in both cases we got the same result. In a new test, we tried one simulation without dust and got the same result. As far as I can guess, it seems that the influence of dust toward the radiation is not calculated, the dust instead is calculated. As far as I know, the working group has changed some modules in "phys" and "chem", so the analysis of the error may be complicated and tricky, but I hope you can help (maybe modules to search for or way of testing where the error is). We use V3.51 for WRF.

In the attachment you can see our namelist (with activated feedback).

Kind regards, Robert

P.S. Hope my english is good enough to understand the problem :D
P.P.S. Don't know why I can't attach files, so I put the code here:

+++++++++++++++++++++++++++++++++
&time_control
run_days = 10,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2012,
start_month = 10,
start_day = 09,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2012,
end_month = 10,
end_day = 11,
end_hour = 00,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file = .true.,
fine_input_stream = 0,
history_interval = 180,
frames_per_outfile = 3000,
restart = .false.,
restart_interval = 1440,
input_from_file = .true.,
restart = .false.,
auxinput5_interval_m = 0,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput4 = 0
io_form_auxinput5 = 0
debug_level = 0
/

&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
s_we = 1,
s_sn = 1,
e_we = 320,
e_sn = 180,
s_vert = 1,
e_vert = 39,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 30000,
dy = 30000,
grid_id = 1,
parent_id = 0,
i_parent_start = 0,
j_parent_start = 0,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 1,
use_adaptive_time_step = .false.,
step_to_output_time = .true.,
target_cfl = 1.2,
max_step_increase_pct = 5,
starting_time_step = -1,
max_time_step = -1,
min_time_step = -1,
smooth_option = 0,
p_top_requested = 5000,
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,
/

&physics
mp_physics = 2,
progn = 0,
ra_lw_physics = 5,
ra_sw_physics = 5,
radt = 30,
sf_sfclay_physics = 11,
sf_surface_physics = 3,
bl_pbl_physics = 1,
bldt = 0,
cu_physics = 5,
prec_acc_dt = 180.,
cudt = 0,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 24,
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 = .false.,
cu_diag = 1,
/

&fdda
/

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

&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,
nested = .false.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

&chem
kemit = 1,
chem_opt = 300,
bioemdt = 30,
photdt = 30,
chemdt = 0.5,
io_style_emissions = 1,
emiss_inpt_opt = 1,
emiss_opt = 6,
chem_in_opt = 0,
phot_opt = 0,
gas_drydep_opt = 1,
bio_emiss_opt = 0,
gas_bc_opt = 1,
gas_ic_opt = 1,
aer_bc_opt = 1,
aer_ic_opt = 1,
gaschem_onoff = 1,
aerchem_onoff = 1,
wetscav_onoff = 0,
conv_tr_wetscav = 1,
cldchem_onoff = 0,
vertmix_onoff = 1,
chem_conv_tr = 1,
seas_opt = 0,
bg_dust = 0,
dust_opt = 4,
dust_schme = 4,
dmsemis_opt = 1,
biomass_burn_opt = 1,
plumerisefire_frq = 30,
have_bcs_chem = .false.,
aer_op_opt = 1,
aer_ra_feedback = 1,
/

&end
twis91
 
Posts: 7
Joined: Sun Sep 28, 2014 5:55 am

Re: Problem with calculating direct aerosol-effect to radiat

Postby easy » Thu Nov 19, 2015 4:40 am

Do you find the reason now? I think the option of aer_ra_feedback should have some impact to radiation impact. what do you mean you get the same result, I mean, which variables?

Best,

zhen
easy
 
Posts: 10
Joined: Sun Jun 07, 2015 10:26 pm

Re: Problem with calculating direct aerosol-effect to radiat

Postby mirabelle » Sun Dec 06, 2015 3:24 am

Hi
I think it is because the WRF3.5 doesn't use the dust even it has the qndeop varaiable. only the WRF/Chem can use the aerosol system.
mirabelle
 
Posts: 3
Joined: Sat Dec 05, 2015 11:52 pm


Return to Physics

Who is online

Users browsing this forum: Baidu [Spider] and 3 guests