GRHydro: Implement entropy evolution (separated from temp evolution)

From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>

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

1 files changed, 127 insertions, 89 deletions

diff --git a/src/GRHydro_Boundaries.F90 b/src/GRHydro_Boundaries.F90
index c04994c..e9def4a 100644
--- a/src/GRHydro_Boundaries.F90
+++ b/src/GRHydro_Boundaries.F90
@@ -85,6 +85,11 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS)
       call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
     !endif
   endif
+
+  if(evolve_entropy.ne.0) then
+      call SetCartSymGN(ierr, cctkGH, sym, "HydroBase::entropy")
+      call SetCartSymGN(ierr, cctkGH, sym, "GRHydro::entropycons")
+  endif
   
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::atmosphere_mask_real")
   
@@ -94,9 +99,9 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS)
   
   !if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
+        call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[0]")
       else
-	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
+        call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[0]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[0]")
@@ -113,15 +118,15 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS)
   
   !if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
+         call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[1]")
       else
-	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
+         call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[1]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[1]")
   if(evolve_mhd.ne.0) then
      !if (sync_conserved_only .eq. 0) then
-	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
+     call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::Bvec[1]")
      !endif
      call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::Bcons[1]")
   endif
@@ -132,9 +137,9 @@ subroutine GRHydro_InitSymBound(CCTK_ARGUMENTS)
   
   !if (sync_conserved_only .eq. 0) then
       if (general_coordinates .ne. 0) then
-	 call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
+        call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::lvel[2]")
       else
-	 call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
+        call SetCartSymVN(ierr, cctkGH, sym, "HydroBase::vel[2]")
       endif
   !endif
   call SetCartSymVN(ierr, cctkGH, sym, "GRHydro::scon[2]")
@@ -201,31 +206,31 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
          "GRHydro::scon", "Flat")
     if (sync_conserved_only .eq. 0) then
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::w_lorentz", "Flat")
+           "HydroBase::w_lorentz", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::rho", "Flat")
+           "HydroBase::rho", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::press", "Flat")
+           "HydroBase::press", "Flat")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::eps", "Flat")
+           "HydroBase::eps", "Flat")
       if (general_coordinates .ne. 0) then
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "GRHydro::lvel", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lvel", "Flat")
       else
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "HydroBase::vel", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::vel", "Flat")
       endif
     endif
     if(evolve_mhd.ne.0) then
-       if (sync_conserved_only .eq. 0) then
-	 if (general_coordinates .ne. 0) then
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "GRHydro::lBvec", "Flat")
-	 else
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "HydroBase::Bvec", "Flat")
-	 endif
-       endif
+      if (sync_conserved_only .eq. 0) then
+        if (general_coordinates .ne. 0) then
+          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "GRHydro::lBvec", "Flat")
+        else
+          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "HydroBase::Bvec", "Flat")
+        endif
+      endif
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Bcons", "Flat")
        if(clean_divergence.ne.0) then
@@ -236,8 +241,8 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 
     if(evolve_tracer.ne.0) then 
       if (sync_conserved_only .eq. 0) then
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "GRHydro::GRHydro_tracers", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::GRHydro_tracers", "Flat")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::GRHydro_cons_tracers", "Flat")
@@ -245,8 +250,8 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 
     if(evolve_y_e.ne.0) then
       if (sync_conserved_only .eq. 0) then
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "HydroBase::Y_e", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Y_e", "Flat")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Y_e_con", "Flat")
@@ -254,13 +259,21 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 
     if(evolve_temper.ne.0) then
       if (sync_conserved_only .eq. 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::temperature", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "Flat")
       endif
     endif
 
+    if(evolve_entropy.ne.0) then
+      if (sync_conserved_only .eq. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "Flat")
+      endif
+      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+           "GRHydro::entropycons", "Flat")
+    endif
 
   endif
 
@@ -274,44 +287,44 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
     
     if (sync_conserved_only .eq. 0) then
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::w_lorentz", "None")
+           "HydroBase::w_lorentz", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::rho", "None")
+           "HydroBase::rho", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::press", "None")
+           "HydroBase::press", "None")
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	    "HydroBase::eps", "None")
+           "HydroBase::eps", "None")
       if (general_coordinates .ne. 0) then
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "GRHydro::lvel", "None")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lvel", "None")
       else
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "HydroBase::vel", "None")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::vel", "None")
       endif
     endif
     
     if(evolve_mhd.ne.0) then
-       if (sync_conserved_only .eq. 0) then
-	 if (general_coordinates .ne. 0) then
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "GRHydro::lBvec", "None")
-	 else
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "HydroBase::Bvec", "None")
-	 endif
-       endif
-       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-            "GRHydro::Bcons", "None")
-       if(clean_divergence.ne.0) then
+      if (sync_conserved_only .eq. 0) then
+        if (general_coordinates .ne. 0) then
           ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "GRHydro::psidc", "None")
