diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-07-06 18:12:15 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-07-06 18:12:15 +0000 |
| commit | acb151ad12413f3f0ef74476f500002f46777d79 (patch) | |
| tree | c3c13c928ab07e88edb0a5add40e444026fa5a06 /src | |
| parent | b8e52dfff0076cd76576ab1eecd8a0d4abdd57bf (diff) | |
GRHydro_InitData: Reschedule to make things work with multipatch. Spit out a compile-time warning since
poloidal ID B-field routine is not yet adapted to MP.
From: Christian Reisswig <reisswig@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@215 ac85fae7-cede-4708-beff-ae01c7fa1c26
Diffstat (limited to 'src')
| -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 |