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
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
|
C Wrapper for boostrotdata. Calls it and vectorini.
C Sets Cauchy data, lapse and shift, and what else is needed
C in the Bona-Masso formalism, at an initial time.
C $Header$
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
#include "cctk_Functions.h"
subroutine Exact__initialize(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
DECLARE_CCTK_PARAMETERS
DECLARE_CCTK_FUNCTIONS
integer i,j,k
integer nx,ny,nz
CCTK_REAL tt, xx, yy, zz
CCTK_REAL alpjunk, dtalpjunk, axjunk, ayjunk, azjunk,
$ betaxjunk, betayjunk, betazjunk,
$ dtbetaxjunk, dtbetayjunk, dtbetazjunk,
$ bxxjunk, bxyjunk, bxzjunk,
$ byxjunk, byyjunk, byzjunk,
$ bzxjunk, bzyjunk, bzzjunk
CCTK_REAL
$ dxgxxjunk, dxgyyjunk, dxgzzjunk,
$ dxgxyjunk, dxgyzjunk, dxgxzjunk,
$ dygxxjunk, dygyyjunk, dygzzjunk,
$ dygxyjunk, dygyzjunk, dygxzjunk,
$ dzgxxjunk, dzgyyjunk, dzgzzjunk,
$ dzgxyjunk, dzgyzjunk, dzgxzjunk
CCTK_REAL
$ exact_psi,
$ exact_psix, exact_psiy, exact_psiz,
$ exact_psixx, exact_psiyy, exact_psizz,
$ exact_psixy, exact_psiyz, exact_psixz
call CCTK_INFO('setting exact data on slice')
C Set conformal state
if (CCTK_EQUALS(metric_type, "static conformal")) then
conformal_state=1
if (CCTK_EQUALS(conformal_storage,"factor+derivs")) then
conformal_state = 2
else if
$ (CCTK_EQUALS(conformal_storage, "factor+derivs+2nd derivs"))
$ then
conformal_state = 3
end if
end if
C Note I assume time has been initialized to physical time.
C Set data pointwise.
nx = cctk_lsh(1)
ny = cctk_lsh(2)
nz = cctk_lsh(3)
C.omp parallel do private(
C.omp$ i, j, k,
C.omp$ tt, xx, yy, zz,
C.omp$ alpjunk, dtalpjunk, axjunk, ayjunk, azjunk,
C.omp$ betaxjunk, betayjunk, betazjunk,
C.omp$ dtbetaxjunk, dtbetayjunk, dtbetazjunk,
C.omp$ bxxjunk, bxyjunk, bxzjunk,
C.omp$ byxjunk, byyjunk, byzjunk,
C.omp$ bzxjunk, bzyjunk, bzzjunk,
C.omp$ dxgxxjunk, dxgyyjunk, dxgzzjunk,
C.omp$ dxgxyjunk, dxgyzjunk, dxgxzjunk,
C.omp$ dygxxjunk, dygyyjunk, dygzzjunk,
C.omp$ dygxyjunk, dygyzjunk, dygxzjunk,
C.omp$ dzgxxjunk, dzgyyjunk, dzgzzjunk,
C.omp$ dzgxyjunk, dzgyzjunk, dzgxzjunk,
C.omp$ exact_psi,
C.omp$ exact_psix, exact_psiy, exact_psiz,
C.omp$ exact_psixx, exact_psiyy, exact_psizz,
C.omp$ exact_psixy, exact_psiyz, exact_psixz)
do k=1,nz
do j=1,ny
do i=1,nx
tt = cctk_time
xx = x(i,j,k) - cctk_time * shift_add_x
yy = y(i,j,k) - cctk_time * shift_add_y
zz = z(i,j,k) - cctk_time * shift_add_z
C Initialize the psi of exact
C (also to tell the models about the conformal_state)
if (conformal_state .ne. 0) then
exact_psi = 1.0D0
else
exact_psi = 0.0D0
end if
exact_psix = 0.0D0
exact_psiy = 0.0D0
exact_psiz = 0.0D0
exact_psixx = 0.0D0
exact_psiyy = 0.0D0
exact_psizz = 0.0D0
exact_psixy = 0.0D0
exact_psiyz = 0.0D0
exact_psixz = 0.0D0
call Exact__Bona_Masso_data(
$ decoded_exact_model,
$ xx, yy, zz, tt,
$ gxx(i,j,k), gyy(i,j,k), gzz(i,j,k),
$ gxy(i,j,k), gyz(i,j,k), gxz(i,j,k),
$ kxx(i,j,k), kyy(i,j,k), kzz(i,j,k),
$ kxy(i,j,k), kyz(i,j,k), kxz(i,j,k),
$ exact_psi,
$ exact_psix, exact_psiy, exact_psiz,
$ exact_psixx, exact_psiyy, exact_psizz,
$ exact_psixy, exact_psiyz, exact_psixz,
$ dxgxxjunk, dxgyyjunk, dxgzzjunk,
$ dxgxyjunk, dxgyzjunk, dxgxzjunk,
$ dygxxjunk, dygyyjunk, dygzzjunk,
$ dygxyjunk, dygyzjunk, dygxzjunk,
$ dzgxxjunk, dzgyyjunk, dzgzzjunk,
$ dzgxyjunk, dzgyzjunk, dzgxzjunk,
$ alpjunk, dtalpjunk, axjunk, ayjunk, azjunk,
$ betaxjunk, betayjunk, betazjunk,
$ dtbetaxjunk, dtbetayjunk, dtbetazjunk,
$ bxxjunk, bxyjunk, bxzjunk,
$ byxjunk, byyjunk, byzjunk,
$ bzxjunk, bzyjunk, bzzjunk)
C Save the conformal factor if wanted
if (conformal_state .ne. 0) then
psi(i,j,k) = exact_psi
if (conformal_state .gt. 1) then
psix(i,j,k) = exact_psix
psiy(i,j,k) = exact_psiy
psiz(i,j,k) = exact_psiz
if (conformal_state .gt. 2) then
psixx(i,j,k) = exact_psixx
psiyy(i,j,k) = exact_psiyy
psizz(i,j,k) = exact_psizz
psixy(i,j,k) = exact_psixy
psiyz(i,j,k) = exact_psiyz
psixz(i,j,k) = exact_psixz
end if
end if
end if
end do
end do
end do
C Tell the code there is no need to treat the conformal factor
C as a separate field. That is, we have set the physical metric here.
c Commented out in einstein revamp, now Exact does not inherit anything
c about the conformal factor
c Now it does again (see above, knarf)
return
end
|