[Python-modules-commits] [python-pysolar] 15/42: Add aerosol correction factor equations

[Wolfgang Borgert]

This is an automated email from the git hooks/post-receive script.

debacle pushed a commit to tag 0.6
in repository python-pysolar.

commit b421aba8a6ceddba66a082c8efe5f20d7a321402
Author: Brandon Stafford <brandon at pingswept.org>
Date:   Fri Mar 28 21:52:40 2014 -0400

    Add aerosol correction factor equations
---
 Pysolar/rest.py | 12 ++++++++++++
 1 file changed, 12 insertions(+)

diff --git a/Pysolar/rest.py b/Pysolar/rest.py
index 54f095d..32ed050 100644
--- a/Pysolar/rest.py
+++ b/Pysolar/rest.py
@@ -37,6 +37,18 @@ def GetAerosolForwardScatteranceFactor(altitude_deg):
 	Z = 90 - altitude_deg
 	return 1 - math.e ** (-0.6931 - 1.8326 * math.cos(math.radians(Z)))
 
+def GetAerosolScatteringCorrectionFactor(band, ma, tau_a):
+	if band == 1:
+		g0 = (3.715 + 0.368 * ma + 0.036294 * ma ** 2)/(1 + 0.0009391 * ma ** 2)
+		g1 = (-0.164 - 0.72567 * ma + 0.20701 * ma ** 2)/(1 + 0.001901 * ma ** 2)
+		g2 = (-0.052288 + 0.31902 * ma + 0.17871 * ma ** 2)/(1 + 0.0069592 * ma ** 2)
+		return (g0 + g1 * tau_a)/(1 + g2 * tau_a)
+	else:
+		h0 = (3.4352 +  0.65267 * ma + 0.00034328 * ma ** 2)/(1 + 0.034388 * ma ** 1.5)
+		h1 = (1.231 - 1.63853 * ma + 0.20667 * ma ** 2)/(1 + 0.1451 * ma ** 1.5)
+		h2 = (0.8889 - 0.55063 * ma + 0.50152 * ma ** 2)/(1 + 0.14865 * ma ** 1.5)
+		return (h0 + h1 * tau_a)/(1 + h2 * tau_a)
+
 def GetRayleighExtinctionForwardScatteringFraction(air_mass, layer):
 	if layer == 1:
 		return 0.5 * (0.89013 - 0.049558 * air_mass + 0.000045721 * air_mass ** 2)

-- 
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/python-modules/packages/python-pysolar.git