-       endif
+               "GRHydro::lBvec", "None")
+        else
+          ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+               "HydroBase::Bvec", "None")
+        endif
+      endif
+      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+           "GRHydro::Bcons", "None")
+      if(clean_divergence.ne.0) then
+         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+              "GRHydro::psidc", "None")
+      endif
     endif
 
     if(evolve_tracer.ne.0) then 
       if (sync_conserved_only .eq. 0) then 
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "GRHydro::GRHydro_tracers", "None")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::GRHydro_tracers", "None")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::GRHydro_cons_tracers", "None")
@@ -319,8 +332,8 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 
     if(evolve_y_e.ne.0) then
       if (sync_conserved_only .eq. 0) then
-	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "HydroBase::Y_e", "None")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Y_e", "None")
       endif
       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
             "GRHydro::Y_e_con", "None")
@@ -328,11 +341,20 @@ subroutine GRHydro_Boundaries(CCTK_ARGUMENTS)
 
     if(evolve_temper.ne.0) then
       if (sync_conserved_only .eq. 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::temperature", "None")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "None")
+      endif
+    endif
+
+    if(evolve_entropy.ne.0) then
+      if (sync_conserved_only .eq. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "None")
       endif
+      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+           "GRHydro::entropycons", "None")
     endif
 
   endif
@@ -394,11 +416,11 @@ subroutine GRHydro_SelectPrimitiveInitialGuessesBoundaries(CCTK_ARGUMENTS)
     endif
     if(evolve_mhd.ne.0) then
       if (general_coordinates .ne. 0) then
-         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "GRHydro::lBvec", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lBvec", "Flat")
       else
-         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "HydroBase::Bvec", "Flat")
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Bvec", "Flat")
       endif
     endif
 
@@ -419,6 +441,10 @@ subroutine GRHydro_SelectPrimitiveInitialGuessesBoundaries(CCTK_ARGUMENTS)
        !      "HydroBase::entropy", "Flat")
     endif
 
+    if(evolve_entropy.ne.0) then
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "Flat")
+    endif
 
   endif
 
@@ -442,12 +468,11 @@ subroutine GRHydro_SelectPrimitiveInitialGuessesBoundaries(CCTK_ARGUMENTS)
     if(evolve_mhd.ne.0) then
       if (general_coordinates .ne. 0) then
         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-	       "GRHydro::lBvec", "None")
+             "GRHydro::lBvec", "None")
       else
         ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-               "HydroBase::Bvec", "None")
+             "HydroBase::Bvec", "None")
       endif
-       
     endif
 
     !if(evolve_tracer.ne.0) then 
@@ -469,6 +494,11 @@ subroutine GRHydro_SelectPrimitiveInitialGuessesBoundaries(CCTK_ARGUMENTS)
        !      "HydroBase::entropy", "None")
     endif
 
+    if(evolve_entropy.ne.0) then
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "None")
+    endif
+
   end if
 
   if (CCTK_EQUALS(bound,"scalar")) then
@@ -528,24 +558,25 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
              "HydroBase::vel", "Flat")
     endif
+
     if(evolve_mhd.ne.0) then
-	 if (general_coordinates .ne. 0) then
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "GRHydro::lBvec", "Flat")
-	 else
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "HydroBase::Bvec", "Flat")
-	 endif
+      if (general_coordinates .ne. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lBvec", "Flat")
+      else
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Bvec", "Flat")
+      endif
     endif
 
     if(evolve_tracer.ne.0) then 
-    	 ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-    	       "GRHydro::GRHydro_tracers", "Flat")
+      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")
+      ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+           "HydroBase::Y_e", "Flat")
     endif
 
     if(evolve_temper.ne.0) then
@@ -555,6 +586,10 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
              "HydroBase::entropy", "Flat")
     endif
 
+    if(evolve_entropy.ne.0) then
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "Flat")
+    endif
 
   endif
 
@@ -576,20 +611,18 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
     endif
     
     if(evolve_mhd.ne.0) then
-	 if (general_coordinates .ne. 0) then
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "GRHydro::lBvec", "None")
-	 else
-	    ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
-		  "HydroBase::Bvec", "None")
-	 endif
-       
+      if (general_coordinates .ne. 0) then
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "GRHydro::lBvec", "None")
+      else
+        ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::Bvec", "None")
+      endif
     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
@@ -605,6 +638,11 @@ subroutine GRHydro_SelectPrimitiveBoundaries(CCTK_ARGUMENTS)
              "HydroBase::entropy", "None")
     endif
 
+    if(evolve_entropy.ne.0) then
+       ierr = ierr + Boundary_SelectGroupForBC(cctkGH, faces, GRHydro_stencil, -ione, &
+             "HydroBase::entropy", "None")
+    endif
+
   end if
 
   if (CCTK_EQUALS(bound,"scalar")) then