diff options
author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-11-09 01:59:16 +0000 |
---|---|---|
committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2012-11-09 01:59:16 +0000 |
commit | 70eafbdf9a9b0c54dd00d543d3cb0142fa7e9c5c (patch) | |
tree | 3a8e8576557868c025f069f1bb407854459f38a8 | |
parent | 34859874794b371642fa6489b24e6ffc808d1a8f (diff) |
GRHydro_InitData: yet another initialization bug fixed in bondiM.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@182 ac85fae7-cede-4708-beff-ae01c7fa1c26
-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 |