Skip to content

Instantly share code, notes, and snippets.

@spestana

spestana/dah.py

Created January 21, 2022 22:26
Show Gist options
  • Save spestana/4fced6c34c4ae821c67ccfdae8283bf5 to your computer and use it in GitHub Desktop.
Save spestana/4fced6c34c4ae821c67ccfdae8283bf5 to your computer and use it in GitHub Desktop.
Download ZIP
Raw
dah.py
def diurnal_anisotropic_heat(aspect, slope, a_max):
'''
Make a Diurnal Anisotropic Heating (DAH) index (Böhner & Antonic, 2009)
a_max = 202.5 (SSW) from Cristea et al., 2017 in Tuolumne area
'''
DAH = np.cos(np.radians(a_max - aspect)) * np.arctan(np.radians(slope))
return DAH
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment