From 9186bf04c53973cf99251b967d4b08de4aa85e21 Mon Sep 17 00:00:00 2001
From: rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>
Date: Tue, 5 Jun 2012 20:28:24 +0000
Subject: GRHydro: sync those primitives which are needed for initial Con2Prim
 guesses

Patch by Christian Reisswig.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@341 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
---
 src/GRHydro_Boundaries.F90 | 63 ++++++++++++++++++++++++----------------------
 1 file changed, 33 insertions(+), 30 deletions(-)

(limited to 'src/GRHydro_Boundaries.F90')

diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90
index 2c6b9ce..099e195 100644
--- a/src/GRHydro_Boundaries.F90
+++ b/src/GRHydro_Boundaries.F90
@@ -378,14 +378,17 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
 
 !!$Flat boundaries if required  
 
+! The commented out code are those primitives which do not require and explicit
+! initial guess for Con2Prim. The guesses are either not needed or computed from other quantities!
+
   if (CCTK_EQUALS(bound,"flat")) then
     
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::w_lorentz", "Flat")
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::w_lorentz", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "Flat")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::press", "Flat")
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::press", "Flat")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "Flat")
     if (general_coordinates .ne. 0) then
@@ -405,33 +408,33 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
       endif
     endif
 
-    if(evolve_tracer.ne.0) then 
-      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::GRHydro_tracers", "Flat")
-    endif
+    !if(evolve_tracer.ne.0) then 
+    !  ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !         "GRHydro::GRHydro_tracers", "Flat")
+    !endif
 
-    if(evolve_y_e.ne.0) then
-      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Y_e", "Flat")
-    endif
+    !if(evolve_y_e.ne.0) then
+    !  ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !         "HydroBase::Y_e", "Flat")
+    !endif
 
     if(evolve_temper.ne.0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
              "HydroBase::temperature", "Flat")
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "Flat")
+       !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+       !      "HydroBase::entropy", "Flat")
     endif
 
 
   endif
 
   if (CCTK_EQUALS(bound,"none")) then      
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::w_lorentz", "None")
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::w_lorentz", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::rho", "None")
-    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-          "HydroBase::press", "None")
+    !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !      "HydroBase::press", "None")
     ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
           "HydroBase::eps", "None")
     if (general_coordinates .ne. 0) then
@@ -453,23 +456,23 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
        
     endif
 
-    if(evolve_tracer.ne.0) then 
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "GRHydro::GRHydro_tracers", "None")
-      
-    endif
+    !if(evolve_tracer.ne.0) then 
+    !   ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !         "GRHydro::GRHydro_tracers", "None")
+    !  
+    !endif
 
-    if(evolve_y_e.ne.0) then
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::Y_e", "None")
-      
-    endif
+    !if(evolve_y_e.ne.0) then
+    !   ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+    !         "HydroBase::Y_e", "None")
+    !  
+    !endif
 
     if(evolve_temper.ne.0) then
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
              "HydroBase::temperature", "None")
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-             "HydroBase::entropy", "None")
+       !ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+       !      "HydroBase::entropy", "None")
     endif
 
   end if
-- 
cgit v1.2.3