diff options
| -rw-r--r-- | src/GRHydro_Prim2ConM.F90 | 11 |
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 |