From 70eafbdf9a9b0c54dd00d543d3cb0142fa7e9c5c Mon Sep 17 00:00:00 2001 From: rhaas Date: Fri, 9 Nov 2012 01:59:16 +0000 Subject: GRHydro_InitData: yet another initialization bug fixed in bondiM. From: Bruno Coutinho Mundim git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@182 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_BondiM_new.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/src/GRHydro_BondiM_new.F90 b/src/GRHydro_BondiM_new.F90 index 1c61374..346bc19 100644 --- a/src/GRHydro_BondiM_new.F90 +++ b/src/GRHydro_BondiM_new.F90 @@ -201,7 +201,6 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS) - bondi_bsmooth = 1.0d0 do i=1,nx do j=1,ny @@ -216,6 +215,8 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS) yhat = yp/riso zhat = zp/riso roverm = riso/M + + bondi_bsmooth = 1.0d0 if(roverm > ONEmTINY) then rsch = 0.25 * ( 2.*riso + M)**2 / riso -- cgit v1.2.3