blob: 43f5b9170ae6dfa4a982b2e07b6dfcc26133fb11 (
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 DYDB_guts.h
@date Jun 98
@author Gabrielle Allen
@desc
Macro to calculate the first derivatives of the
shift with respect to y
The macro is defined in terms of standard variables in
@seefile DYDB_declare.h
@enddesc
@@*/
#ifndef DYDB_GUTS
#define DYDB_GUTS
#ifdef FCODE
DYDB_OO2DY = 1D0/(2D0*DYDB_DY)
if (admmacros_advectiony .eq. 0) then /* CENTER */
DYDB_DYDBX = DYDB_OO2DY*(DYDB_BX_JP - DYDB_BX_JM)
else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */
DYDB_DYDBX = (DYDB_BX_JP - DYDB_BX)/DYDB_DY
else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */
DYDB_DYDBX = (DYDB_BX - DYDB_BX_JM)/DYDB_DY
else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */
DYDB_DYDBX = - 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BX - 4.0D0*DYDB_BX_JP + DYDB_BX_JPP)
else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */
DYDB_DYDBX = 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BX - 4.0D0*DYDB_BX_JM + DYDB_BX_JMM)
end if
if (admmacros_advectiony .eq. 0) then /* CENTER */
DYDB_DYDBY = DYDB_OO2DY*(DYDB_BY_JP - DYDB_BY_JM)
else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */
DYDB_DYDBY = (DYDB_BY_JP - DYDB_BY)/DYDB_DY
else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */
DYDB_DYDBY = (DYDB_BY - DYDB_BY_JM)/DYDB_DY
else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */
DYDB_DYDBY = - 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BY - 4.0D0*DYDB_BY_JP + DYDB_BY_JPP)
else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */
DYDB_DYDBY = 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BY - 4.0D0*DYDB_BY_JM + DYDB_BY_JMM)
end if
if (admmacros_advectiony .eq. 0) then /* CENTER */
DYDB_DYDBZ = DYDB_OO2DY*(DYDB_BZ_JP - DYDB_BZ_JM)
else if (admmacros_advectiony .eq. 1) then /* UPWIND1 */
DYDB_DYDBZ = (DYDB_BZ_JP - DYDB_BZ)/DYDB_DY
else if (admmacros_advectiony .eq. -1) then /* UPWIND1 */
DYDB_DYDBZ = (DYDB_BZ - DYDB_BZ_JM)/DYDB_DY
else if (admmacros_advectiony .eq. 2) then /* UPWIND2 */
DYDB_DYDBZ = - 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BZ - 4.0D0*DYDB_BZ_JP + DYDB_BZ_JPP)
else if (admmacros_advectiony .eq. -2) then /* UPWIND2 */
DYDB_DYDBZ = 0.5D0/DYDB_DY \
*(3.0D0*DYDB_BZ - 4.0D0*DYDB_BZ_JM + DYDB_BZ_JMM)
end if
#endif
#ifdef CCODE
DYDB_OO2DY = 1/(2*cctkGH->cctk_delta_space[1]);
DYDB_DYDBX = DYDB_OO2DY*(DYDB_BX_JP - DYDB_BX_JM);
DYDB_DYDBY = DYDB_OO2DY*(DYDB_BY_JP - DYDB_BY_JM);
DYDB_DYDBZ = DYDB_OO2DY*(DYDB_BZ_JP - DYDB_BZ_JM);
#endif
#endif
|