c  Mahc_Evolve: Numerical Evolution of the General Relativistic Hydro Equations
c  Copyright (C) 2001 Mark Miller

/*@@
  @file      Mahc_Boundary.F
  @date      December 1999
  @author    Mark Miller
  @desc
     Boundary condition for MAHC hydro variables.
  @enddesc
@@*/
#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"


/*@@
  @routine    Mahc_Boundary
  @date       December 1999
  @author     Mark Miller
  @desc
     Boundary condition for MAHC hydro variables.
  @enddesc
  @calls
  @calledby
  @history
  @endhistory
@@*/

      subroutine Mahc_Boundary(GH,name)
      implicit none
      DECLARE_CCTK_PARAMETERS
      CCTK_POINTER :: GH
      CHARACTER*(*) :: name

      integer, dimension(3) :: sw = 2
      integer :: CCTK_Equals,ierr

      CCTK_REAL, PARAMETER :: zero = 0.0
      CCTK_REAL, PARAMETER :: one  = 1.0

c     Apply symmetry boundary conditions.
      call CartSymGN(ierr,GH,name)

c     Apply flat boundaries.
      if (CCTK_EQUALS(bound,"flat")) then
         call BndFlatGN(ierr,GH,sw,name)
      endif

      return
      end subroutine Mahc_Boundary