User: rhaas Date: 2012/11/08 07:59 PM
Modified: /trunk/src/ GRHydro_BondiM_new.F90
Log: GRHydro_InitData: Fix B-field smoothing in BondiM routine.
From: Bruno Coutinho Mundim bcmsma@astro.rit.edu
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_BondiM_new.F90 Delta lines: +9 -10 =================================================================== --- trunk/src/GRHydro_BondiM_new.F90 2012-11-09 01:59:11 UTC (rev 180) +++ trunk/src/GRHydro_BondiM_new.F90 2012-11-09 01:59:14 UTC (rev 181) @@ -65,6 +65,7 @@ CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck CCTK_REAL :: uiso, uisocheck, vcheck, ucheck2, vcheck2, xhat,yhat, zhat, xp, yp, zp CCTK_REAL :: f,df,ddf,a,b,c,rsm,roverm,dudr,uisocheck2,auiso,buiso + CCTK_REAL :: bondi_bsmooth
character(400) :: debug_message
@@ -199,11 +200,9 @@ nz = cctk_lsh(3)
- write(debug_message,'(a,3f22.14)') 'Lower bound coordinates:',x(1,1,1),y(1,1,1),z(1,1,1) - call CCTK_INFO(debug_message) - write(debug_message,'(a,3f22.14)') 'Upper bound coordinates:',x(nx,ny,nz),y(nx,ny,nz),z(nx,ny,nz) - call CCTK_INFO(debug_message) - + + bondi_bsmooth = 1.0d0 + do i=1,nx do j=1,ny do k=1,nz @@ -238,7 +237,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 + bondi_bsmooth = 8.0d0*riso**3 endif utmp = Kval * rho(i,j,k)**( gam ) / (gam - 1.)
@@ -259,10 +258,10 @@
det=SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))
- Bvecx(i,j,k) = bondi_bmag*M**2*xhat/sqrt(det)/riso**2 - Bvecy(i,j,k) = bondi_bmag*M**2*yhat/sqrt(det)/riso**2 - Bvecz(i,j,k) = bondi_bmag*M**2*zhat/sqrt(det)/riso**2 - + Bvecx(i,j,k) = bondi_bsmooth*bondi_bmag*M**2*xhat/sqrt(det)/riso**2 + Bvecy(i,j,k) = bondi_bsmooth*bondi_bmag*M**2*yhat/sqrt(det)/riso**2 + Bvecz(i,j,k) = bondi_bsmooth*bondi_bmag*M**2*zhat/sqrt(det)/riso**2 + call Prim2ConGenM(GRHydro_eos_handle,gxx(i,j,k),gxy(i,j,k), & gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k), & det, dens(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k), &