From e3fbc24c806b834739036fb24d50880cac598c1d Mon Sep 17 00:00:00 2001
From: cott <cott@ac85fae7-cede-4708-beff-ae01c7fa1c26>
Date: Wed, 18 Aug 2010 18:42:06 +0000
Subject: * make work with updated EOS_Omni interface

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@113 ac85fae7-cede-4708-beff-ae01c7fa1c26
---
 interface.ccl                     | 6 ++++--
 src/GRHydro_ReadConformalData.F90 | 4 +++-
 src/GRHydro_SimpleWave.F90        | 3 ++-
 3 files changed, 9 insertions(+), 4 deletions(-)

diff --git a/interface.ccl b/interface.ccl
index 8da93f2..70912bd 100644
--- a/interface.ccl
+++ b/interface.ccl
@@ -62,8 +62,9 @@ CCTK_INT FUNCTION EOS_Omni_GetHandle(CCTK_STRING IN name)
 USES FUNCTION EOS_Omni_GetHandle
 
 
-void FUNCTION EOS_Omni_press(CCTK_INT IN eoshandle,         \
+void FUNCTION EOS_Omni_press(CCTK_INT IN eoshandle,      \
                              CCTK_INT IN havetemp,       \
+			     CCTK_REAL IN rf_precision,  \
                              CCTK_INT IN npoints,        \
                              CCTK_REAL IN ARRAY rho,     \
                              CCTK_REAL INOUT ARRAY eps,  \
@@ -75,8 +76,9 @@ void FUNCTION EOS_Omni_press(CCTK_INT IN eoshandle,         \
 
 USES FUNCTION EOS_Omni_press
 
-void FUNCTION EOS_Omni_EpsFromPress(CCTK_INT IN eoshandle,         \
+void FUNCTION EOS_Omni_EpsFromPress(CCTK_INT IN eoshandle,      \
                              CCTK_INT IN havetemp,              \
+			     CCTK_REAL IN rf_precision,         \
                              CCTK_INT IN npoints,               \
                              CCTK_REAL IN ARRAY rho,            \
                              CCTK_REAL INOUT ARRAY eps,         \
diff --git a/src/GRHydro_ReadConformalData.F90 b/src/GRHydro_ReadConformalData.F90
index aaf6cf3..c492d3c 100644
--- a/src/GRHydro_ReadConformalData.F90
+++ b/src/GRHydro_ReadConformalData.F90
@@ -61,6 +61,7 @@ subroutine GRHydro_ReadConformalData(CCTK_ARGUMENTS)
   integer :: keytemp = 0
   integer :: anyerr = 0
   integer :: keyerr(1) = 0
+  real*8  :: rf_precision = 1.0d-10
   real*8  :: xpress(1) = 0.0d0
   real*8  :: xeps(1) = 0.0d0
   real*8  :: xtemp(1) = 0.0d0
@@ -155,7 +156,8 @@ subroutine GRHydro_ReadConformalData(CCTK_ARGUMENTS)
     do j=1,cctk_lsh(2)
       do k=1,cctk_lsh(3)
 #ifdef USE_EOS_OMNI
-        call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,rho(i,j,k),xeps,xtemp,xye,&
+        call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,rf_precision,&
+             n,rho(i,j,k),xeps,xtemp,xye,&
              press(i,j,k),eps(i,j,k),keyerr,anyerr)
 
 #else
diff --git a/src/GRHydro_SimpleWave.F90 b/src/GRHydro_SimpleWave.F90
index 85e191a..2991ff4 100644
--- a/src/GRHydro_SimpleWave.F90
+++ b/src/GRHydro_SimpleWave.F90
@@ -66,6 +66,7 @@ subroutine GRHydro_SimpleWave(CCTK_ARGUMENTS)
   real*8  :: xeps(1) = 0.0d0
   real*8  :: xtemp(1) = 0.0d0
   real*8  :: xye(1) = 0.0d0
+  real*8  :: rf_precision = 1.0d-10
 ! end EOS Omni vars
 #endif
   
@@ -150,7 +151,7 @@ subroutine GRHydro_SimpleWave(CCTK_ARGUMENTS)
       w_lorentz(i,1,1) = 1.d0
 #if USE_EOS_OMNI
       xeps = 1.0d0
-  call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,n,&
+  call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,rf_precision,n,&
        rho(i,1,1),xeps,xtemp,xye,press(i,1,1),keyerr,anyerr)
 #else
       press(i,1,1) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,1,1), 1.0d0)  
-- 
cgit v1.2.3