diff options
Diffstat (limited to 'src/AHFinder_gau.F')
| -rw-r--r-- | src/AHFinder_gau.F | 16 |
1 files changed, 6 insertions, 10 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F index 72c4359..3261bc7 100644 --- a/src/AHFinder_gau.F +++ b/src/AHFinder_gau.F @@ -7,7 +7,7 @@ c Find gaussian curvature of surface. The gaussian c curvature is defined as -R/2, where R is the Ricci c scalar of the induced 2-geometry of the surface. c -c As an extra `goody', this routine also integrates +c As an extra "goody", this routine also integrates c the equatorial circumference of the horizon and c two polar circumferences at phi=0 and phi=pi/2. c @enddesc @@ -38,7 +38,6 @@ c Note that including cactus.h will also include AHFinder.h integer gxx_index,gyy_index,gzz_index,gxy_index,gxz_index,gyz_index integer ahfgauss_index integer handle,x_index,y_index,z_index - CCTK_INT ierror CCTK_REAL maxval(3),minval(3) INTEGER istat @@ -303,8 +302,7 @@ c Note that including cactus.h will also include AHFinder.h dphi = 0.25D0*dphi else if (refx) then dphi = half*dphi - print *, 'This combination of symmetries has not' - print *, 'been properly implemented yet!' + call CCTK_INFO('This combination of symmetries has not been properly implemented yet!') else if (refy) then dphi = half*dphi end if @@ -485,10 +483,10 @@ c Note that including cactus.h will also include AHFinder.h ! Synchronize. -#ifdef MPI + call CCTK_SyncGroup(cctkGH,"ahfinder::ahfinder_gauss") ! call synconefunc(ahfgauss) -#endif + ! ************************************* ! *** FIND INTERPOLATING POINTS *** @@ -530,7 +528,7 @@ c Note that including cactus.h will also include AHFinder.h ! 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 +! l s even if I do not want some terms. This is ! because in order to use the recursion relations I ! need all polynomials. @@ -593,14 +591,12 @@ c Note that including cactus.h will also include AHFinder.h ! below). #ifdef MPI - call mpi_allreduce(error1,i,1,MPI_INTEGER, . MPI_SUM,MPI_COMM_WORLD,ierror) error1 = i call mpi_allreduce(error2,i,1,MPI_INTEGER, . MPI_SUM,MPI_COMM_WORLD,ierror) error2 = i - #endif ! If there was an error then return from the subroutine @@ -818,7 +814,7 @@ c Note that including cactus.h will also include AHFinder.h ! This was the old calculation of the gaussian curvature. ! It is correct as far as I know, but it depends on taking ! second derivatives of the interpolated metric, and this -! seems not to behave well numerically. I didn't want to +! seems not to behave well numerically. I did not want to ! delete the code since it might be useful in the future, ! so I just jump over it. |