registry.chem variables are not saving to wrfout in v3.7.1

Posts related to the installation of WRF Chemistry

registry.chem variables are not saving to wrfout in v3.7.1

Postby mike_wrf » Wed Jan 27, 2016 6:01 am

Hi All

I am adding variables to wrfout using iofield_list.txt.
This used to work in v3.4.1 and save the wet dep variables. But does not work in v3.7.1.

So name list option is:
iofields_filename = "iofield_list.txt","iofield_list.txt",

iofield_list.txt. look like:
+:h:0:wd_so4_sc,wd_no3_sc,wd_so4_cu,wd_no3_cu

registry.chem looks like:
state real wd_so4_sc ij misc 1 - rdu "wd_so4_sc" "SO4 surface wet deposition, accumulated (Sc)" "mmol/m2"
state real wd_no3_sc ij misc 1 - rdu "wd_no3_sc" "NO3 surface wet deposition, accumulated (Sc)" "mmol/m2"

This used to work in v3.4.1 and save the wet dep variables. But does not work in v3.7.1.

rsl.error.0000 says:
med_initialdata_input: calling input_input
W A R N I N G : Unable to modify mask for wd_so4_sc. Variable not found. File: iofield_list.txt at line 1

So, I don't get why this should cause an error because it is not in the input, but in the output? And It worked with the "0" option in iofields_list.txt previously (where 0 = input and output).

Any ideas as to why it works in one version (v3.4.1) and not v3.7.1?
mike_wrf
 
Posts: 7
Joined: Wed Aug 25, 2010 4:11 am

Re: registry.chem variables are not saving to wrfout in v3.7

Postby mike_wrf » Thu Mar 10, 2016 5:04 am

The answer to my problem seems to be with the registry.chem file in v3.7.1.
Thanks to Ravan Ahmadov from NOAA for feedback through the wrf chem help email address!

The wd* variables were included in the chem array of 202. This needs to be fixed.
i.e. line starting with;
package mozart_mosaic_4bin_aq_kpp chem_opt==202
in registry.chem.
I removed the wd_* variables from this line and recompiled wrfchem.
I am not using chem_opt=202.
After the recompile, the wd_ variables are saving to my wrfout again :D
mike_wrf
 
Posts: 7
Joined: Wed Aug 25, 2010 4:11 am

Re: No data into wrfchemi file.

Postby wrfchemuser » Thu Apr 28, 2016 9:14 pm

Hi, I'm having some difficulty with something different. I'm currently running wrfchem and I'm trying to get my wrfchemi files from the NEI11 database. I've been able to successfully complete both getting the intermediate binary files and the wrfchemi netcdf file but I've noticed that there is no data and only global attributes are within the file when I perform a ncview/ncdump command. My namelist.input options are described below:

&time_control
run_days = 0,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2012, 2000, 2000,
start_month = 07, 01, 01,
start_day = 02, 24, 24,
start_hour = 00, 12, 12,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2012, 2000, 2000,
end_month = 07, 01, 01,
end_day = 02, 25, 25,
end_hour = 11, 12, 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 = 5000,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
debug_level = 0
!io_form_auxinput2 = 2
io_form_auxinput5 = 2
!auxinput5_inname = "wrfchemi_d<domain>_<date>"
auxinput5_interval = 3600
auxinput1_inname = "met_em.d<domain>.<date>"
/

&domains
time_step = 15,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 75, 112, 94,
e_sn = 70, 97, 91,
e_vert = 30, 30, 30,
p_top_requested = 5000,
num_metgrid_levels = 40,
num_metgrid_soil_levels = 4,
dx = 2500, 10000, 3333.33,
dy = 2500, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 31, 30,
j_parent_start = 1, 17, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 0
/

&physics
mp_physics = 3, 3, 3,
ra_lw_physics = 1, 1, 1,
ra_sw_physics = 1, 1, 1,
radt = 30, 30, 30,
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, 1, 0,
cudt = 5, 5, 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0, 0, 0,
cu_diag = 1,
cu_rad_feedback = .true,
cugd_avedx = 3,
/

&fdda
/

&dynamics
!w_damping = 0,
diff_opt = 1, 1, 1,
km_opt = 4, 4, 4,
diff_6th_opt = 0, 0, 0,
!diff_6th_factor = 0.12, 0.12, 0.12,
base_temp = 290.
damp_opt = 0,
zdamp = 5000., 5000., 5000.,
dampcoef = 0.2, 0.2, 0.2
khdif = 0, 0, 0,
!kvdif = 0, 0, 0,
non_hydrostatic = .true., .true., .true.,
!moist_adv_opt = 1, 1, 1,
!scalar_adv_opt = 1, 1, 1,
chem_adv_opt = 2,
moist_adv_opt = 2,
scalar_adv_opt = 2,
tke_adv_opt = 2,
/

&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
chem_opt = 1, 1,
chem_in_opt = 0, 1,
chemdt = 0., 3.,
bio_emiss_opt = 1, 3,
bioemdt = 2., 0.4,
phot_opt = 1, 2,
photdt = 30., 30.,
emiss_opt = 3, 13,
emiss_inpt_opt = 1, 102,
kemit = 19,
io_style_emissions = 2,
!frames_per_emissfile = 12,
biomass_burn_opt = 0, 0,
gas_bc_opt = 1, 1,
gas_ic_opt = 1, 1,
aer_bc_opt = 1, 1,
aer_ic_opt = 1, 1,
gaschem_onoff = 1, 1,
aerchem_onoff = 1, 1,
wetscav_onoff = -10, -10,
cldchem_onoff = 0, 0,
vertmix_onoff = 1, 1,
chem_conv_tr = 1, 0,
aer_ra_feedback = 0, 1,
have_bcs_chem = .true., .false.,
have_bcs_upper = .false., .false.,
ne_area = 140,
/

Is there something I need to change so that when I run convert_emiss.exe on the NEI11 emissions, it will give me the appropriate amount of data within the file?
wrfchemuser
 
Posts: 2
Joined: Thu Apr 14, 2016 5:08 pm


Return to Installation

Who is online

Users browsing this forum: No registered users and 4 guests