Scheduled Downtime
On Tuesday 24 October 2023 @ 5pm MT the forums will be in read only mode in preparation for the downtime. On Wednesday 25 October 2023 @ 5am MT, this website will be down for maintenance and expected to return online later in the morning.
Normal Operations
The forums are back online with normal operations. If you notice any issues or errors related to the forums, please reach out to help@ucar.edu

greenhouse gases chemistry running error

L

ljdiao@gmail_com

Member
I am trying to run WACCM with ghg chemistry. But I got the error message as (shr_sys_abort)  ERROR: chem_surfvals_readnl:: ERROR reading namelist How can I fixes this problem? BTW, for ghg chemistry, is there only the greenhouse gases in the chemistry? thanks,
 
E

ep2764@columbia_edu

Member
I don't know too much about running with ghg chem, but the error you are getting is a namelist error. This is one of the earliest steps in running the model, so debugging shouldn't take too long.What changes did you make to the namelist and/or to the source code? Don't know if I can help, but I can try! -Ethan
 
mmills

mmills

CSEG and Liaisons
Staff member
The easiest way to set up ghg chemistry is by using a specified chemistry (SC-WACCM) compset. We have one set up for present-day conditions (FWSC) in CESM1.0.5, 1.1.1, and 1.2. If you are using 1.2, you can create custom compsets for other conditions using an XML schema that is described when you type "create_newcase -list compsets" under $srcdir/scripts.The ghg chemistry is described in chem_mech.in and chem_mech.doc in models/atm/cam/src/chemistry/pp_waccm_ghg. It includes the following species:      Solution species
  1. CH4             
  2. N2O             
  3. CFC11 (CFCl3)
  4. CFC12 (CF2Cl2)
  5. H2O             
     Invariant species
  1. M               
  2. N2              
  3. O2              
     Column integrals
  1. O3
  2. O2
The following species are prescribed by the ozone_forcing_file and waccm_forcing_file in the namelist, which is monthly zonal means:
Code: 
        float CO2(time, lev, lat) ;

Code: 
                CO2:units = "mol/mol" ;

Code: 
                CO2:long_name = "CO2 concentration" ;

Code: 
                CO2:cell_methods = "time: mean" ;

Code: 
        float H(time, lev, lat) ;

Code: 
                H:units = "mol/mol" ;

Code: 
                H:long_name = "H concentration" ;

Code: 
                H:cell_methods = "time: mean" ;

Code: 
        float NO(time, lev, lat) ;

Code: 
                NO:units = "mol/mol" ;

Code: 
                NO:long_name = "NO concentration" ;

Code: 
                NO:cell_methods = "time: mean" ;

Code: 
        float O(time, lev, lat) ;

Code: 
                O:units = "mol/mol" ;

Code: 
                O:long_name = "O concentration" ;

Code: 
                O:cell_methods = "time: mean" ;

Code: 
        float O2(time, lev, lat) ;

Code: 
                O2:units = "mol/mol" ;

Code: 
                O2:long_name = "O2 concentration" ;

Code: 
                O2:cell_methods = "time: mean" ;

Code: 
        float O3(time, lev, lat) ;

Code: 
                O3:units = "mol/mol" ;

Code: 
                O3:long_name = "O3 concentration" ;

Code: 
                O3:cell_methods = "time: mean" ;

Code: 
        float QRS_TOT(time, lev, lat) ;

Code: 
                QRS_TOT:units = "K/s" ;

Code: 
                QRS_TOT:long_name = "Merged SW heating: QRS+QCP+QRS_EUV+QRS_CO2NIR+QRS_AUR+QTHERMAL" ;

Code: 
                QRS_TOT:cell_methods = "time: mean" ;
  
 
You must log in or register to reply here.
Top