diff options
Diffstat (limited to 'src/GRHydro_BvecfromAvec.F90')
-rw-r--r-- | src/GRHydro_BvecfromAvec.F90 | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/GRHydro_BvecfromAvec.F90 b/src/GRHydro_BvecfromAvec.F90 index ea0a58d..30da906 100644 --- a/src/GRHydro_BvecfromAvec.F90 +++ b/src/GRHydro_BvecfromAvec.F90 @@ -53,7 +53,7 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) implicit none CCTK_REAL :: Az_y, Ay_z, Ax_z, Az_x, Ay_x, Ax_y - CCTK_REAL :: det, isdet + CCTK_REAL :: sdet, isdet CCTK_REAL :: dx, dy, dz, idx, idy, idz CCTK_INT :: i, j, k, nx, ny, nz, GRHydro_UseGeneralCoordinates @@ -102,7 +102,7 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) call CCTK_WARN(1,"Bvec from Avec start."); - !$OMP PARALLEL DO PRIVATE( det, isdet, local_spatial_order, & + !$OMP PARALLEL DO PRIVATE( sdet, isdet, local_spatial_order, & !$OMP Az_y, Ay_z, Ax_z, Az_x, Ay_x, Ax_y) do k = GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil+1 do j = GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil+1 @@ -129,8 +129,8 @@ subroutine GRHydro_BvecfromAvec(CCTK_ARGUMENTS) Ax_y = DIFF_Y_4(Avecx) end if - det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) - isdet = 1.d0/sqrt(det) + sdet = sdetg(i,j,k) + isdet = 1.d0/sdet Bvecx(i,j,k) = isdet*( Az_y - Ay_z ) Bvecy(i,j,k) = isdet*( Ax_z - Az_x ) |