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