diff options
| author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-11-09 01:53:00 +0000 |
|---|---|---|
| committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-11-09 01:53:00 +0000 |
| commit | 373bf59d5c9a3744030ab72825cd3bd75f3f4aba (patch) | |
| tree | 6a75d6767e5be78f89c8920869d932bf9a50e3af /src/GRHydro_EoSChangeGamma.F90 | |
| parent | 264d86a805094d26866db74aae895c5dbb11e421 (diff) | |
GRHydro: remove usage of EOS_Module, poly_gamma_ini
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@431 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_EoSChangeGamma.F90')
| -rw-r--r-- | src/GRHydro_EoSChangeGamma.F90 | 20 |
1 files changed, 8 insertions, 12 deletions
diff --git a/src/GRHydro_EoSChangeGamma.F90 b/src/GRHydro_EoSChangeGamma.F90 index fb075e2..52c2d3d 100644 --- a/src/GRHydro_EoSChangeGamma.F90 +++ b/src/GRHydro_EoSChangeGamma.F90 @@ -39,8 +39,6 @@ subroutine GRHydro_EoSChangeGamma(CCTK_ARGUMENTS) - USE EOS_Omni_Module, only: press_gf, inv_rho_gf, poly_k_cgs, rho_gf - implicit none ! save memory when MP is not used @@ -95,8 +93,8 @@ subroutine GRHydro_EoSChangeGamma(CCTK_ARGUMENTS) end if local_Gamma = 1.0d0 + xpress/xeps - press = press_gf * poly_k_cgs * & - (rho * inv_rho_gf)**local_Gamma + press = poly_k * & + rho**local_Gamma eps = press / (rho * (local_Gamma - 1.d0)) !!$ Change the pressure and specific internal energy @@ -277,8 +275,6 @@ end subroutine GRHydro_EoSChangeK subroutine GRHydro_EoSChangeGammaK_Shibata(CCTK_ARGUMENTS) - USE EOS_Omni_Module, only: press_gf, inv_rho_gf, poly_k_cgs, rho_gf - implicit none ! save memory when MP is not used @@ -293,7 +289,7 @@ subroutine GRHydro_EoSChangeGammaK_Shibata(CCTK_ARGUMENTS) CCTK_INT :: i, j, k CCTK_REAL :: det - CCTK_REAL :: local_Gamma, local_k, eos_k_initial_cgs + CCTK_REAL :: local_Gamma, local_k, eos_k_initial CCTK_REAL, dimension(cctk_lsh(1),cctk_lsh(2),cctk_lsh(3)) :: Q @@ -340,13 +336,13 @@ subroutine GRHydro_EoSChangeGammaK_Shibata(CCTK_ARGUMENTS) local_Gamma = 1.0d0 + xpress/xeps local_K = xpress - eos_k_initial_cgs = initial_k * rho_gf**initial_Gamma / press_gf + eos_k_initial = initial_k - press = (local_Gamma - 1.d0) / (initial_Gamma - 1.0d0 ) * press_gf * eos_k_initial_cgs * & - (rho * rho_gf) ** initial_Gamma + press = (local_Gamma - 1.d0) / (initial_Gamma - 1.0d0 ) * eos_k_initial * & + rho ** initial_Gamma - eps = press_gf * eos_k_initial_cgs * & - (rho * inv_rho_gf) ** initial_Gamma / & + eps = eos_k_initial * & + rho ** initial_Gamma / & (rho * (initial_Gamma - 1.0d0)) do k = 1, cctk_lsh(3) |