* remove GZ patch system code

* test suites pass



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@223 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45

1 files changed, 0 insertions, 156 deletions

diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90
index c2e418f..9e95adc 100644
--- a/src/GRHydro_Boundaries.F90
+++ b/src/GRHydro_Boundaries.F90
@@ -257,159 +257,3 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 end subroutine GRHydro_Boundaries
 
 
- /*@@
-   @routine    GRHydro_OutflowBoundaries
-   @date       Tue May 17 16:58:06 2005
-   @author     Ian Hawke
-   @desc 
-   Set outflow boundaries over only part of the domain.
-   This is designed to be used with GZPatchSystem.
-   @enddesc 
-   @calls     
-   @calledby   
-   @history 
- 
-   @endhistory 
-
-@@*/
-
-subroutine GRHydro_OutflowBoundaries(CCTK_ARGUMENTS)
-  
-  implicit none
-  
-  DECLARE_CCTK_ARGUMENTS
-  DECLARE_CCTK_PARAMETERS
-  DECLARE_CCTK_FUNCTIONS
-  
-  integer, dimension(3) :: sw
-  integer :: ierr
-  integer :: i,j,k
-  
-  CCTK_REAL, dimension(3) :: posn
-
-  CCTK_REAL :: det,psi4pt
-
-  sw = GRHydro_stencil
-
-  if (r(1,1,1) < r(1,1,cctk_lsh(3))) then
-
-    if (cctk_bbox(6) .ne. 0) then
-
-      do k = cctk_lsh(3) - sw(3), cctk_lsh(3)
-        do j = 1, cctk_lsh(2)
-          do i = 1, cctk_lsh(1)
-            
-            posn(1) = x(i,j,k)
-            posn(2) = y(i,j,k)
-            posn(3) = z(i,j,k)
-            
-            if (dot_product(outflowboundary_normal, posn) > 0.d0) then
-              
-              rho(i,j,k)        = rho(i,j,k-1)
-              velx(i,j,k)       = velx(i,j,k-1)
-              vely(i,j,k)       = vely(i,j,k-1)
-              velz(i,j,k)       = velz(i,j,k-1)
-              eps(i,j,k)        = eps(i,j,k-1)
-              press(i,j,k)      = press(i,j,k-1)
-              w_lorentz(i,j,k)  = w_lorentz(i,j,k-1)
-              if(evolve_mhd.ne.0) then
-                 Bvecx(i,j,k)       = Bvecx(i,j,k-1)
-                 Bvecy(i,j,k)       = Bvecy(i,j,k-1)
-                 Bvecz(i,j,k)       = Bvecz(i,j,k-1)
-                 if(clean_divergence.ne.0) then
-                    psidc(i,j,k)    = psidc(i,j,k-1)
-                 endif
-              endif
-  
-              psi4pt = 1.0d0
-              det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\
-                   gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))
-              
-              if(evolve_mhd.ne.0) then
-                 call prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),&
-                      psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),&
-                      psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),&
-                      det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),&
-                      tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),&
-                      rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),&
-                      eps(i,j,k),press(i,j,k),w_lorentz(i,j,k))
-              else
-                 call prim2con(GRHydro_eos_handle,psi4pt*gxx(i,j,k),&
-                      psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),&
-                      psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),&
-                      det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),&
-                      tau(i,j,k),rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),&
-                      eps(i,j,k),press(i,j,k),w_lorentz(i,j,k))
-              endif
-            end if
-          
-          end do
-        end do
-      end do
-      
-    end if
-  
-  else
-
-    if (cctk_bbox(5) .ne. 0) then
-
-      do k = sw(3), 1, -1
-        do j = 1, cctk_lsh(2)
-          do i = 1, cctk_lsh(1)
-            
-            posn(1) = x(i,j,k)
-            posn(2) = y(i,j,k)
-            posn(3) = z(i,j,k)
-            
-            if (dot_product(outflowboundary_normal, posn) > 0.d0) then
-              
-              rho(i,j,k)        = rho(i,j,k+1)
-              velx(i,j,k)       = velx(i,j,k+1)
-              vely(i,j,k)       = vely(i,j,k+1)
-              velz(i,j,k)       = velz(i,j,k+1)
-              eps(i,j,k)        = eps(i,j,k+1)
-              press(i,j,k)      = press(i,j,k+1)
-              w_lorentz(i,j,k)  = w_lorentz(i,j,k+1)
-              
-              if(evolve_mhd.ne.0) then
-                 Bvecx(i,j,k)       = Bvecx(i,j,k+1)
-                 Bvecy(i,j,k)       = Bvecy(i,j,k+1)
-                 Bvecz(i,j,k)       = Bvecz(i,j,k+1)
-                 if(clean_divergence.ne.0) then
-                    psidc(i,j,k)    = psidc(i,j,k+1)
-                 endif
-              endif
-
-              psi4pt = 1.0d0
-              det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\
-                   gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))
-              if(evolve_mhd.ne.0) then
-                 call prim2conM(GRHydro_eos_handle,psi4pt*gxx(i,j,k),&
-                      psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),&
-                      psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),&
-                      det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),&
-                      tau(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),&
-                      rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),&
-                      eps(i,j,k),press(i,j,k),w_lorentz(i,j,k))
-              else
-                 call prim2con(GRHydro_eos_handle,psi4pt*gxx(i,j,k),&
-                      psi4pt*gxy(i,j,k),psi4pt*gxz(i,j,k),&
-                      psi4pt*gyy(i,j,k),psi4pt*gyz(i,j,k),psi4pt*gzz(i,j,k),&
-                      det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),&
-                      tau(i,j,k),rho(i,j,k),velx(i,j,k),vely(i,j,k),velz(i,j,k),&
-                      eps(i,j,k),press(i,j,k),w_lorentz(i,j,k))
-              endif
-            end if
-          
-          end do
-        end do
-      end do
-      
-    end if
-
-  end if
-  
-end subroutine GRHydro_OutflowBoundaries
-
-
-