[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 583)

[rhaas at tapir.caltech.edu]

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)