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

[rhaas at tapir.caltech.edu]

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