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, \