[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 180)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2012/11/08 07:59 PM

Modified:
 /trunk/src/
  GRHydro_BondiM_new.F90

Log:
 GRHydro_InitData: Force mag field to go to zero at the puncture.
 
 From: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_BondiM_new.F90
Delta lines: +1 -0
===================================================================
--- trunk/src/GRHydro_BondiM_new.F90	2012-11-09 01:59:09 UTC (rev 179)
+++ trunk/src/GRHydro_BondiM_new.F90	2012-11-09 01:59:11 UTC (rev 180)
@@ -238,6 +238,7 @@
                  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.)