blob: 951bb4b4987a3a2acc339cb1be8c8b64f4703232 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
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
|