diff options
author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-11-15 15:00:36 +0000 |
---|---|---|
committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-11-15 15:00:36 +0000 |
commit | c4ba23db8d3690fde5907126d4e071be68216a44 (patch) | |
tree | 53ffd23384068f7b71acb6f94927fef70cd78b7b /src | |
parent | 76ed6b202981f27e241f5dff900a2c325dad9451 (diff) |
minor changes. Cleaning up a bit.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@259 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src')
-rw-r--r-- | src/AHFinder_pow.F | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/src/AHFinder_pow.F b/src/AHFinder_pow.F index 862f63e..149a474 100644 --- a/src/AHFinder_pow.F +++ b/src/AHFinder_pow.F @@ -21,12 +21,12 @@ implicit none - integer NCOM + logical error - CCTK_REAL, allocatable, dimension(:) :: PCOM,XICOM + integer NCOM CCTK_REAL XMIN - logical error + CCTK_REAL, allocatable, dimension(:) :: PCOM,XICOM ! *************** @@ -74,8 +74,8 @@ CCTK_REAL FUNC,DEL,FP,FPTT,FRET,FTOL,T,TOL1,ZEPS - CCTK_REAL, dimension(N) :: P,PT,PTT,XIT - CCTK_REAL, dimension(N,N) :: XI + CCTK_REAL, dimension(1:N) :: P,PT,PTT,XIT + CCTK_REAL, dimension(1:N,1:N) :: XI character*200 logf @@ -238,9 +238,10 @@ CCTK_REAL zero CCTK_REAL dx,dy,dz - CCTK_REAL, dimension(N) :: P,XI + CCTK_REAL, dimension(1:N) :: P,XI CCTK_REAL F1DIM + EXTERNAL F1DIM @@ -250,14 +251,13 @@ zero = 0.0d0 - allocate(PCOM(N),XICOM(N)) - PCOM = zero - XICOM= zero + allocate(PCOM(1:N),XICOM(1:N)) ! ************************* ! *** START ROUTINE *** ! ************************* + dx = cctk_delta_space(1) dy = cctk_delta_space(2) dz = cctk_delta_space(3) @@ -330,7 +330,7 @@ CCTK_REAL XX,FUNC - CCTK_REAL, dimension(NCOM) :: XT + CCTK_REAL, dimension(1:NCOM) :: XT ! ************************* @@ -585,7 +585,7 @@ ! Test for done here. TOL1 = 0.5D0*TOL*(dabs(FP)+dabs(FX)) + ZEPS - if (dabs(FP-FX).le.TOL1)goto 3 + if (dabs(FP-FX).le.TOL1) goto 3 ! Construct a trial parabolic fit. @@ -689,7 +689,7 @@ CCTK_REAL zero,half,one CCTK_REAL aux - CCTK_REAL, dimension(N) :: P + CCTK_REAL, dimension(1:N) :: P character*200 logf |