diff options
Diffstat (limited to 'src/GRHydro_SourceM.F90')
-rw-r--r-- | src/GRHydro_SourceM.F90 | 10 |
1 files changed, 4 insertions, 6 deletions
diff --git a/src/GRHydro_SourceM.F90 b/src/GRHydro_SourceM.F90 index cd734eb..4b1321c 100644 --- a/src/GRHydro_SourceM.F90 +++ b/src/GRHydro_SourceM.F90 @@ -76,7 +76,7 @@ subroutine SourceTermsM(CCTK_ARGUMENTS) CCTK_INT :: i, j, k, nx, ny, nz CCTK_REAL :: one, two, half CCTK_REAL :: t00, t0x, t0y, t0z, txx, txy, txz, tyy, tyz, tzz - CCTK_REAL :: sqrtdet, det, uxx, uxy, uxz, uyy, uyz, uzz + CCTK_REAL :: sqrtdet, uxx, uxy, uxz, uyy, uyz, uzz CCTK_REAL :: shiftx, shifty, shiftz, velxshift, velyshift, velzshift CCTK_REAL :: vlowx, vlowy, vlowz CCTK_REAL :: dx_betax, dx_betay, dx_betaz, dy_betax, dy_betay,& @@ -254,7 +254,7 @@ subroutine SourceTermsM(CCTK_ARGUMENTS) !$OMP PARALLEL DO PRIVATE(i, j, k, local_spatial_order,& !$OMP localgxx,localgxy,localgxz,localgyy,localgyz,localgzz,& - !$OMP det,sqrtdet,shiftx,shifty,shiftz,& + !$OMP sqrtdet,shiftx,shifty,shiftz,& !$OMP dx_betax,dx_betay,dx_betaz,dy_betax,dy_betay,dy_betaz,& !$OMP dz_betax,dz_betay,dz_betaz,velxshift,velyshift,velzshift,& !$OMP vlowx,vlowy,vlowz,Bvecxlow,Bvecylow,Bveczlow, & @@ -290,10 +290,8 @@ subroutine SourceTermsM(CCTK_ARGUMENTS) localgyz = g23(i,j,k) localgzz = g33(i,j,k) - det = SPATIAL_DETERMINANT(localgxx, localgxy, localgxz,\ - localgyy, localgyz, localgzz) - sqrtdet = sqrt(det) - call UpperMetric(uxx, uxy, uxz, uyy, uyz, uzz, det, localgxx,& + sqrtdet = sdetg(i,j,k) + call UpperMetric(uxx, uxy, uxz, uyy, uyz, uzz, sqrtdet*sqrtdet, localgxx,& localgxy, localgxz, localgyy, localgyz, localgzz) |