Removed divisions by zero when calculating coefficients for the spherical harmonic recursion relations, that were causing problems on Decs.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@240 89daf98e-ef62-4674-b946-b8ff9de2216c

1 files changed, 7 insertions, 3 deletions

diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F
index 8324a3e..f34d656 100644
--- a/src/AHFinder_fun.F
+++ b/src/AHFinder_fun.F
@@ -91,13 +91,17 @@
       if (nonaxi) then
          do l=1,lmax
             dsqrt_factor(l,1) = dsqrt(2*dble(2*l+1)/dble(l*(l+1)))
-            lm_factor1(l,1)=dble(2*l-1)/dble(l-1)
-            lm_factor2(l,1)=dble(l)/dble(l-1)
             do m=2,l
                dsqrt_factor(l,m)=dsqrt_factor(l,m-1)/dsqrt(dble((l-m+1)*(l+m)))
+           end do
+         end do
+         do l=2,lmax
+            lm_factor1(l,1)=dble(2*l-1)/dble(l-1)
+            lm_factor2(l,1)=dble(l)/dble(l-1)
+            do m=2,l-1
                lm_factor1(l,m)=dble(2*l-1)/dble(l-m)
                lm_factor2(l,m)=dble(l+m-1)/dble(l-m)
-            end do
+           end do
          end do
       end if