blob: 0933cd79dd3dbd3064e6f972363283bdf5ef8fa3 (
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
|
C $Header$
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"
#include "cctk_Functions.h"
#include "Exact.inc"
subroutine Exact__slice_initialize(CCTK_ARGUMENTS)
implicit none
DECLARE_CCTK_ARGUMENTS
c #define-ing the symbol EXACT_NO_F90 will turn this subroutine into a no-op
#ifndef EXACT_NO_F90
DECLARE_CCTK_PARAMETERS
C Initialize position of slice in exact solution spacetime.
C For now the spatial coordinates are identified with those of the
C exact solution, and the exact solution time is a Gaussian.
slicex = x
slicey = y
slicez = z
slicet = slice_gauss_ampl*exp(-(x**2 + y**2 + z**2)
$ /slice_gauss_width**2)
C Calculate Cauchy data and dx^A/dt.
call Exact__slice_data(CCTK_ARGUMENTS)
C This tells the code that what we have set is the physical metric,
C not a conformally rescaled one.
c Commented out in einstein revamp, now Exact doesnot inherit anything
c about the conformal factor
c conformal_state = 0
c psi = 1.0D0
#endif
return
end
|