diff options
Diffstat (limited to 'src/GRHydro_PoloidalMagFieldM.F90')
-rw-r--r-- | src/GRHydro_PoloidalMagFieldM.F90 | 10 |
1 files changed, 9 insertions, 1 deletions
diff --git a/src/GRHydro_PoloidalMagFieldM.F90 b/src/GRHydro_PoloidalMagFieldM.F90 index 37b4377..c809001 100644 --- a/src/GRHydro_PoloidalMagFieldM.F90 +++ b/src/GRHydro_PoloidalMagFieldM.F90 @@ -85,7 +85,8 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) write(*,*)'GRHydro_InitData: Setting up initial poloidal magnetic field' - +! Initialize to zero + Bvec = 0.0d0 do i=2,nx-1 do j=2,ny-1 @@ -106,6 +107,7 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) !! write(*,*)'delPcut = ', delPcut !! write(*,*)'maxP_Pcut = ', maxP_Pcut +#warning "This algorithm does only work on Cartesian grids!!" rho_dx = 0.5d0*(rho(i+1,j,k)-rho(i-1,j,k))/dx rho_dy = 0.5d0*(rho(i,j+1,k)-rho(i,j-1,k))/dy rho_dz = 0.5d0*(rho(i,j,k+1)-rho(i,j,k-1))/dz @@ -183,6 +185,12 @@ subroutine GRHydro_PoloidalMagFieldM(CCTK_ARGUMENTS) if (clean_divergence .ne. 0) then psidcrhs =0.0 endif + + !Bvec = 0 + !lBvec = 0 + !Bcons = 0 + !Avec = 0 + return end subroutine GRHydro_PoloidalMagFieldM |