User: rhaas Date: 2014/04/15 02:50 PM

Modified: /trunk/src/ GRHydro_Prim2ConM.F90

Log: GRHydro: no longer compute speed of sound twice in F90 MHD code

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_Prim2ConM.F90 Delta lines: +0 -11 =================================================================== --- trunk/src/GRHydro_Prim2ConM.F90 2014-04-15 19:50:08 UTC (rev 640) +++ trunk/src/GRHydro_Prim2ConM.F90 2014-04-15 19:50:11 UTC (rev 641) @@ -85,17 +85,6 @@ 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