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
69
70
71
72
73
74
75
|
# Schedule definitions for thorn AHFinder
#c/*@@
#c @date July 1999
#c @author Lars Nerger
#c @@*/
schedule AHFinder_SetSym at CCTK_WRAGH
{
LANG: Fortran
OPTIONS: global
} "Set symmetries for AHFinder grid functions"
schedule AHFinder_Startup at CCTK_STARTUP
{
LANG: C
OPTIONS: global
} "Register AHFinder as an IO Method"
schedule AHFinder_ParamCheck at PARAMCHECK
{
LANG: Fortran
OPTIONS: global
} "Check for physical or conformal metric"
schedule AHFinder_InitOutput at CCTK_INITIAL
{
LANG: Fortran
} "Create output files, write headers"
if (ahf_active)
{
STORAGE: hole1_bounds,hole2_bounds,hole3_bounds,ahf_centroid
STORAGE: out_scalars
STORAGE: ahfinder_gaussian
if (ahf_lmax > 0)
{
STORAGE: out_1d_legen, out_2d_legen
}
if (ahf_persists)
{
STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid
STORAGE: ahfmask
schedule ahfinder before driftcorrect at CCTK_ANALYSIS
{
LANG: Fortran
STORAGE: triggergroup
TRIGGERS: triggergroup
} "Call apparent horizon finder with persisting grid functions"
}
else if (CCTK_Equals (ahf_mask,"off"))
{
schedule ahfinder before driftcorrect at CCTK_ANALYSIS
{
LANG: Fortran
STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid
STORAGE: ahfmask, triggergroup
TRIGGERS:triggergroup
} "Call apparent horizon finder"
}
else
{
STORAGE: ahfmask
STORAGE: hole1_bounds,hole2_bounds,hole3_bounds
schedule ahfinder before driftcorrect at CCTK_ANALYSIS
{
LANG: Fortran
STORAGE: ahfindergrid,ahfinderexp,ahfgradient,ahfinder_gauss,find3grid, triggergroup
TRIGGERS: triggergroup
} "Call apparent horizon finder with persisting mask"
}
}
|