Hi,
I have observed big differences in all-sky surface and TOA radiative forcings between WRF only simulation and WRF-Chem simulation with aer_rad_feedback=0. Differences are close to 10-12 Wm-2. However clear sky radiative forcings are similar which implies difference in cloud simulation. What could be the possible reason for this difference as otherwise both simulation use the same physics options? Are clouds simulated slightly differently with Chem than in WRF only simulation? Shouldn't such pair of simulations provide similar model outputs?
I have observed big differences in all-sky surface and TOA radiative forcings between WRF only simulation and WRF-Chem simulation with aer_rad_feedback=0. Differences are close to 10-12 Wm-2. However clear sky radiative forcings are similar which implies difference in cloud simulation. What could be the possible reason for this difference as otherwise both simulation use the same physics options? Are clouds simulated slightly differently with Chem than in WRF only simulation? Shouldn't such pair of simulations provide similar model outputs?
Code:
&physics
mp_physics = 10,10,10,10,
progn = 1,1,1,1,
ra_lw_physics = 4,4,4,4,
ra_sw_physics = 4,4,4,4,
radt = 10,10,10,10,
sf_sfclay_physics = 5,5,5,5,
sf_surface_physics = 2,2,2,2,
bl_pbl_physics = 6,6,6,6,
bldt = 0,0,0,0,
cu_physics = 3,3,3,0,
shcu_physics = 0,0,0,0,
cudt = 0,0,0,0,
cugd_avedx = 1,
ishallow = 0,
cu_diag = 1,1,1,0,
shcu_aerosols_opt = 2,2,2,2,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 3,
num_land_cat = 21,
sf_urban_physics = 0,0,0,0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .true.,.true,.false.,.false.,
sst_update = 1,
sst_skin = 1,
usemonalb = .true.,
rdmaxalb = .true.,
rdlai2d = .true.,
slope_rad = 1,1,1,1,
topo_shading = 1,1,1,1,
shadlen = 25000,
/