diff options
-rw-r--r-- | src/GRHydro_BondiM_new.F90 | 3 |
1 files changed, 2 insertions, 1 deletions
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 |