[cig-commits] r4179 - in long/3D/Gale/trunk/src/Underworld: . Rheology/src

[walter at geodynamics.org]

Author: walter
Date: 2006-08-01 01:55:59 -0700 (Tue, 01 Aug 2006)
New Revision: 4179

Modified:
   long/3D/Gale/trunk/src/Underworld/
   long/3D/Gale/trunk/src/Underworld/Rheology/src/Arrhenius.c
Log:
 r428 at earth:  boo | 2006-08-01 01:53:43 -0700
  r408 at earth (orig r279):  PatrickSunter | 2006-07-25 19:49:52 -0700
  Set a very small tolerance on the activationVolume in the
  Arrhenius rheology, since its a double.
  
 



Property changes on: long/3D/Gale/trunk/src/Underworld
___________________________________________________________________
Name: svk:merge
   - 9570c393-cf10-0410-b476-9a651db1e55a:/cig:427
c24a034b-ab11-0410-afe6-cfe714e2959e:/trunk:278
   + 9570c393-cf10-0410-b476-9a651db1e55a:/cig:428
c24a034b-ab11-0410-afe6-cfe714e2959e:/trunk:279

Modified: long/3D/Gale/trunk/src/Underworld/Rheology/src/Arrhenius.c
===================================================================
--- long/3D/Gale/trunk/src/Underworld/Rheology/src/Arrhenius.c	2006-08-01 08:55:56 UTC (rev 4178)
+++ long/3D/Gale/trunk/src/Underworld/Rheology/src/Arrhenius.c	2006-08-01 08:55:59 UTC (rev 4179)
@@ -174,7 +174,7 @@
 	FeVariable_InterpolateWithinElement( temperatureField, lElement_I, xi, &temperature );
 	
 	/* If activationVolume is 0 there is no need to calculate the depth of the particle see viscosity line below. */
-	if( activationVolume != 0.0 ) {
+	if( activationVolume > (0.0 + 1e-12 )  ) {
 		/* Calculate Depth */
 		height = self->blockGeometry->max[ J_AXIS ];