diff options
author | lars <lars@89daf98e-ef62-4674-b946-b8ff9de2216c> | 1999-07-21 09:50:23 +0000 |
---|---|---|
committer | lars <lars@89daf98e-ef62-4674-b946-b8ff9de2216c> | 1999-07-21 09:50:23 +0000 |
commit | c5c68fbff7f9b27f7bb7439b6fcc77879823dcd9 (patch) | |
tree | f8ddb8690186d9961b73172736a300fcdf8047c3 /src/AHFinder_find3.F | |
parent | 8fafcd075208081ec07cb07155b40939889cf234 (diff) |
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
Diffstat (limited to 'src/AHFinder_find3.F')
-rw-r--r-- | src/AHFinder_find3.F | 68 |
1 files changed, 68 insertions, 0 deletions
diff --git a/src/AHFinder_find3.F b/src/AHFinder_find3.F new file mode 100644 index 0000000..951bb4b --- /dev/null +++ b/src/AHFinder_find3.F @@ -0,0 +1,68 @@ +c/*@@ +c @file AHFinder_find3.F +c @date March 1999 +c @author Lars Nerger +c @desc +c Muliplicate gridfunctions for final grid function +c @enddesc +c@@*/ + +c Note that including cactus.h will also include AHFinder.h +!#include "cactus.h" +#include "cctk.h" +#include "cctk_parameters.h" +#include "cctk_arguments.h" + +! To output single gridfunctions when searching for 3 horizons +! the gridfunctions are here multiplicated. + + subroutine AHFinder_find3(CCTK_FARGUMENTS,mtype) + + use AHFinder_dat + + implicit none + + DECLARE_CCTK_FARGUMENTS + DECLARE_CCTK_PARAMETERS + + integer mtype + + REAL one + +! Description of variables +! +! mtype Type of surface: see AHFinder.F + + +! ********************** +! *** DEFINE one *** +! ********************** + + one = 1.0D0 + + +! ************************* +! *** START ROUTINE *** +! ************************* + + if ((mfind.eq.0).and.(mtype.eq.1)) then + ahfgrid3 = ahfgrid + ahf_exp3 = ahf_exp + else if ((mfind.eq.0).and.(mtype.ne.1)) then + ahfgrid3 = one + ahf_exp3 = one + else if ((mfind.eq.1).and.(mtype.eq.1)) then + ahfgrid3 = ahfgrid*ahfgrid3 + ahf_exp3 = ahf_exp*ahf_exp3 + else if ((mfind.eq.2).and.(mtype.eq.1)) then + ahfgrid3 = ahfgrid*ahfgrid3 + ahf_exp3 = ahf_exp*ahf_exp3 + end if + + +! *************** +! *** END *** +! *************** + + end subroutine AHFinder_find3 + |