diff options
Diffstat (limited to 'src/GRHydro_ReadConformalData.F90')
-rw-r--r-- | src/GRHydro_ReadConformalData.F90 | 30 |
1 files changed, 25 insertions, 5 deletions
diff --git a/src/GRHydro_ReadConformalData.F90 b/src/GRHydro_ReadConformalData.F90 index 5829d80..4513f47 100644 --- a/src/GRHydro_ReadConformalData.F90 +++ b/src/GRHydro_ReadConformalData.F90 @@ -51,9 +51,23 @@ subroutine GRHydro_ReadConformalData(CCTK_ARGUMENTS) CCTK_INT :: i,j,k,handle,ierr - CCTK_REAL :: eos_k, eos_gamma, rho_min, det - - +! CCTK_REAL :: eos_k, eos_gamma + CCTK_REAL :: rho_min, det + +#if USE_EOS_OMNI +! begin EOS Omni vars + integer :: n = 1 + integer :: poly_eoskey = 0 + integer :: keytemp = 0 + integer :: anyerr = 0 + integer :: keyerr(1) = 0 + real*8 :: xpress = 0.0d0 + real*8 :: xeps = 0.0d0 + real*8 :: xtemp = 0.0d0 + real*8 :: xye = 0.0d0 + poly_eoskey = GRHydro_poly_eoskey +! end EOS Omni vars +#endif ! only gxx has been read in; copy it into gyy and gzz as well gyy = gxx @@ -131,8 +145,8 @@ subroutine GRHydro_ReadConformalData(CCTK_ARGUMENTS) ! if (handle < 0) call CCTK_WARN(0,"For this hack you need to compile with EOS_Polytrope") - eos_k = EOS_Pressure(GRHydro_eos_handle, 1.0, 1.0) - eos_gamma = 1.0 + eos_k / EOS_SpecificIntEnergy(GRHydro_eos_handle,1.0,1.0) +! eos_k = EOS_Pressure(GRHydro_eos_handle, 1.0, 1.0) +! eos_gamma = 1.0 + eos_k / EOS_SpecificIntEnergy(GRHydro_eos_handle,1.0,1.0) ! press = eos_k * rho**eos_gamma @@ -140,7 +154,13 @@ subroutine GRHydro_ReadConformalData(CCTK_ARGUMENTS) do i=1,cctk_lsh(1) 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,press(i,j,k),eps(i,j,k),keyerr,anyerr) + +#else eps(i,j,k) = EOS_SpecificIntEnergy(GRHydro_eos_handle,rho(i,j,k),press(i,j,k)) +#endif end do end do end do |