Hello,
I am using WRF-4. 5 to get some radiation variables like shortwave flux and I got the variables like ACSWUPT for normal and ACSWUPTC for clear sky shortwave flux at top in wrf output. However, there is no difference in normal and clear sky variables that I got from the output. Are their values supposed to be same in WRF? Or I have to change something in the registry file to get the different output. I am using RRTMG scheme (=4) for sw and lw both. I have also tested with with CAM (radiation scheme=3) for sw and lw physics.Again the result is same. Both output for radiation ASWUPT and ASWUPTC and others radiation values are same. Please clarify this.
Here is my namelist.input:
&time_control
!run_days = 0,
run_hours = 72,
!run_minutes = 0,
!run_seconds = 0,
start_year = 2018, 2018,
start_month = 11, 11,
start_day = 14, 14,
start_hour = 18, 18,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2018, 2018,
end_month = 11, 11,
end_day = 17, 17,
end_hour = 18, 18,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,
history_interval = 180, 60,
frames_per_outfile = 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360,
io_form_auxinput4 = 2,
auxinput1_inname = "met_em.d<domain>.<date>",
force_use_old_data = .true.,
/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 414,667,
e_sn = 484,661,
e_vert = 53, 53,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 9000,3000,
dy = 9000,3000,
grid_id = 1,2,
parent_id = 1,1,
i_parent_start = 1,116,
j_parent_start = 1,50,
parent_grid_ratio = 1,3,
parent_time_step_ratio = 1,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
mp_physics = 10, 10,
ra_lw_physics = 3, 3,
ra_sw_physics = 3, 3,
radt = 9, 9,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 1, 0,
cu_diag = 1,
cudt = 0, 0,
ishallow = 0,
isfflx = 1,
ifsnow = 1,
icloud = 0,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 1, 1
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
do_radar_ref = 1,
lightning_option = 2, 2,
cu_rad_feedback = .true.,
num_land_cat = 24,
sst_update = 1,
/
&fdda
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.01, 0.01,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 2, 2,
scalar_adv_opt = 2, 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,
hybrid_opt = 0,
use_theta_m = 0,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I am using WRF-4. 5 to get some radiation variables like shortwave flux and I got the variables like ACSWUPT for normal and ACSWUPTC for clear sky shortwave flux at top in wrf output. However, there is no difference in normal and clear sky variables that I got from the output. Are their values supposed to be same in WRF? Or I have to change something in the registry file to get the different output. I am using RRTMG scheme (=4) for sw and lw both. I have also tested with with CAM (radiation scheme=3) for sw and lw physics.Again the result is same. Both output for radiation ASWUPT and ASWUPTC and others radiation values are same. Please clarify this.
Here is my namelist.input:
&time_control
!run_days = 0,
run_hours = 72,
!run_minutes = 0,
!run_seconds = 0,
start_year = 2018, 2018,
start_month = 11, 11,
start_day = 14, 14,
start_hour = 18, 18,
start_minute = 00, 00,
start_second = 00, 00,
end_year = 2018, 2018,
end_month = 11, 11,
end_day = 17, 17,
end_hour = 18, 18,
end_minute = 00, 00,
end_second = 00, 00,
interval_seconds = 21600
input_from_file = .true.,.true.,
history_interval = 180, 60,
frames_per_outfile = 1000, 1000,
restart = .false.,
restart_interval = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
auxinput4_inname = "wrflowinp_d<domain>",
auxinput4_interval = 360,
io_form_auxinput4 = 2,
auxinput1_inname = "met_em.d<domain>.<date>",
force_use_old_data = .true.,
/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 2,
e_we = 414,667,
e_sn = 484,661,
e_vert = 53, 53,
p_top_requested = 5000,
num_metgrid_levels = 32,
num_metgrid_soil_levels = 4,
dx = 9000,3000,
dy = 9000,3000,
grid_id = 1,2,
parent_id = 1,1,
i_parent_start = 1,116,
j_parent_start = 1,50,
parent_grid_ratio = 1,3,
parent_time_step_ratio = 1,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
mp_physics = 10, 10,
ra_lw_physics = 3, 3,
ra_sw_physics = 3, 3,
radt = 9, 9,
sf_sfclay_physics = 1, 1,
sf_surface_physics = 2, 2,
bl_pbl_physics = 1, 1,
bldt = 0, 0,
cu_physics = 1, 0,
cu_diag = 1,
cudt = 0, 0,
ishallow = 0,
isfflx = 1,
ifsnow = 1,
icloud = 0,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 1, 1
mp_zero_out = 2,
mp_zero_out_thresh = 1.e-12,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
do_radar_ref = 1,
lightning_option = 2, 2,
cu_rad_feedback = .true.,
num_land_cat = 24,
sst_update = 1,
/
&fdda
/
&dynamics
rk_ord = 3,
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 0, 0,
diff_6th_factor = 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000.,
dampcoef = 0.01, 0.01,
khdif = 0, 0,
kvdif = 0, 0,
non_hydrostatic = .true., .true.,
moist_adv_opt = 2, 2,
scalar_adv_opt = 2, 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,
hybrid_opt = 0,
use_theta_m = 0,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true., .false.,
nested = .false., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/