src/macro/DA_guts.h


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

/*@@
  @header   DA_guts.h
  @date     Jul 98
  @author   Gabrielle Allen
  @desc
  Macro to calculate all first spatial derivative of lapse
  @enddesc
@@*/

#ifndef DA_GUTS
#define DA_GUTS

#ifdef FCODE 

      DA_OO2DX = 1D0/(2D0*DA_DX)
      DA_OO2DY = 1D0/(2D0*DA_DY)
      DA_OO2DZ = 1D0/(2D0*DA_DZ)

      DA_DXDA = DA_OO2DX*(DA_A_IP - DA_A_IM)
      DA_DYDA = DA_OO2DY*(DA_A_JP - DA_A_JM)
      DA_DZDA = DA_OO2DZ*(DA_A_KP - DA_A_KM)

#endif

#ifdef CCODE

      DA_OO2DX = 1/(2*DA_DX);
      DA_OO2DY = 1/(2*DA_DY);
      DA_OO2DZ = 1/(2*DA_DZ);

      DA_DXDA = DA_OO2DX*(DA_A_IP - DA_A_IM);
      DA_DYDA = DA_OO2DY*(DA_A_JP - DA_A_JM);
      DA_DZDA = DA_OO2DZ*(DA_A_KP - DA_A_KM);

#endif

#endif