diff options
| author | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-05 13:38:47 +0000 |
|---|---|---|
| committer | miguel <miguel@89daf98e-ef62-4674-b946-b8ff9de2216c> | 2000-05-05 13:38:47 +0000 |
| commit | 352260d6792d2594fb65b26f9973568a255476c7 (patch) | |
| tree | d7ec13a70b17402213fc209fa0f1cec6b966422d /src/AHFinder_int.F | |
| parent | 34e9611eb9a48c55ae5e5f872c42afb5b07e95a1 (diff) | |
Putting back "integer" istead of CCTK_INT in most places. I have been
advised by Gerd that it is better to keep "integer" for things that
are not communicated.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinder/trunk@96 89daf98e-ef62-4674-b946-b8ff9de2216c
Diffstat (limited to 'src/AHFinder_int.F')
| -rw-r--r-- | src/AHFinder_int.F | 36 |
1 files changed, 20 insertions, 16 deletions
diff --git a/src/AHFinder_int.F b/src/AHFinder_int.F index 82ff1bf..3a00d5c 100644 --- a/src/AHFinder_int.F +++ b/src/AHFinder_int.F @@ -25,12 +25,13 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS - CCTK_INT i,j,l,m - CCTK_INT npt,npp - CCTK_INT l_ntheta,l_nphi,theta0,phi0 - CCTK_INT npoints - CCTK_INT ierror,auxi - CCTK_INT handle + integer i,j,l,m + integer npt,npp + integer l_ntheta,l_nphi,theta0,phi0 + integer npoints + integer ierror + + CCTK_INT handle,auxi CCTK_REAL LEGEN CCTK_REAL theta,phi,xp,yp,zp,rp @@ -136,14 +137,6 @@ ! ! da Area element. -! Get the reduction handle for the sum operation. - - call CCTK_ReductionArrayHandle(reduction_handle,"sum") - - if (reduction_handle.lt.0) then - call CCTK_WARN(1,"Cannot get reduction handle for SUM operation.") - endif - ! ************************** ! *** DEFINE NUMBERS *** @@ -157,12 +150,23 @@ three = 3.0D0 four = 4.0D0 - pi = 3.141592654D0 + pi = acos(-one) myproc = CCTK_MyProc(cctkGH) nprocs = CCTK_nProcs(cctkGH) +! **************************************** +! *** GET REDUCTION HANDLE FOR SUM *** +! **************************************** + + call CCTK_ReductionArrayHandle(reduction_handle,"sum") + + if (reduction_handle.lt.0) then + call CCTK_WARN(1,"Cannot get reduction handle for SUM operation.") + end if + + ! ********************************* ! *** INITIALIZE PARAMETERS *** ! ********************************* @@ -323,7 +327,7 @@ ! tells my on which "theta" bin I am, while myproc/npt ! tells me on which "phi" bin I am. - i = mod(int(myproc),int(npt)) + i = mod(myproc,npt) j = myproc/npt theta0 = i*l_ntheta |