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
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
# Interface definition for thorn EHFinder
# $Header$
implements: ehfinder
inherits: grid admbase coordgauge staticconformal spacemask boundary
USES INCLUDE: Boundary.h
USES INCLUDE: MoL.h
public:
CCTK_REAL level_set TYPE=GF TIMELEVELS=2
{
f
} "The scalar level set function that defines the null surface"
CCTK_REAL slevel_set TYPE=GF TIMELEVELS=1
{
sf
} "Source for the level set function"
CCTK_REAL dlevel_set TYPE=GF TIMELEVELS=1
{
dfx, dfy, dfz, dfsq
} "Derivatives of the level set function"
CCTK_REAL ftmp_set TYPE=GF TIMELEVELS=1
{
ftmp, sftmp
} "temporary variables used in pde re-parametrization"
CCTK_REAL level_bak TYPE=GF TIMELEVELS=1
{
fbak
} "Temporary placeholder for the level set during re-parametrization"
CCTK_INT eh_mask_all TYPE=GF TIMELEVELS=1
{
eh_mask, tm_mask
} "Masks to define active cells"
CCTK_INT mask_bak TYPE=GF TIMELEVELS=1
{
eh_mask_bak
} "Temporary placeholder for the mask during re-parametrization"
CCTK_INT surface_index TYPE=GF TIMELEVELS=1
{
sc
}
CCTK_INT surface_counter TYPE=SCALAR
CCTK_INT points_counter TYPE=SCALAR
CCTK_INT more_surfaces TYPE=SCALAR
CCTK_INT more_points TYPE=SCALAR
CCTK_INT integrate_counter TYPE=SCALAR
CCTK_REAL surface_arrays TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
{
ctheta, cphi, rsurf, sintheta, costheta, sinphi, cosphi,
drdtheta, drdphi, da, gtt, gtp, gpp, weights
} "Grid arrays for points on the surface"
CCTK_REAL surface_tmp_arrays TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
{
drsurf, interp_x, interp_y, interp_z, f_interp,
dfdx_interp, dfdy_interp, dfdz_interp
} "Temporary grid arrays for finding points on the surface"
CCTK_REAL center_arrays TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=3 DISTRIB=CONSTANT
{
center
} "The cartesian location of the center for the spherical coordinate system"
CCTK_INT surface_int_array TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
{
n_since_last_reduction
} "Temporary integer grid array for finding points on the surface"
CCTK_REAL interp_metric_arrays TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
{
gxxi, gxyi, gxzi, gyyi, gyzi, gzzi, psii
} "Arrays for holding the interpolated metric and conformal factor"
CCTK_REAL integrate_tmp_array TYPE=ARRAY DIM=2 TIMELEVELS=1 SIZE=ntheta,nphi GHOSTSIZE=n_array_ghosts,n_array_ghosts DISTRIB=DEFAULT
{
int_tmp
}
CCTK_INT rep_mask TYPE=GF TIMELEVELS=1
CCTK_INT re_param_control_pde TYPE=SCALAR
CCTK_INT re_param_control_approx TYPE=SCALAR
#CCTK_REAL eh_area TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number GHOSTSIZE=0 DISTRIB=DEFAULT
CCTK_REAL eh_area TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT
#CCTK_REAL eh_area[maximum_surface_number] TYPE=SCALAR
#CCTK_REAL eh_area2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number GHOSTSIZE=0 DISTRIB=DEFAULT #"The area of the event horizon(s)"
CCTK_REAL eh_area2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT #"The area of the event horizon(s)"
#CCTK_REAL eh_area2[maximum_surface_number] TYPE=SCALAR
#CCTK_REAL eh_area3 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT
CCTK_REAL eh_centroid TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT
{
eh_centroid_x, eh_centroid_y, eh_centroid_z
} "The centroid of the event horizon(s)"
CCTK_REAL eh_centroid2 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT
{
eh_centroid2_x, eh_centroid2_y, eh_centroid2_z
} "The centroid of the event horizon(s)"
#CCTK_REAL eh_centroid3 TYPE=ARRAY DIM=1 TIMELEVELS=1 SIZE=maximum_surface_number DISTRIB=CONSTANT
#{
# eh_centroid3_x, eh_centroid3_y, eh_centroid3_z
#} "The centroid of the event horizon(s)"
|