diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-03 10:55:47 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-03 10:55:47 +0000 |
commit | cd0d8c667d3faa3a5ea1450ee626b07dcfaa3c74 (patch) | |
tree | 0aff9b7f9ef75cd3970ba052668ae69c93a1b569 | |
parent | 77a660590f49e7b3693eb797651603f8e4c26a3b (diff) |
Minor changes.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@88 89daf98e-ef62-4674-b946-b8ff9de2216c
-rw-r--r-- | src/AHFinder_fun.F | 39 |
1 files changed, 5 insertions, 34 deletions
diff --git a/src/AHFinder_fun.F b/src/AHFinder_fun.F index 4eac146..e386d1d 100644 --- a/src/AHFinder_fun.F +++ b/src/AHFinder_fun.F @@ -23,8 +23,6 @@ CCTK_INT i,j,k CCTK_INT l,m - CCTK_INT istat - CCTK_INT reduction_handle CCTK_REAL LEGEN CCTK_REAL xp,yp,zp,rp @@ -39,32 +37,20 @@ ! *** DEFINE NUMBERS *** ! ************************** - zero = 0.0D0 - half = 0.5D0 - one = 1.0D0 - two = 2.0D0 - pi = acos(-one) + zero = 0.0D0 + half = 0.5D0 + one = 1.0D0 + two = 2.0D0 + pi = acos(-one) halfpi = half*pi twopi = two*pi -! *************************** -! *** GRID PARAMETERS *** -! *************************** - - nx = cctk_lsh(1) - ny = cctk_lsh(2) - nz = cctk_lsh(3) - - ! ********************************* ! *** FIND HORIZON FUNCTION *** ! ********************************* - rhmax = zero - rhmin = 1.0D10 - do k=1,nz do j=1,ny do i=1,nx @@ -146,25 +132,10 @@ STOP end if -! Maximum and minimum radius. - - if (rhmax.lt.aux1) rhmax = aux1 - if (rhmin.gt.aux1) rhmin = aux1 - end do end do end do - call CCTK_ReductionArrayHandle(reduction_handle,"maximum") - call CCTK_ReduceLocScalar(ierr,cctkGH,-1,reduction_handle, - . rhmax,red_tmp,CCTK_VARIABLE_REAL) - rhmax = aux1 - - call CCTK_ReductionArrayHandle(reduction_handle,"minimum") - call CCTK_ReduceLocScalar(ierr,cctkGH,-1,reduction_handle, - . rhmin,red_tmp,CCTK_VARIABLE_REAL) - rhmin = aux1 - ! *************** ! *** END *** |