* modify reconstruction:

  * get rid of nonsensical re-computation of the conserved tracer
  * remove support for shift_state = 0 and get rid of unnecessary 
    assignments.




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

1 files changed, 6 insertions, 8 deletions

diff --git a/src/GRHydro_Reconstruct.F90 b/src/GRHydro_Reconstruct.F90
index ca93a4c..09cbf16 100644
--- a/src/GRHydro_Reconstruct.F90
+++ b/src/GRHydro_Reconstruct.F90
@@ -1,7 +1,7 @@
  /*@@
    @file      GRHydro_Reconstruct.F90
    @date      Sat Jan 26 02:13:25 2002
-   @author    
+   @author    Bruno Mundim, Josh Faber, Christian D. Ott
    @desc 
    Wrapper routine to perform the reconstruction.
    @enddesc 
@@ -17,7 +17,7 @@
  /*@@
    @routine    Reconstruction
    @date       Sat Jan 26 02:13:47 2002
-   @author     Luca Baiotti, Ian Hawke
+   @author     Luca Baiotti, Ian Hawke, Christian D. Ott
    @desc 
    A wrapper routine to do reconstruction. Currently just does
    TVD on the primitive variables.
@@ -41,6 +41,10 @@ subroutine Reconstruction(CCTK_ARGUMENTS)
   CCTK_INT  :: i,j,k
   CCTK_REAL :: local_min_tracer
 
+  ! set things up
+  if (shift_state .eq. 0) then
+     call CCTK_WARN(0,"This code no longer supports shift_state = 0");
+  endif
 
   if (CCTK_EQUALS(recon_method,"tvd")) then
     ! this handles MHD and non-MHD
@@ -116,12 +120,6 @@ subroutine Reconstruction(CCTK_ARGUMENTS)
     call CCTK_WARN(0,"Variable type to reconstruct not recognized.")
   end if
 
- if(evolve_tracer.ne.0) then
-    ! Call the conserved tracer routine in any case because (accord. to
-    ! Christian Ott) this is the only way this works
-    call Prim2ConservativeTracer(CCTK_PASS_FTOF)
- end if
-
   return
 
 end subroutine Reconstruction