diff options
Diffstat (limited to 'src/GRHydro_PoloidalMagFieldM.F90')
| -rw-r--r-- | src/GRHydro_PoloidalMagFieldM.F90 | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/src/GRHydro_PoloidalMagFieldM.F90 b/src/GRHydro_PoloidalMagFieldM.F90 index df74d75..47ce06c 100644 --- a/src/GRHydro_PoloidalMagFieldM.F90 +++ b/src/GRHydro_PoloidalMagFieldM.F90 @@ -58,6 +58,7 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) CCTK_INT :: i, j, k, nx, ny, nz CCTK_REAL :: det + CCTK_REAL :: sdet CCTK_REAL :: dx,dy,dz CCTK_REAL :: rhofac, delPcut, maxP_Pcut CCTK_REAL :: Aphi, Ax, Ay, Az @@ -112,12 +113,14 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) Ay_dx = Aphi + x(i,j,k)*Aphi_dx Ay_dz = x(i,j,k)*Aphi_dz - Bvecx(i,j,k) = alp(i,j,k)*Ay_dz - Bvecy(i,j,k) = -alp(i,j,k)*Ax_dz - Bvecz(i,j,k) = alp(i,j,k)*(Ax_dy-Ay_dx) 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)) - + sdet = sqrt(det) + + Bvecx(i,j,k) = Ay_dz/sdet + Bvecy(i,j,k) = - Ax_dz/sdet + Bvecz(i,j,k) = (Ax_dy-Ay_dx)/sdet + if (CCTK_EQUALS(GRHydro_eos_type,"Polytype")) then call Prim2ConPolyM(GRHydro_eos_handle,gxx(i,j,k),gxy(i,j,k),& gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k),& |