Removed ' in comments for preprocessor, removed ifdef MPIs, fixed

couple of compiler complaints


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

1 files changed, 6 insertions, 4 deletions

diff --git a/src/AHFinder_calcsigma.F b/src/AHFinder_calcsigma.F
index 1cfaf9b..a3b84ab 100644
--- a/src/AHFinder_calcsigma.F
+++ b/src/AHFinder_calcsigma.F
@@ -74,7 +74,9 @@ c Note that including cactus.h will also include AHFinder.h
 
 !     Calculate quantities in spherical coordinates.
       er = sqrt(xp**2 + yp**2 + zp**2)
-      if (er .eq. 0.d0) stop 'r=0 in AHFinder_calcsigma'
+      if (er .eq. 0.d0) then
+         call CCTK_WARN(0,'r=0 in AHFinder_calcsigma')
+      end if
      
 C     ********PATCH******************************************
 C     Keep away from the z-axis, artificially.
@@ -111,7 +113,7 @@ C     Other l>0.
 C     Initialize m = l:
 C     c1 = sqrt(2l+1) * sqrt(2l!) / 2^l / l! 
 C     Y_ll = c1 * sin(theta)^l 
-C     Y_ll' = c1 * l * sin(theta)^(l-1) * costheta
+C     Y_llp = c1 * l * sin(theta)^(l-1) * costheta
          m = l
          prod = - prod * sqrt(dble(2*l-1) / dble(2*l)) 
          c1 = sqrt(dble(2*l+1)) * prod
@@ -122,7 +124,7 @@ C     Recursion relations for the other m:
 C     c2 = sqrt[(2l+1)(2l-1)/(l^2-m^2)]
 C     c3 = sqrt[((l-1)^2-m^2)(2l=1)/(l^2-m^2)/(2l-3)]
 C     Y_lm = c2 Y_l-1,m cos(theta) + c3 Y_l-2,m
-C     Y_lm' = c2 Y_l-1,m' cos(theta) + c2 Y_l-1,m + c3 Y_l-2,m'
+C     Y_lmp = c2 Y_l-1,mp cos(theta) + c2 Y_l-1,m + c3 Y_l-2,mp
          m = l-1
          c2 = sqrt(dble((2*l+1)*(2*l-1)) 
      $        / dble(l**2-m**2)) 
@@ -168,7 +170,7 @@ C     to the basis cos/sin(mphi).
          end do
       else
 C     See above for a patch that avoids ever getting here.
-         stop 'ylmder2 cannot handle axis'
+          call CCTK_WARN(0,'ylmder2 cannot handle axis')
 C     For x = y = 0 the following are not defined.
 !         cosphi = 7.d77
 !         sinphi = 7.d77