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

/*@@
  @file      Mahc_corr_apply_bound.F
  @date      December 1999
  @author    Mark Miller
  @desc
     Apply the boundary update for the corrector step.
  @enddesc
@@*/

#include "cctk.h"
#include "cctk_Parameters.h"
#include "cctk_Arguments.h"


/*@@
  @routine    Mahc_corr_apply_bound
  @date       December 1999
  @author     Mark Miller
  @desc
     Apply the boundary update for the corrector step.
  @enddesc
  @calls
  @calledby
  @history
  @endhistory
@@*/

      subroutine Mahc_corr_apply_bound(CCTK_FARGUMENTS)
      implicit none
      DECLARE_CCTK_FARGUMENTS
      DECLARE_CCTK_PARAMETERS
      integer order,CCTK_Equals

      order = -1
      if (CCTK_EQUALS(hydro_order,'fluxlim')) then
         order = 2
      endif
      if (CCTK_EQUALS(hydro_order,'upwind')) then
         order = 1
      endif
c     if we only want 1st order, then we are done!
      if(order .lt. 2) then
         return
      endif

      call Mahc_Boundary(cctkGH,"gen_rel_perfect_fluid::Mahc_evol_vars")
      call Mahc_Boundary(cctkGH,"gen_rel_perfect_fluid::Mahc_prim_vars")

      return
      end