fix fac in swegn96

keurfonluu · Feb 7, 2021 · 949e787 · 949e787
1 parent 10dcb4c
commit 949e787
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/disba/_cps/_swegn96.py b/disba/_cps/_swegn96.py
@@ -505,9 +505,9 @@ def shfunc(omega, wvno, d, a, b, rho):
     for i in range(1, mmax):
         if b[i] > 0.01:
             ext += exl[i - 1]
-            fac = numpy.exp(ext) if ext < 80.0 else 0.0
-            uu[i] /= fac
-            tt[i] /= fac
+            fac = 1.0 / numpy.exp(ext) if ext < 80.0 else 0.0
+            uu[i] *= fac
+            tt[i] *= fac
         else:
             uu[i] = 0.0
             tt[i] = 0.0