diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:13 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:13 +0000 |
| commit | 7c39d4964e8b048681a402347b2dc7c4d57f8fee (patch) | |
| tree | 5b2958e2f36fdc0366e27dcd56c0e05e3ff1340d /src | |
| parent | c5d0dacc66b08443898399c0eb14ac973bbce544 (diff) | |
GRHydro_InitData: Add basic vector potential support
Initial Avec constrained to poloidal at the moment.
From: Tanja Bode <tanja.bode@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@195 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src')
| -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 |