diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-11-09 01:59:11 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-11-09 01:59:11 +0000 |
| commit | 08d4b0f807cb5d7c0a9c6577cbc612b7262ed508 (patch) | |
| tree | 9d5a5cb2b52b5c53e71c516b0844e134d89a05dc | |
| parent | a0d787cd49f3782f29b231ae26f51ed1d1c0642a (diff) | |
GRHydro_InitData: Force mag field to go to zero at the puncture.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@180 ac85fae7-cede-4708-beff-ae01c7fa1c26
| -rw-r--r-- | src/GRHydro_BondiM_new.F90 | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/src/GRHydro_BondiM_new.F90 b/src/GRHydro_BondiM_new.F90 index 325e158..9d4b661 100644 --- a/src/GRHydro_BondiM_new.F90 +++ b/src/GRHydro_BondiM_new.F90 @@ -238,6 +238,7 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS) rho(i,j,k) = rhocheck+drhodr*riso*(riso-M)/M else rho(i,j,k) = (rhocheck-drhodr*M/4.0)*(1.-cos(2.*M_PI*riso/M))/2.0 + bondi_bmag = bondi_bmag*8.0d0*riso**3 endif utmp = Kval * rho(i,j,k)**( gam ) / (gam - 1.) |