aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2014-04-15 19:50:11 +0000
committerrhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45>2014-04-15 19:50:11 +0000
commitc493c3fd7b27c831678b67077d4b5a830f4f5ee0 (patch)
tree52b6462854db712a32f4582b5a28d822864ff96a
parent2e0ea99d145009d33b397bee59e861a6ccd3425b (diff)
GRHydro: no longer compute speed of sound twice in F90 MHD code
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@641 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r--src/GRHydro_Prim2ConM.F9011
1 files changed, 0 insertions, 11 deletions
diff --git a/src/GRHydro_Prim2ConM.F90 b/src/GRHydro_Prim2ConM.F90
index e330ef0..978fede 100644
--- a/src/GRHydro_Prim2ConM.F90
+++ b/src/GRHydro_Prim2ConM.F90
@@ -85,17 +85,6 @@ subroutine primitive2conservativeM(CCTK_ARGUMENTS)
pBprim = loc(Bvec)
end if
- ! currently only the C++ code uses the speed-of-sound grid function, however
- ! in the MHD code, both C++ and F90 code share the same prim2con routine (this
- ! one). GRHydro_SpeedOfSound is a misnomer since it computes both cs2 and
- ! pressure so we call it first so that the actual prim2con below overwrites
- ! the pressur again. This is certainly not the "right" way of doing this but
- ! it lets us have a working version for now.
- ! TODO: extend prim2con in GRHydro_Prim2Con.c to handle MHD.
- if (use_cxx_code .ne. 0) then
- call GRHydro_SpeedOfSound(cctkGH)
- end if
-
! constraint transport needs to be able to average fluxes in the directions
! other that flux_direction, which in turn need the primitives on interfaces