aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authormiguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c>2000-05-03 10:55:47 +0000
committermiguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c>2000-05-03 10:55:47 +0000
commitcd0d8c667d3faa3a5ea1450ee626b07dcfaa3c74 (patch)
tree0aff9b7f9ef75cd3970ba052668ae69c93a1b569
parent77a660590f49e7b3693eb797651603f8e4c26a3b (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.F39
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 ***