RIT MHD development:

 Add #include "cctk_Functions.h" and DECLARE_CCTK_FUNCTIONS 
to GRHydro_C2P2CM.F90.
 
 Solve several issues of rank mismatch in the EOS Omni 
interface.

 Initialize the velocities at GRHydro_C2P2CM.F90 since
GRHydro_Con2PrimM_pt.c uses past values as initial guess
for its Newton-Raphson procedure. These unitialized values
were causing C2P2CM to fail from time to time, since it
would access memory positions with very large numerical
values, i.e. v>>1. C2P2CM test always passes now.



git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@118 ac85fae7-cede-4708-beff-ae01c7fa1c26

1 files changed, 13 insertions, 13 deletions

diff --git a/src/GRHydro_ShockTubeM.F90 b/src/GRHydro_ShockTubeM.F90
index 96b0e7c..186cf69 100644
--- a/src/GRHydro_ShockTubeM.F90
+++ b/src/GRHydro_ShockTubeM.F90
@@ -118,22 +118,22 @@ subroutine GRHydro_shocktubeM(CCTK_ARGUMENTS)
           epsl = 1500.d0
           epsr = 1.5d-2
        else if (CCTK_EQUALS(shock_case,"Balsara1")) then
-          rhol = 1.d0
+          rhol = 1.0d0
           rhor = 0.125d0
-          velxl = 0.d0
-          velxr = 0.d0
-          velyl = 0.d0
-          velyr = 0.d0
-          velzl = 0.d0
-          velzr = 0.d0
+          velxl = 0.0d0
+          velxr = 0.0d0
+          velyl = 0.0d0
+          velyr = 0.0d0
+          velzl = 0.0d0
+          velzr = 0.0d0
           bvcxl=0.5d0
           bvcxr=0.5d0
-          bvcyl=1.d0
-          bvcyr=-1.d0
-          bvczl=0.d0
-          bvczr=0.d0
-          epsl = 1.0/rhol
-          epsr = 0.1/rhor
+          bvcyl=1.0d0
+          bvcyr=-1.0d0
+          bvczl=0.0d0
+          bvczr=0.0d0
+          epsl = 1.0d0/rhol
+          epsr = 0.1d0/rhor
        else if (CCTK_EQUALS(shock_case,"Balsara2")) then
           rhol = 1.d0
           rhor = 1.d0