diff options
author | cott <cott@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2010-08-18 18:42:06 +0000 |
---|---|---|
committer | cott <cott@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2010-08-18 18:42:06 +0000 |
commit | e3fbc24c806b834739036fb24d50880cac598c1d (patch) | |
tree | e85e96720b9495032b58cdf2c429d2f0a7b83f7d | |
parent | 8fc217e9d9b586969b3b2379ae40cb1741242983 (diff) |
* 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
-rw-r--r-- | interface.ccl | 6 | ||||
-rw-r--r-- | src/GRHydro_ReadConformalData.F90 | 4 | ||||
-rw-r--r-- | 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) |