diff options
| author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-01-22 11:57:51 +0000 |
|---|---|---|
| committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2001-01-22 11:57:51 +0000 |
| commit | f753cf5efd24c88ceed4e015f4fbd7c47f140231 (patch) | |
| tree | 4c0bb32ba20f8e0b602bad10f6d9fa047d5ba455 /src/AHFinder.F | |
| parent | 3bee885b5616fe2b641a37ea68c3514972175a87 (diff) | |
Making the circumferences visible to all routines, plus getting rid
of goto statements.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@170 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder.F')
| -rw-r--r-- | src/AHFinder.F | 26 |
1 files changed, 12 insertions, 14 deletions
diff --git a/src/AHFinder.F b/src/AHFinder.F index 4fd6250..15522c4 100644 --- a/src/AHFinder.F +++ b/src/AHFinder.F @@ -38,7 +38,6 @@ CCTK_REAL atime CCTK_REAL rmx,r0 CCTK_REAL intexp_h,intexp2_h,intarea_h - CCTK_REAL circ_eq,meri_p1,meri_p2 CCTK_REAL aux,expin,expout,small1,small2 CCTK_REAL c0_temp(0:18) CCTK_REAL maxval(3),minval(3) @@ -58,7 +57,7 @@ ! ! horizon Horizon found? ! -! gauss Ouput Gaussian curvature? +! gauss Output Gaussian curvature? ! ! ahfso How often we want to find a horizon? ! ahfafter After how many iterations start looking for horizons? @@ -69,11 +68,7 @@ ! ! rmx Radius of largest sphere that fits in the grid. ! -! circ_eq Equatorial circumference. -! meri_p1 Length of meridian at phi=0. -! meri_p2 Length of meridian at phi=pi/2. -! -! mtype Type of surface. +! mtype Type of surface: ! 0. No surface found. ! 1. Integral of H^2 is positive out and negative in. ! 2. Integral of H^2 is negative out and positive in. @@ -419,22 +414,25 @@ zc = ahf_zc end if +! Check if the center of the expansion is +! out of the grid. + if ((xc.gt.xmx).or.(xc.lt.xmn)) then write(*,*) - write(*,*) 'xc is out of the grid, giving up.' - goto 5 + write(*,*) 'xc is out of the grid, setting it to zero.' + xc = 0.0D0 end if if ((yc.gt.ymx).or.(yc.lt.ymn)) then write(*,*) - write(*,*) 'yc is out of the grid, giving up.' - goto 5 + write(*,*) 'yc is out of the grid, setting it to zero.' + yc = 0.0D0 end if if ((zc.gt.zmx).or.(zc.lt.zmn)) then write(*,*) - write(*,*) 'zc is out of the grid, giving up.' - goto 5 + write(*,*) 'zc is out of the grid, setting it to zero.' + zc = 0.0D0 end if @@ -930,7 +928,7 @@ call AHFinder_int(CCTK_FARGUMENTS) if (ahf_gaussout.ne.0) then - call AHFinder_gau(CCTK_FARGUMENTS,circ_eq,meri_p1,meri_p2) + call AHFinder_gau(CCTK_FARGUMENTS) else circ_eq = 0.0D0 meri_p1 = 0.0D0 |