Skip to content

Commit

Permalink
ENH: Add documentation for galvez_davison_index
Browse files Browse the repository at this point in the history
  • Loading branch information
alakocy-pec authored and kgoebber committed Sep 21, 2023
1 parent d95e789 commit c96683d
Show file tree
Hide file tree
Showing 4 changed files with 29 additions and 15 deletions.
1 change: 1 addition & 0 deletions docs/_templates/overrides/metpy.calc.rst
Original file line number Diff line number Diff line change
Expand Up @@ -76,6 +76,7 @@ Soundings
critical_angle
cross_totals
el
galvez_davison_index
k_index
lcl
lfc
Expand Down
4 changes: 4 additions & 0 deletions docs/api/references.rst
Original file line number Diff line number Diff line change
Expand Up @@ -76,6 +76,10 @@ References
Services and Supporting Research, 2003.
`FCM-R19-2003 <../_static/FCM-R19-2003-WindchillReport.pdf>`_, 75 pp.
.. [Galvez2015] Galvez, J. M. and Michel Davison, 2015. “The Gálvez-Davison Index for Tropical
Convection.” `GDI_Manuscript_V20161021
<https://www.wpc.ncep.noaa.gov/international/gdi/GDI_Manuscript_V20161021.pdf>`_.
.. [Galway1956] Galway, J. G., 1956: The Lifted Index as a Predictor of Latent Instability.
*American Meteorology Society*,
doi:`10.1175/1520-0477-37.10.528
Expand Down
1 change: 1 addition & 0 deletions examples/calculations/Sounding_Calculations.py
Original file line number Diff line number Diff line change
Expand Up @@ -96,6 +96,7 @@ def effective_layer(p, t, td, h, height_layer=False):
# Compute common sounding index parameters
ctotals = mpcalc.cross_totals(p, T, Td)
kindex = mpcalc.k_index(p, T, Td)
gdi = mpcalc.galvez_davison_index(p, T, Td)
showalter = mpcalc.showalter_index(p, T, Td)
total_totals = mpcalc.total_totals_index(p, T, Td)
vert_totals = mpcalc.vertical_totals(p, T)
Expand Down
38 changes: 23 additions & 15 deletions src/metpy/calc/thermo.py
Original file line number Diff line number Diff line change
Expand Up @@ -4339,7 +4339,11 @@ def galvez_davison_index(pressure, temperature, dewpoint, vertical_dim=0):
"""
Calculate GDI from the pressure temperature and dewpoint.
GDI formula derived from [Galvez2015](https://www.wpc.ncep.noaa.gov/international/gdi/GDI_Manuscript_V20150910.pdf)_:
Calculation of the GDI relies on temperatures and mixing ratios at 950,
850, 700, and 500 hPa. These four levels define three layers: A) Boundary,
B) Trade Wind Inversion (TWI), C) Mid-Troposphere.
GDI formula derived from [Galvez2015]_:
.. math:: GDI = CBI + MWI + II + TC
Expand All @@ -4350,20 +4354,24 @@ def galvez_davison_index(pressure, temperature, dewpoint, vertical_dim=0):
* :math:`II` is the Inversion Index
* :math:`TC` is the Terrain Correction [optional]
Calculation of the GDI relies on temperatures and mixing ratios at 950,
850, 700, and 500 hPa. These four levels define three layers: A) Boundary,
B) Trade Wind Inversion (TWI), C) Mid-Troposphere.
----------------------------------------------------------------------------------
GDI Value | Expected Convective Regime
----------------------------------------------------------------------------------
>=45 | Scattered to widespread thunderstorms likely.
35 to 45 | Scattered thunderstorms and/or scattered to widespread rain showers.
25 to 35 | Isolated to scattered thunderstorms and/or scattered showers.
15 to 25 | Isolated thunderstorms and/or isolated to scattered showers.
5 to 10 | Isolated to scattered showers.
<5 | Strong TWI likely, light rain possible.
----------------------------------------------------------------------------------
.. list-table:: GDI Values & Corresponding Convective Regimes
:widths: 15 75
:header-rows: 1
* - GDI Value
- Expected Convective Regime
* - >=45
- Scattered to widespread thunderstorms likely.
* - 35 to 45
- Scattered thunderstorms and/or scattered to widespread rain showers.
* - 25 to 35
- Isolated to scattered thunderstorms and/or scattered showers.
* - 15 to 25
- Isolated thunderstorms and/or isolated to scattered showers.
* - 5 to 10
- Isolated to scattered showers.
* - <5
- Strong TWI likely, light rain possible.
Parameters
----------
Expand Down

0 comments on commit c96683d

Please sign in to comment.