User: rhaas Date: 2013/10/29 10:19 AM
Modified: /trunk/src/ GRHydro_Wrappers.F90
Log: GRHydro: correct usage of sqrt(detg) in prim2con wrappers
found by Ian Hawke
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Wrappers.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_Wrappers.F90 2013-10-28 16:36:15 UTC (rev 582) +++ trunk/src/GRHydro_Wrappers.F90 2013-10-29 15:19:01 UTC (rev 583) @@ -193,7 +193,7 @@
sdetg = sqrt(det)
- call prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, dens,& + call prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg, dens,& sx, sy, sz, tau, rho, velx, vely, velz, epsilon,& press, w_lorentz) end subroutine @@ -219,7 +219,7 @@
sdetg = sqrt(det)
- call prim2conpolytype(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, dens,& + call prim2conpolytype(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg, dens,& sx, sy, sz, tau, rho, velx, vely, velz, epsilon,& press, w_lorentz) end subroutine @@ -318,7 +318,7 @@
sdetg = sqrt(det)
- call prim2conpolytypeM(handle, gxx, gxy, gxz, gyy, gyz, gzz, det,& + call prim2conpolytypeM(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg,& dens, sx, sy, sz, tau, Bconsx, Bconsy, Bconsz,& rho, velx, vely, velz, epsilon, press, Bvecx,& Bvecy, Bvecz, w_lorentz)