[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 113)

[cott at tapir.caltech.edu]

User: cott
Date: 2010/08/18 01:42 PM

Modified:
 /trunk/
  interface.ccl
 /trunk/src/
  GRHydro_ReadConformalData.F90, GRHydro_SimpleWave.F90

Log:
 * make work with updated EOS_Omni interface

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_ReadConformalData.F90
Delta lines: +3 -1
===================================================================
--- trunk/src/GRHydro_ReadConformalData.F90	2010-08-18 18:15:02 UTC (rev 112)
+++ trunk/src/GRHydro_ReadConformalData.F90	2010-08-18 18:42:06 UTC (rev 113)
@@ -61,6 +61,7 @@
   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 @@
     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

File [modified]: GRHydro_SimpleWave.F90
Delta lines: +2 -1
===================================================================
--- trunk/src/GRHydro_SimpleWave.F90	2010-08-18 18:15:02 UTC (rev 112)
+++ trunk/src/GRHydro_SimpleWave.F90	2010-08-18 18:42:06 UTC (rev 113)
@@ -66,6 +66,7 @@
   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 @@
       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)  

Directory: /trunk/
==================

File [modified]: interface.ccl
Delta lines: +4 -2
===================================================================
--- trunk/interface.ccl	2010-08-18 18:15:02 UTC (rev 112)
+++ trunk/interface.ccl	2010-08-18 18:42:06 UTC (rev 113)
@@ -62,8 +62,9 @@
 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 @@
 
 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,         \