diff options
Diffstat (limited to 'src/GRHydro_C2P2C.F90')
| -rw-r--r-- | src/GRHydro_C2P2C.F90 | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/src/GRHydro_C2P2C.F90 b/src/GRHydro_C2P2C.F90 index 5808c1b..a2523a8 100644 --- a/src/GRHydro_C2P2C.F90 +++ b/src/GRHydro_C2P2C.F90 @@ -34,10 +34,12 @@ subroutine c2p2c(CCTK_ARGUMENTS) DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS +#if !USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" +#endif integer didit,i,j,k,nx,ny,nz CCTK_REAL det @@ -50,6 +52,21 @@ subroutine c2p2c(CCTK_ARGUMENTS) CCTK_INT C2P_failed logical epsnegative +#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 + call CCTK_WARN(1,"This test works only with Ideal_Fluid EoS") nx = cctk_lsh(1) @@ -90,8 +107,16 @@ subroutine c2p2c(CCTK_ARGUMENTS) epsnegative = .false. GRHydro_rho_min = 1.e-10 +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + GRHydro_rho_min,1.0d0,xtemp,xye,pmin,keyerr,anyerr) + + call EOS_Omni_EpsFromPress(poly_eoskey,keytemp,n,& + GRHydro_rho_min,xeps,xtemp,xye,pmin,epsmin,keyerr,anyerr) +#else pmin = EOS_Pressure(GRHydro_eos_handle, GRHydro_rho_min, 1.0d0) epsmin = EOS_SpecificIntEnergy(GRHydro_eos_handle, GRHydro_rho_min, pmin) +#endif C2P_failed = 0 write(*,*) 'C2P2C test: initial values.' |