noresm:aerosolbudgetsnoresm2

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision 		Previous revision
Next revision Both sides next revision
noresm:aerosolbudgetsnoresm2 [2015-03-25 15:49:06]
michaels 		noresm:aerosolbudgetsnoresm2 [2015-12-16 20:52:47]
michaels
Line 1: Line 1:
-===Configuring a run with aerosol budgets (NorESM2) ===+===Configuring a run with more aerosol diagnostics in (NorESM2) ===
  
 In both the default CAM5-aerosol packages (MAM3,MAM7) and the Oslo-aerosol packages, the budget terms can be taken out using a variable in the namelist :  In both the default CAM5-aerosol packages (MAM3,MAM7) and the Oslo-aerosol packages, the budget terms can be taken out using a variable in the namelist : 
Line 8: Line 8:
 / /
 </file> </file>
 +
 +Two more diagnostics are useful:
 +  * Enable estimates multiple calls to radiation which are necessary for effective radiative forcing estimates
 +  * Enable diagnostics for AEROCOM
 +
 +
 +To enable this, take the file cam/src/physics/cam_oslo$ vim preprocessorDefinitions.h and copy it to your SourceMods/src.cam folder
 +
 +Change both preprocessor definitions to true
 +
 +<file>
 +#define AEROCOM
 +#define AEROFFL
 +</file>
 +
 +The AEROCOM-token turns on diagnostics needed for AEROCOM
 +The AEROFFL-token tells the model to do additional radiation-diagnostics for aerosol indirect effect
  
 ===Fields produced in monthly average files when running with budgets activated=== ===Fields produced in monthly average files when running with budgets activated===
Line 24: Line 41:
 | DF_{Tracer}             | dry deposition flux (gas tracers)                      | output with history_aerosol with CAM-Oslo only| | DF_{Tracer}             | dry deposition flux (gas tracers)                      | output with history_aerosol with CAM-Oslo only|
 | WD_A_{Tracer}           | wet deposititon flux (gas tracers)                     | output with history_aerosol with CAM-Oslo only| | WD_A_{Tracer}           | wet deposititon flux (gas tracers)                     | output with history_aerosol with CAM-Oslo only|
-| Tracer_{CLXF}           | 3D-emissions ("external forcing")                      | output with history_aerosol with CAM-Oslo only|+| {Tracer}_CLXF           | 3D-emissions ("external forcing")                      | output with history_aerosol with CAM-Oslo only| 
 +| {Tracer}_clcoagTend     | loss of tracer due to coagulation with cloud droplets  | output with history_aerosol with CAM-Oslo only|
  
 Note: Since 3D-emissions and and gas washout rates are included in the term GS_{Tracer} in the mozart chemistry solver, the individual terms can be found like this (example for SO2):  Note: Since 3D-emissions and and gas washout rates are included in the term GS_{Tracer} in the mozart chemistry solver, the individual terms can be found like this (example for SO2): 
Line 30: Line 48:
 ncap2 -O -s GS_ONLY_SO2=GS_SO2-WD_A_SO2-SO2_CLXF infile.nc outfile.nc ncap2 -O -s GS_ONLY_SO2=GS_SO2-WD_A_SO2-SO2_CLXF infile.nc outfile.nc
 </file> </file>
 +
 +More info on SO2 budgets (see /models/atm/cam/tools/diagnostics/ncl/ModIvsModII/ for scripts with info on all tracers): 
 +
 +GS_SO2 contains the SO2 budget terms for all that goes on in the chemistry-routine, which is \\
 +1) Gas phase chemistry, 2) Wet deposition, and 3) 3D-emissions.\\
 +Gas phase chemistry is both production from DMS (GS_DMS) and loss through OH (GL_OH) \\
 +For calculations of net loss, e.g. used to calculate SO2 life-times, we're interested in the \\
 +loss through OH from the chemistry-term (GL_OH).\\
 +GS_SO2 = GL_OH + SO2_CLXF - WD_A_SO2 - GS_DMS*64/62 \\
 +or \\
 +GL_OH = GS_SO2 - SO2_CLXF + WD_A_SO2 + GS_DMS*64/62 \\
 +
 +Estimating chemical loss w.r.t. S (instead of SO2 or DMS), for comparison with CAM4-Oslo numbers:\\
 +net chemial loss gas phase = (GS_SO2/1.998 - SO2_CLXF + WD_A_SO2)/1.998  + GS_DMS/1.938 \\
 +net chemical loss = net chemial loss gas phase + AQ_SO2/1.998 \\
 +
 +Finally, total net loss (used to calculate life-time = -load/(net loss), where load = cb_SO2/1.998):\\
 +net loss = \\
 +- WD_A_SO2/1.998    ;wet deposition in kg/m2/sec (positive in output file) \\
 +- DF_SO2/1.998      ;dry deposition in kg/m2/sec (positive in output file) \\
 ++ AQ_SO2/1.998      ;wet phase production of SO4 in kg/m2/ses (negative in output file) \\
 ++ (GS_SO2 - SO2_CLXF + WD_A_SO2)/1.998 + GS_DMS/1.938  ; net chemical loss gas phase \\
  
 ===Looking at the aerosol budgets (CAM-Oslo only) === ===Looking at the aerosol budgets (CAM-Oslo only) ===
