Process Analysis Diagnostics

From Geos-chem
Revision as of 00:37, 14 September 2013 by Barronh (Talk | contribs) (New page: <big><big><strong>Process Analysis Diagnostics</strong></big></big> 1. Overview - What is Process Analysis? 2. Details - No really, what is it? 3. Development Status - What features are a...)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Process Analysis Diagnostics

1. Overview - What is Process Analysis? 2. Details - No really, what is it? 3. Development Status - What features are available? 4. Installation and Application - How can I use it? 5. Example Outputs

Overview

Process-based Analysis examines the change in each species due to each process and reaction. Models predict atmospheric state, which in a time-series can be used to create net-change of each species. What this cannot tell us, is which processes led to that change. To supplement state (or concentration, GEOS-Chem has long archived emissions and employed advanced diagnostics to predict gross chemical production or loss. Process Analysis goes a step further archiving grid-cell budgets for each species, and decomposing gross production/loss into individual reaction contributions.

Details

The budget is built from process contributions and reaction throughput. Process contributions are referred to as Integrated Process Rates (IPR) and can be calculated because of operator splitting. Each process is distinctly called, and state before and after the call can be differenced to see the process contribution.

Reaction contributions are more complex to calculate than IPR. Reaction contributions are constructed from integrated reaction rates (IRR). IRR are calculated according to equation 1. The rate coefficient (<math>k</math>) is evaluated by calcrate.f and then SMVGEAR or KPP solve over short timesteps (<math>t_s</math>). To calculate the IRR, we must add code to SMVGEAR or KPP that finds successful steps and archives equation 1.

<math>IRR_{rxn} = \Delta t_s\cdot k_{rxn} \cdot \prod_{rct}{\([C_{rct}]\)}</math> <math>\Delta C_{i,irr} = \alpha_{i,rxn}IRR_{rxn}</math>

Development Status

Installation and Application