blob: 295da7b5c16a6908980489787dc3ebb2292002d6 (
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
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
|
/*@@
@header DXDB_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the first derivatives of the
shift with respect to x
The macro is defined in terms of standard variables in
@seefile DXDB_declare.h
@enddesc
@@*/
#ifndef DXDB_GUTS
#define DXDB_GUTS
#ifdef FCODE
DXDB_OO2DX = 1D0/(2D0*DXDB_DX)
if (admmacros_advectionx .eq. 0) then /* CENTER */
DXDB_DXDBX = DXDB_OO2DX*(DXDB_BX_IP - DXDB_BX_IM)
else if (admmacros_advectionx .eq. 1) then /* UPWIND1 */
DXDB_DXDBX = (DXDB_BX_IP - DXDB_BX)/DXDB_DX
else if (admmacros_advectionx .eq. -1) then /* UPWIND1 */
DXDB_DXDBX = (DXDB_BX - DXDB_BX_IM)/DXDB_DX
else if (admmacros_advectionx .eq. 2) then /* UPWIND2 */
DXDB_DXDBX = - 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BX - 4.0D0*DXDB_BX_IP + DXDB_BX_IPP)
else if (admmacros_advectionx .eq. -2) then /* UPWIND2 */
DXDB_DXDBX = 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BX - 4.0D0*DXDB_BX_IM + DXDB_BX_IMM)
end if
if (admmacros_advectionx .eq. 0) then /* CENTER */
DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM)
else if (admmacros_advectionx .eq. 1) then /* UPWIND1 */
DXDB_DXDBY = (DXDB_BY_IP - DXDB_BY)/DXDB_DX
else if (admmacros_advectionx .eq. -1) then /* UPWIND1 */
DXDB_DXDBY = (DXDB_BY - DXDB_BY_IM)/DXDB_DX
else if (admmacros_advectionx .eq. 2) then /* UPWIND2 */
DXDB_DXDBY = - 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BY - 4.0D0*DXDB_BY_IP + DXDB_BY_IPP)
else if (admmacros_advectionx .eq. -2) then /* UPWIND2 */
DXDB_DXDBY = 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BY - 4.0D0*DXDB_BY_IM + DXDB_BY_IMM)
end if
if (admmacros_advectionx .eq. 0) then /* CENTER */
DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM)
else if (admmacros_advectionx .eq. 1) then /* UPWIND1 */
DXDB_DXDBZ = (DXDB_BZ_IP - DXDB_BZ)/DXDB_DX
else if (admmacros_advectionx .eq. -1) then /* UPWIND1 */
DXDB_DXDBZ = (DXDB_BZ - DXDB_BZ_IM)/DXDB_DX
else if (admmacros_advectionx .eq. 2) then /* UPWIND2 */
DXDB_DXDBZ = - 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BZ - 4.0D0*DXDB_BZ_IP + DXDB_BZ_IPP)
else if (admmacros_advectionx .eq. -2) then /* UPWIND2 */
DXDB_DXDBZ = 0.5D0/DXDB_DX \
*(3.0D0*DXDB_BZ - 4.0D0*DXDB_BZ_IM + DXDB_BZ_IMM)
end if
#endif
#ifdef CCODE
DXDB_OO2DX = 1/(2*cctkGH->cctk_delta_space[0]);
DXDB_DXDBX = DXDB_OO2DX*(DXDB_BX_IP - DXDB_BX_IM);
DXDB_DXDBY = DXDB_OO2DX*(DXDB_BY_IP - DXDB_BY_IM);
DXDB_DXDBZ = DXDB_OO2DX*(DXDB_BZ_IP - DXDB_BZ_IM);
#endif
#endif
|