diff options
Diffstat (limited to 'src/GRHydro_SimpleWave.F90')
| -rw-r--r-- | src/GRHydro_SimpleWave.F90 | 24 |
1 files changed, 23 insertions, 1 deletions
diff --git a/src/GRHydro_SimpleWave.F90 b/src/GRHydro_SimpleWave.F90 index 829170f..d1f10d4 100644 --- a/src/GRHydro_SimpleWave.F90 +++ b/src/GRHydro_SimpleWave.F90 @@ -46,14 +46,31 @@ subroutine GRHydro_SimpleWave(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS +#ifndef USE_EOS_OMNI #ifdef _EOS_BASE_INC_ #undef _EOS_BASE_INC_ #endif #include "EOS_Base.inc" - +#endif + CCTK_INT :: i, j, k, nx, ny, nz CCTK_REAL :: dr, k1, k2, k3, k4, in_data, old_data, source_data, new_data, c_0, det, pi +#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_INFO("Setting initial data for a simple wave as Anile Miller Motta") call CCTK_WARN(1, "The simple-wave initial-data routine works only for unigrid and on one node.") @@ -133,7 +150,12 @@ subroutine GRHydro_SimpleWave(CCTK_ARGUMENTS) eps(i,1,1) = rho_abs_min**(1.d0/3.d0) velx(i,1,1) = 0.d0 w_lorentz(i,1,1) = 1.d0 +#if USE_EOS_OMNI + call EOS_Omni_press(poly_eoskey,keytemp,n,& + rho(i,1,1),1.0d0,xtemp,xye,press(i,1,1),keyerr,anyerr) +#else press(i,1,1) = EOS_Pressure(GRHydro_polytrope_handle, rho(i,1,1), 1.0d0) +#endif ! polytrope only (initial data) end if |