This commit was generated by cvs2svn to compensate for changes in r2, which

included commits to RCS files with non-trunk default branches.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@3 89daf98e-ef62-4674-b946-b8ff9de2216c

1 files changed, 142 insertions, 0 deletions

diff --git a/src/AHFinder_dat.F b/src/AHFinder_dat.F
new file mode 100644
index 0000000..efcc6ae
--- /dev/null
+++ b/src/AHFinder_dat.F
@@ -0,0 +1,142 @@
+c/*@@
+c  @file      AHFinder_dat.F
+c  @date      April 1998
+c  @author    Miguel Alcubierre
+c  @desc 
+c             F90 module containing data to be shared by
+c             the different subroutines.
+c
+c             For the meaning of the spectral components of
+c             the different flows see:
+c
+c             C. Gundlach, Phys. Rev. D 57, 863 (1998).
+c
+c  @enddesc 
+c@@*/
+#include "cctk.h"
+      module AHFinder_dat
+
+      implicit none
+
+      logical firstfun,firstleg,firstint
+      logical nonaxi,refx,refy,refz,cartoon
+      logical offset,wander
+      logical minarea,flow,find_trapped_surface
+      logical logfile,verbose,veryver,guessverbose
+      logical interror,guessold
+      logical find3,sloppy,inner,guess_absmin,manual_guess
+
+      CCTK_INT lmax,stepx,stepy,stepz
+      CCTK_INT ntheta,nphi
+      CCTK_INT nfile
+      CCTK_INT interror1,interror2
+      CCTK_INT mfind,ahf_ncall
+      CCTK_INT nx,ny,nz
+      integer nprocs,myproc
+
+      
+      CCTK_REAL dx,dy,dz
+      CCTK_REAL xc,yc,zc
+      CCTK_REAL xc_old,yc_old,zc_old
+      CCTK_REAL intexp,intexp2,intexpdel2,intarea
+
+      CCTK_REAL hw,cw,nw
+      CCTK_REAL inside_min_count,inside_min_neg_count
+      CCTK_REAL dx,dy,dz
+
+      CCTK_REAL, allocatable, dimension (:)   :: c0,c0_old
+      CCTK_REAL, allocatable, dimension (:,:) :: cc,cs,cc_old,cs_old
+
+      CCTK_REAL, allocatable, dimension (:)   :: hflow0,cflow0,nflow0
+      CCTK_REAL, allocatable, dimension (:,:) :: hflowc,cflowc,nflowc
+      CCTK_REAL, allocatable, dimension (:,:) :: hflows,cflows,nflows
+
+
+      character*200 filestr
+
+      data mfind / 0 /
+      data ahf_ncall / 0 /
+
+      save mfind,ahf_ncall
+      save c0,cc,cs
+
+!     Description of variables:
+!
+!     firstfun   First call to function FUNC?
+!     firstleg   First call to function LEGEN?
+!     firstint   First call to subroutine AHFinder_int?
+!
+!     nonaxi     Is the data non-axisymmetric?
+!
+!     refx       Reflection symmetry x->-x?
+!
+!     refy       Reflection symmetry y->-y?
+!
+!     refz       Reflection symmetry z->-z?
+!
+!     cartoon    Are we using cartoon_2d?
+!
+!     inner     Look for inner horizon?
+!
+!     sloppy    Sloppy initial guess?
+!
+!     offset     Is the center offset from the origin?
+!     wander     Do we allow the center to wander?
+!
+!     minarea    Do we want to minimize area?
+!
+!     flow       Use fast flow algorithm instead of minimization?
+!
+!     logfile    Write log file?
+!
+!     verbose    Write messages to screen?
+!
+!     veryver    Write lots of messages to screen?
+!
+!     guessverbose   Write info on initial guess points?
+!
+!     lmax       Maximum number of terms in theta expansion.
+!
+!     stepx      1 if (refx=1), 0 otherwise.
+!     stepy      1 if (refy=1), 0 otherwise.
+!     stepz      1 if (refz=1), 0 otherwise.
+!
+!     ntheta     Number of subdivisions in theta.
+!     nphi       Number of subdivisions in phi.
+!
+!     xc         Coordinate x of centre of surface.
+!     yc         Coordinate y of centre of surface.
+!     zc         Coordinate z of centre of surface.
+!
+!     intexp     Integral of expansion over the surface.
+!     intexp2    Integral of expansion squared over the surface.
+!     intarea    Area of the surface.
+!
+!     c0(l)      Coefficients of LEGEN(l,0,cos(theta))
+!     cc(l,m)    Coefficients of LEGEN(l,m,cos(theta))*cos(m*phi)
+!     cs(l,m)    Coefficients of LEGEN(l,m,cos(theta))*sin(m*phi)
+!
+!     hflow*     Spectral components for `H' flow.
+!     cflow*     Spectral components for `C' flow.
+!     nflow*     Spectral components for `N' flow.
+!
+!     interror1  Different from zero if radius is negative.
+!     interror2  Different from zero if outside computational domain.
+!     interror   True if either interror1 or interro2 are non-zero.
+!
+!     inside_min_count        Number of elements in integral (should always
+!                             be ntheta*nphi)
+!     inside_min_neg_count    Number of elements in integral that have
+!                             negative expansion.
+!
+!     find3      Look for 3 horizons?
+!     mfind
+
+
+!     ***************
+!     ***   END   ***
+!     ***************
+
+      end module AHFinder_dat
+
+