src/macro/DYDB_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
39

/*@@
  @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)
    
      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

#ifdef CCODE

      DYDB_OO2DY = 1/(2*GH->dy0);
    
      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