diff options
-rw-r--r-- | src/AHFinder_fun.F | 14 |
1 files changed, 1 insertions, 13 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F index d09455f..a6a426e 100644 --- a/src/AHFinder_fun.F +++ b/src/AHFinder_fun.F @@ -100,13 +100,12 @@ ! Non-axisymmetric terms. Notice how the sum over M is first. ! This will allow me to use the recursion relations to avoid ! having to start from scratch every time. Also, I sum over -! all l's even if I don't want some terms. This is because +! all l even if I do not want some terms. This is because ! in order to use the recursion relations I need all polynomials. if (nonaxi) then do m=1,lmax aux2 = dble(m)*phi - sina = sin(aux2) cosa = cos(aux2) do l=m,lmax @@ -121,17 +120,6 @@ ahfgrid(i,j,k) = rp - aux1 -! Check for NaNs. - - if (ahfgrid(i,j,k).ge.zero) then - else if (ahfgrid(i,j,k).lt.zero) then - else - if (myproc.eq.0) then - write(*,*) 'NaN in ahfgrid at point',i,j,k - end if - STOP - end if - end do end do end do |