diff options
Diffstat (limited to 'src/GRHydro_PoloidalMagFieldM.F90')
-rw-r--r-- | src/GRHydro_PoloidalMagFieldM.F90 | 18 |
1 files changed, 17 insertions, 1 deletions
diff --git a/src/GRHydro_PoloidalMagFieldM.F90 b/src/GRHydro_PoloidalMagFieldM.F90 index 52e9ca8..37b4377 100644 --- a/src/GRHydro_PoloidalMagFieldM.F90 +++ b/src/GRHydro_PoloidalMagFieldM.F90 @@ -29,6 +29,9 @@ Keisuke Taniguchi - Phys. Rev. D 78, 024012 (2008) #define Bconsx(i,j,k) Bcons(i,j,k,1) #define Bconsy(i,j,k) Bcons(i,j,k,2) #define Bconsz(i,j,k) Bcons(i,j,k,3) +#define Avecx(i,j,k) Avec(i,j,k,1) +#define Avecy(i,j,k) Avec(i,j,k,2) +#define Avecz(i,j,k) Avec(i,j,k,3) /*@@ @routine GRHydro_PoloidalMagFieldM @@ -56,7 +59,7 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS - CCTK_INT :: i, j, k, nx, ny, nz + CCTK_INT :: i, j, k, nx, ny, nz, set_Avec CCTK_REAL :: det CCTK_REAL :: sdet CCTK_REAL :: dx,dy,dz @@ -75,8 +78,15 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) dy = CCTK_DELTA_SPACE(2) dz = CCTK_DELTA_SPACE(3) + set_Avec = 0 + if ( CCTK_EQUALS(Bvec_evolution_method,"GRHydro_Avec") ) then + set_Avec = 1 + end if + write(*,*)'GRHydro_InitData: Setting up initial poloidal magnetic field' + + do i=2,nx-1 do j=2,ny-1 do k=2,nz-1 @@ -123,6 +133,12 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) Bvecy(i,j,k) = - Ax_dz/sdet Bvecz(i,j,k) = (Ax_dy-Ay_dx)/sdet + if ( set_Avec.gt.0 ) then + Avecx(i,j,k) = Ax + Avecy(i,j,k) = Ay + Avecz(i,j,k) = Az + end if + !Bvecx(i,j,k) = 0.0d0 !Bvecy(i,j,k) = 0.0d0 !Bvecz(i,j,k) = 0.00000001/sdet |