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@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.)