From a25f8ba2499d399d44916f12607167e597076298 Mon Sep 17 00:00:00 2001 From: rhaas Date: Mon, 17 Sep 2012 16:27:26 +0000 Subject: GRHydro_InitData: relax on the tiny value hard coded. From: Bruno Coutinho Mundim git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@165 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_Bondi_new.F90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/GRHydro_Bondi_new.F90 b/src/GRHydro_Bondi_new.F90 index 63243eb..a2eabcd 100755 --- a/src/GRHydro_Bondi_new.F90 +++ b/src/GRHydro_Bondi_new.F90 @@ -52,7 +52,7 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS) CCTK_INT :: i, j, k, nx, ny, nz, imin, jb,N_points CCTK_REAL :: ONEmTINY, tiny - PARAMETER (N_points=2000,ONEmTINY=0.999999d0,tiny=1.0d-14) + PARAMETER (N_points=2000,ONEmTINY=0.999999d0,tiny=1.0d-12) CCTK_REAL :: M, Msq, Mdot, rs, gam, rmin_bondi, rmax_bondi, cs_sq,cs,vs_sq,vs,rhos,gtemp,hs, Kval, Qdot CCTK_REAL :: logrmin,dlogr,rhotmp,utmp,vtmp,rspher CCTK_REAL :: r_bondi(N_points), logr_bondi(N_points), rho_bondi(N_points), u_bondi(N_points), v_bondi(N_points) @@ -224,8 +224,8 @@ subroutine GRHydro_Bondi_Iso(CCTK_ARGUMENTS) ! if(riso.gt.1.014d0.and.riso.lt.1.015) then if(abs(z(i,j,k)).lt.tiny.and. & - x(i,j,k).gt.6.7d0.and.x(i,j,k).lt.6.9.and. & - y(i,j,k).gt.1.34d0.and.y(i,j,k).lt.1.38 ) then + x(i,j,k).gt.6.7d0.and.x(i,j,k).lt.6.9d0.and. & + y(i,j,k).gt.1.34d0.and.y(i,j,k).lt.1.38d0 ) then write(debug_message,'(a,15f22.14)') 'Point to check:', & x(i,j,k),y(i,j,k),z(i,j,k),riso,gxx(i,j,k),dens(i,j,k),& tau(i,j,k),sx(i,j,k),sy(i,j,k),sz(i,j,k),rho(i,j,k),eps(i,j,k),& -- cgit v1.2.3