modify the chemistry scheme

[ljdiao@gmail_com]

I am running waccm_mozart but with only O2_1S species and related reactions deleted. When running, I got the error message as "ERROR: GETFIL: FAILED to get WACCM_FORCING_FILE"It doesn't make sense to me. Why is needs waccm_forcing_file? How to fix the problem then? thanks,

 

[mmills]

Please tell us which CESM code base you are using, which compset and resolution you are using, and what steps you took to modify the default setup.

 

[santos]

It looks like there is WACCM chemistry code (mo_waccm_hrates.F90) that will only calculate heating rates if certain species are present, and O2_1S is one of those. If these heating rates cannot be calculated, then WACCM will attempt to read forcing data from a file, which is normally only present for "waccm_ghg" (specified chemistry).I'm not sure about the appropriate way to deal with this; maybe a scientist can comment.

 

[mmills]

If the research goal is to study the atmosphere in the absence of O2_1S (including heating effects), you might try creating a waccm_forcing_file in which values of O2_1S are zero.For a model of such a file, look at inputdata/atm/waccm/ub/ghg_forcing_2000_c110321.nc, which is used for specified chemistry (waccm_ghg) runs. It does not include O2_1S, but you could create one that does.Mike Mills
WACCM Liaison
NCAR Foothills Lab

 

[ljdiao@gmail_com]

Hi,  I am using cesm 1.2.0. I built waccm compset B_2000_WACCM_CN with resolution f19_g16.  And then I used./xmlchange -file env_build.xml -id CAM_CONFIG_OPTS -val "-phys cam4 -chem waccm_mozart -usr_mech_infile /path/chem_mech.in"command to modify the chemistry, where the chem_mech.in file is modified from default waccm_mozart chem_mech.in file with only O2_1S related things deleted. thanks,

 

[santos]

Hmm. That seems like the easiest thing to do, but does that require you to specify all the constituents, or can you just force O2_1S to zero?