diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-10-14 08:35:29 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-10-14 08:35:29 +0000 |
commit | 17788d73436659db4c7731458edfc5c73b1ee4f1 (patch) | |
tree | fafc9ca154112702f28b34764e3eaf04224e8665 /src/GRHydro_EoSChangeGamma.F90 | |
parent | 7a689eb93032bc54ddd3d7935188baef49b0f759 (diff) |
correct formula used to compute eps in EoSChangeK
this only shows up if we use a general eos (not a polytype) during evolution
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@282 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_EoSChangeGamma.F90')
-rw-r--r-- | src/GRHydro_EoSChangeGamma.F90 | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/src/GRHydro_EoSChangeGamma.F90 b/src/GRHydro_EoSChangeGamma.F90 index 432b469..2820c1b 100644 --- a/src/GRHydro_EoSChangeGamma.F90 +++ b/src/GRHydro_EoSChangeGamma.F90 @@ -195,7 +195,7 @@ subroutine GRHydro_EoSChangeK(CCTK_ARGUMENTS) end if press = local_k * rho**local_gamma - eps = (local_gamma - 1.d0) * local_k * rho**local_gamma + eps = local_k / (local_gamma - 1.d0) * rho**(local_gamma-1) call Primitive2ConservativePolyCells(CCTK_ARGUMENTS) |