Difference between revisions of "Additional diagnostics"
(vertical wind derivation from standard diagnostics) |
|||
Line 26: | Line 26: | ||
Air_Density [Kg/m3] * Box_Height [m] / Airmass[Kg] = 1 / surface [m2] | Air_Density [Kg/m3] * Box_Height [m] / Airmass[Kg] = 1 / surface [m2] | ||
− | so you can get the box surface from the grid description (see [http:// | + | so you can get the box surface from the grid description (see [http://acmg.seas.harvard.edu/gamap/doc/by_alphabet/gamap_c.html#CTM_BOXSIZE%20%20(FUNCTION) ctm_boxsize]) and replace (1) with: |
dP/dt = Air_Mass_Flux [Kg/s] * 9.8 [m/s2] / surface [m2] (2) | dP/dt = Air_Mass_Flux [Kg/s] * 9.8 [m/s2] / surface [m2] (2) |
Latest revision as of 13:48, 7 July 2009
Additional diagnostics can be obtained by combining the standard GEOS-Chem outputs. This page presents some of them.
Vertical wind
Since dP/dt is a conventional unit in atmospheric chemistry, we can proceed as follows to convert a vertical flux into vertical velocity.
The parameters used are ND45 (mixing ratio), ND26 (upward flux), and ND68 (box height, air mass, air density).
The air density must be converted from [molec/m3] to [Kg/m3]:
Air_Density [Kg/m3] = Air_Density [molec/m3] * 28.9e3 / Avogadro_Number
To calculate the vertical velocity you must choose one tracer that you included in the ND26 diagnostic. For example, if we use CO (molecular weight 28):
Air_Mass_Flux [Kg/s] = ND26_Up_Flux(CO) [Kg/s] / ND45_Mixing_Ratio(CO) * 28.9 / 28 Vertical_velocity [m/s] = Air_Mass_Flux [Kg/s] * Box_Height [m]/Airmass[Kg]
or in the preferred dP/dt format:
dP/dt = Air_Mass_Flux [Kg/s] * Box_Height [m]/Airmass[Kg] * Air_Density [Kg/m3] * 9.8 [m/s2] (1)
Note that
Air_Density [Kg/m3] * Box_Height [m] / Airmass[Kg] = 1 / surface [m2]
so you can get the box surface from the grid description (see ctm_boxsize) and replace (1) with:
dP/dt = Air_Mass_Flux [Kg/s] * 9.8 [m/s2] / surface [m2] (2)
By using (2), you do not need ND68 anymore.
--phs 17:19, 15 January 2009 (EST)