Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
noresm:aerosolbudgetsnoresm2 [2015-03-25 12:09:19] alfg |
noresm:aerosolbudgetsnoresm2 [2022-05-31 09:29:32] (current) |
||
|---|---|---|---|
| Line 1: | Line 1: | ||
| - | ===Configuring a run with aerosol | + | ===Configuring a run with more aerosol |
| 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: | ||
| / | / | ||
| </ | </ | ||
| + | |||
| + | 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/ | ||
| + | |||
| + | Change both preprocessor definitions to true | ||
| + | |||
| + | < | ||
| + | #define AEROCOM | ||
| + | #define AEROFFL | ||
| + | </ | ||
| + | |||
| + | 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} | | DF_{Tracer} | ||
| | WD_A_{Tracer} | | WD_A_{Tracer} | ||
| - | | Tracer_{CLXF} | 3D-emissions (" | + | | {Tracer}_CLXF | 3D-emissions (" |
| + | | {Tracer}_clcoagTend | ||
| 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 | ||
| </ | </ | ||
| + | |||
| + | More info on SO2 budgets (see / | ||
| + | |||
| + | GS_SO2 contains the SO2 budget terms for all that goes on in the chemistry-routine, | ||
| + | 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/ | ||
| + | or \\ | ||
| + | GL_OH = GS_SO2 - SO2_CLXF + WD_A_SO2 + GS_DMS*64/ | ||
| + | |||
| + | 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/ | ||
| + | net chemical loss = net chemial loss gas phase + AQ_SO2/ | ||
| + | |||
| + | Finally, total net loss (used to calculate life-time = -load/(net loss), where load = cb_SO2/ | ||
| + | net loss = \\ | ||
| + | - WD_A_SO2/ | ||
| + | - DF_SO2/ | ||
| + | + AQ_SO2/ | ||
| + | + (GS_SO2 - SO2_CLXF + WD_A_SO2)/ | ||
| ===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) === |
| + | |||
| + | * Make a local copy (on Linux) of the directory models/ | ||
| + | * Assuming that you have produced output data from 4 simulations: | ||
| + | * 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: ./ | ||
| + | * Furthermore, | ||
| + | * - download htm template files from ftp:// | ||
| + | * - edit general model info (only) in ModIvsModII.htm, | ||
| + | * - 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:// | ||
| + | |||
| + | ===Preparing output for AEROCOM | ||
| + | |||
| + | To prepare output so that it is processed automatically by the aerocom tools, use the script located at **models/ | ||
| + | [[http:// | ||
| + | |||
| + | The script requires < | ||
| + | |||
| + | < | ||
| + | |||
| + | < | ||
| + | of the form NorESM-CAM5_svn{RevisionNumber}_YYMMDD{initials}_Freetext. \\ | ||
| + | Example: " | ||
| + | Initials AG: Alf Grini, AK: Alf Kirkevåg, DO: Dirk Olivie... | ||
| + | |||
| + | The script creates files named like | ||
| + | |||
| + | “aerocom3_< | ||
| - | To prepare output for aerocom, use the script located at models/ | + | < |
| - | aerocom3.${MODELNAME}_${EXPERIMENTNAME}.SVNR{RevisionNumber}.monthly.${VariableName}.nc | + | < |
| + | < | ||
| + | < | ||
| + | < | ||
| + | < | ||
| - | Copy these files to the catalogue XXXX and they will be picked up by scripts which show them on the aerocom webpage | + | Note that VerticalCoordinateType is dependent |
| + | The script copies files on norstore into | ||
| + | **/ | ||