From 08d4b0f807cb5d7c0a9c6577cbc612b7262ed508 Mon Sep 17 00:00:00 2001 From: rhaas Date: Fri, 9 Nov 2012 01:59:11 +0000 Subject: GRHydro_InitData: Force mag field to go to zero at the puncture. From: Bruno Coutinho Mundim git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@180 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_BondiM_new.F90 | 1 + 1 file changed, 1 insertion(+) 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.) -- cgit v1.2.3