[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 380) - Commits

4 Jul 2012

User: rhaas
Date: 2012/07/05 12:36 AM
Modified:
 /trunk/src/
  GRHydro_CalcUpdate.F90
Log:
 * change comparison of double with integer to double with double
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_CalcUpdate.F90
Delta lines: +2 -2
===================================================================

--- trunk/src/GRHydro_CalcUpdate.F90	2012-07-05 05:36:38 UTC (rev 379)
+++ trunk/src/GRHydro_CalcUpdate.F90	2012-07-05 05:36:40 UTC (rev 380)
@@ -358,8 +358,8 @@
! Eliminate non-radial fluid velocities to obtain pseudo 1D scheme
                rnorm = (x(i,j,k)**2 + y(i,j,k)**2 + z(i,j,k)**2)
-               if (rnorm.eq.0) then
-                  rnormI = 0
+               if (rnorm.lt.1.0d-10) then
+                  rnormI = 0.0d0
                else
                   rnormI = 1.0d0/rnorm
                endif

    