diff options
Diffstat (limited to 'src/AHFinder_fun.F')
-rw-r--r-- | src/AHFinder_fun.F | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F index e565782..0af446c 100644 --- a/src/AHFinder_fun.F +++ b/src/AHFinder_fun.F @@ -99,7 +99,7 @@ do l=1,lmax - dsqrt_factor2(l) = c0(l)*dsqrt(dble(2*l+1)) + dsqrt_factor2(l) = c0(l)*sqrt(dble(2*l+1)) lm_factor1(l,0) = dble(2*l-1)/dble(l) lm_factor2(l,0) = dble(l-1)/dble(l) @@ -113,11 +113,11 @@ do l=1,lmax - dsqrt_factor(l,1) = dsqrt(2.0D0*dble(2*l+1)/dble(l*(l+1))) + dsqrt_factor(l,1) = sqrt(2.0D0*dble(2*l+1)/dble(l*(l+1))) do m=2,l dsqrt_factor(l,m) = dsqrt_factor(l,m-1) - . / dsqrt(dble((l-m+1)*(l+m))) + . / sqrt(dble((l-m+1)*(l+m))) end do end do @@ -237,7 +237,7 @@ ! Find rp. - rp = dsqrt(xp**2 + yp**2 + zp**2) + rp = sqrt(xp**2 + yp**2 + zp**2) ! Monopole term. @@ -281,7 +281,7 @@ ! if (nonaxi) then - aux2 = dsqrt((one-cost)*(one+cost)) + aux2 = sqrt((one-cost)*(one+cost)) ! This general loop can only be used until l=lmax-2 since ! the recursion relations requires 2 starting values. |