Line 37: Line 77:
   * Run the script budgets.sh to create a pdf-file (output.pdf)   * Run the script budgets.sh to create a pdf-file (output.pdf)
  
-===Preparing output for AEROCOM ===+===Making ncl plots of often used aerosol and cloud fields, including ERFs, for two model versions (CAM-Oslo only) ===
  
-To prepare output for aerocom, use the script located at models/atm/cam/tools/aerocom/ in the svn repository. +  * Make a local copy (on Linux) of the directory models/atm/cam/tools/diagnostics/ncl/ModIvsModII 
 +  * Assuming that you have produced output data from 4 simulations: two different model versions, each with PD and PI emissions, and all run with #define AEROCOM & AEROFFL:    
 +  * In ModIvsModII.csh (note: read the header info): 
 +  * - edit model info for the first model (shown to the left in the plots): modelI = CAM4-Oslo or modelI = CAM5-Oslo ? 
 +  * - provide paths and partial file names of the model data (PD and PI) for Model I (CAM4-Oslo or CAM5-Oslo) and Model II (must be CAM5-Oslo) 
 +  * - choose desired plot format (plotf=ps, eps, pdf or png) 
 +  * Run the script: ./ModIvsModII.csh 
 +  * Furthermore, to display the plots in an organized form by use of a web browser (only possible if the chosen plot format is png):     
 +  * - download htm template files from ftp://ftp.met.no/projects/noresmatm/upload/NorESM2Diagnostics/ModIvsModII/htm-templates/ 
 +  * - edit general model info (only) in ModIvsModII.htm, and manually cut and paste the mass budget numbers from the script output into this file  
 +  * - copy all png (plots) and htm files to the desired output (common) directory 
 +  * - open ModIvsModII.htm in your browser: hyper-links to all other htm files, including plots, are found here 
 +  * Example: ftp://ftp.met.no/projects/noresmatm/upload/NorESM2Diagnostics/ModIvsModII/revision610inclSOA-Nudged_1984-12to1985-11_vs_CAM4-Oslo/ModIvsModII.htm 
 + 
 +===Preparing output for AEROCOM analysis=== 
 + 
 +To prepare output for the aerocom tools, use the script located at **models/atm/cam/tools/aerocom/** in the svn repository. 
  
 The script requires <ModelName>_<ExperimentName> and <Period> as input. The script requires <ModelName>_<ExperimentName> and <Period> as input.
 +
 +<ModelName>_<ExperimentName> is the identifier under which the plots appear on the AeroCom webinterface.
 +
  
 The script creates files named like The script creates files named like
Line 53: Line 112:
 <Period> ⇒ eg “2008”, “2010”, "9999" \\ <Period> ⇒ eg “2008”, “2010”, "9999" \\
 <Frequency> ⇒ “timeinvariant”,”hourly”, “daily”, “monthly”, “sat1000”, “sat1330”, “sat2200”, “sat0130” \\ <Frequency> ⇒ “timeinvariant”,”hourly”, “daily”, “monthly”, “sat1000”, “sat1330”, “sat2200”, “sat0130” \\
 +
 +Note that VerticalCoordinateType is dependent on the variable!! It is not a question about "vertical coordinate type used in model simulations"!
  
 The script copies files on norstore into The script copies files on norstore into
-/projects/NS2345K/CAM-Oslo/DO_AEROCOM/<ModelName>_<ExperimentName>/renamed/+**/projects/NS2345K/CAM-Oslo/DO_AEROCOM/<ModelName>_<ExperimentName>/renamed/** 
  
  • noresm/aerosolbudgetsnoresm2.txt
  • Last modified: 2022-05-31 09:29:32
  • (external edit)