Hi Everyone,
I was wondering if you could help me with linking Carbon Tracker emissions to WRF-Chem v3.9.1.1 tracers. I have the wrfchemi emissions files and they are being read in without error (I’ve checked the rsl.out and rsl.errror files). I have used the tracers in the chem_trace2 package in the registry as:
The emissions are setup as:
With the package:
However, when In check my wrfout files, while the E_TRA1 - E_TRA10 variables have non-zero values, there tracers are not impacted by them and seem to run without the tracers.
I’ve also attached my nameless file for the run, I’m not sure if I made errors there.
I’m really stuck at this problem, as I’ve not been able to figure out how to link emissions to the tracers. Any advice would be very helpful.
Thanks,
Arkayan
namelist file:
I was wondering if you could help me with linking Carbon Tracker emissions to WRF-Chem v3.9.1.1 tracers. I have the wrfchemi emissions files and they are being read in without error (I’ve checked the rsl.out and rsl.errror files). I have used the tracers in the chem_trace2 package in the registry as:
Code:
package chem_trace2 chem_opt==14 - chem:tracer_1,tracer_2,tracer_3,tracer_4,tracer_5,tracer_6,tracer_7,tracer_8,tracer_9,tracer_10The emissions are setup as:
Code:
# Anthropogenic emissions
#
state real - i+jf emis_ant - - - - "Anthropogenic Emissions" ""
state real e_tra1 i+jf emis_ant 1 Z i5h "E_TRA1" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra2 i+jf emis_ant 1 Z i5h "E_TRA2" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra3 i+jf emis_ant 1 Z i5h "E_TRA3" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra4 i+jf emis_ant 1 Z i5h "E_TRA4" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra5 i+jf emis_ant 1 Z i5h "E_TRA5" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra6 i+jf emis_ant 1 Z i5h "E_TRA6" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra7 i+jf emis_ant 1 Z i5h "E_TRA7" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra8 i+jf emis_ant 1 Z i5h "E_TRA8" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra9 i+jf emis_ant 1 Z i5h "E_TRA9" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1"
state real e_tra10 i+jf emis_ant 1 Z i5h "E_TRA10" "CO2 EMISSIONS (Sib Crop input)" "mol km^-2 hr^-1”With the package:
Code:
package sibco2 emiss_opt==12 - emis_ant:e_tra1,e_tra2,e_tra3,e_tra4,e_tra5,e_tra6,e_tra7,e_tra8,e_tra9,e_tra10I’ve also attached my nameless file for the run, I’m not sure if I made errors there.
I’m really stuck at this problem, as I’ve not been able to figure out how to link emissions to the tracers. Any advice would be very helpful.
Thanks,
Arkayan
namelist file:
Code:
&time_control
run_days = 05,
run_hours = 00,
run_minutes = 00,
run_seconds = 00,
start_year = 2016, 2016, 2016, 2016,
start_month = 07, 07, 07, 07,
start_day = 29, 29, 29, 29,
start_hour = 12, 12, 12, 12,
start_minute = 00, 00, 00, 00,
start_second = 00, 00, 00, 00,
end_year = 2016, 2016, 2a016, 2016,
end_month = 08, 08, 08, 08,
end_day = 03, 03, 03, 03,
end_hour = 12, 12, 12, 12,
end_minute = 00, 00, 00, 00,
end_second = 00, 00, 00, 00,
interval_seconds = 21600
history_interval = 60, 60, 2, 2,
frames_per_outfile = 1, 1, 1, 1,
input_from_file = .true., .true., .true., .true.,
fine_input_stream = 0, 0, 0, 0,
auxinput1_inname = 'met_em.d<domain>.<date>',
restart = .false.,
!restart_interval = 360,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
auxinput12_inname = 'wrf_chem_input',
io_form_auxinput12 = 2,
io_form_auxinput4 = 2,
auxinput4_inname = 'wrflowinp_d<domain>',
auxinput4_interval = 360, 360, 360, 360
auxinput5_inname = 'wrfchemi_d<domain>_<date>'
auxinput5_interval_m = 60, 60, 60, 60
io_form_auxinput5 = 2,
debug_level = 300,
io_form_auxinput2 = 2
/
&domains
time_step = 40,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
s_we = 1, 1, 1,
e_we = 250, 283, 274,
s_sn = 1, 1, 1,
e_sn = 185, 226, 235,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
e_vert = 51, 51, 51, 51,
eta_levels = 1.0000,
0.9975,0.9937,0.9884,0.9821,0.9758,0.9695,0.9632,
0.9568,0.9505,0.9442,0.9368,0.9295,0.9221,0.9147,
0.9074,0.9000,0.8926,0.8853,0.8779,0.8705,0.8632,
0.8558,0.8474,0.8389,0.8263,0.8084,0.7895,0.7684,
0.7474,0.7211,0.6947,0.6579,0.6211,0.5842,0.5474,
0.5105,0.4737,0.4337,0.3895,0.3463,0.3053,0.2674,
0.2326,0.2000,0.1684,0.1389,0.1105,0.0758,0.0368
0.0000
interp_type = 2
extrap_type = 2
t_extrap_type = 2
lowest_lev_from_sfc = .false.
use_levels_below_ground = .true.
use_surface = .true.
use_tavg_for_tsk = .false.
lagrange_order = 1
force_sfc_in_vinterp = 1
zap_close_levels = 500
sfcp_to_sfcp = .false.
adjust_heights = .false.
smooth_cg_topo = .false.
dx = 27000, 9000, 3000,
dy = 27000, 9000, 3000,
grid_id = 1, 2, 3,
parent_id = 1, 1, 2,
parent_grid_ratio = 1, 3, 3,
i_parent_start = 1, 64, 158,
j_parent_start = 1, 45, 70,
parent_time_step_ratio = 1, 3, 3,
feedback = 0,
smooth_option = 0,
/
&physics
mp_physics = 8, 8, 8, 8,
progn = 0,
ra_lw_physics = 4, 4, 4, 4,
ra_sw_physics = 4, 4, 4, 4,
radt = 5, 5, 5, 5,
swrad_scat = 1
shadlen = 25000.
slope_rad = 0
topo_shading = 0
sf_sfclay_physics = 5, 5, 5, 5,
bl_pbl_physics = 5, 5, 5, 5,
bldt = 0, 0, 0, 0,
iz0tlnd = 1,
sf_urban_physics = 0, 0, 1, 1,
pxlsm_smois_init = 0, 0, 0, 0,
fractional_seaice = 0
topo_wind = 0, 0, 0, 0,
grav_settling = 0
sf_surface_physics = 2, 2, 2, 2,
num_soil_layers = 4,
surface_input_source = 1,
num_land_cat = 28,
num_soil_cat = 16,
usemonalb = .true.
rdmaxalb = .true.
rdlai2d = .false.
sst_update = 1
isfflx = 1,
ifsnow = 0,
icloud = 1,
cu_physics = 3, 3, 0, 0,
cugd_avedx = 1,
cu_diag = 1, 1, 0, 0,
prec_acc_dt = 180.,
cudt = 0, 0, 0, 0,
cu_rad_feedback = .true., .true., .false.,.false.,
ishallow = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
/
&dfi_control
dfi_opt = 0
/
&fdda
/
&dynamics
rk_ord = 3
c_s = 0.25
w_damping = 1
diff_opt = 2, 2, 2, 2
km_opt = 4, 4, 4, 4
damp_opt = 3,
zdamp = 5000., 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2, 0.2,
gwd_opt = 0,
khdif = 0, 0, 0, 0,
kvdif = 0, 0, 0, 0,
smdiv = 0.1
emdiv = 0.01
epssm = 0.1
do_avgflx_em = 1, 1, 1, 1,
do_avgflx_cugd = 1,
mix_full_fields = .false.
mix_isotropic = 0, 0, 0, 0,
mix_upper_bound = 0.1, 0.1, 0.1, 0.1,
diff_6th_opt = 2, 2, 2, 2,
diff_6th_factor = 0.12, 0.12, 0.12, 0.12,
base_temp = 290.
non_hydrostatic = .true., .true., .true., .true.,
moist_adv_opt = 2, 2, 2, 2,
chem_adv_opt = 2, 2, 2, 2,
scalar_adv_opt = 2, 2, 2, 2,
tke_adv_opt = 2, 2, 2, 2,
time_step_sound = 4, 4, 4, 4,
h_mom_adv_order = 5, 5, 5, 5,
v_mom_adv_order = 3, 3, 3, 3,
h_sca_adv_order = 5, 5, 5, 5,
v_sca_adv_order = 3, 3, 3, 3,
tke_drag_coefficient = 0., 0., 0., 0.,
tke_heat_flux = 0., 0., 0., 0.,
top_lid = .false., .false., .false., .false.,
/
&bdy_control
spec_bdy_width = 7,
spec_zone = 1,
relax_zone = 6,
spec_exp = 0.0,
specified = .true.,
nested = .false., .true., .true., .true.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&diags
output_velvar = 1, 1, 1
/
&chem
chem_opt = 14, 14, 14, 14,
chem_in_opt = 1, 1, 1, 1,
have_bcs_chem = .true., .false., .false., .false.,
have_bcs_tracer = .false., .false., .false., .false.,
!emi_inname = "wrfchemi_d<domain>_<date>",
!input_chem_inname = "wrfchemi_d<domain>_<date>",
gaschem_onoff = 0, 0, 0, 0,
aerchem_onoff = 0, 0, 0, 0,
vertmix_onoff = 1, 1, 1, 1,
chem_conv_tr = 1, 1, 0, 0,
gas_drydep_opt = 0, 0, 0, 0,
aer_drydep_opt = 0, 0, 0, 0,
diagnostic_chem = 0, 0, 0, 0,
chemdt = 0.033, 0.033, 0.033, 0.033,
bioemdt = 30, 30, 30, 30,
emiss_inpt_opt = 1, 1, 1, 1,
emiss_opt = 12, 12, 12, 12,
kemit = 1,
io_style_emissions = 2,
aircraft_emiss_opt = 0,
bio_emiss_opt = 0, 0, 0, 0,
phot_opt = 0, 0, 0, 0,
photdt = 30, 30, 30, 30,
wetscav_onoff = 0, 0, 0, 0,
cldchem_onoff = 0, 0, 0, 0,
conv_tr_aqchem = 0, 0, 0, 0,
conv_tr_wetscav = 0, 0, 0, 0,
seas_opt = 0,
dust_opt = 0,
dmsemis_opt = 0,
biomass_burn_opt = 0, 0, 0, 0,
plumerisefire_frq = 30, 30, 30, 30,
gas_bc_opt = 16, 16, 16, 16,
gas_ic_opt = 16, 16, 16, 16,
aer_bc_opt = 1, 1, 1, 1,
aer_ic_opt = 1, 1, 1, 1,
aer_ra_feedback = 0, 0, 0, 0,
opt_pars_out = 0,
/