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

walter at geodynamics.org walter at geodynamics.org
Fri Nov 6 13:20:56 PST 2009


Author: walter
Date: 2009-11-06 13:20:56 -0800 (Fri, 06 Nov 2009)
New Revision: 15938

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Underworld/Rheology/src/NonNewtonian.c
Log:
 r2691 at dante:  boo | 2009-11-06 13:20:03 -0800
 Fix a factor of two error in NonNewtonian



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2689
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2691

Modified: long/3D/Gale/trunk/src/Underworld/Rheology/src/NonNewtonian.c
===================================================================
--- long/3D/Gale/trunk/src/Underworld/Rheology/src/NonNewtonian.c	2009-11-04 22:26:19 UTC (rev 15937)
+++ long/3D/Gale/trunk/src/Underworld/Rheology/src/NonNewtonian.c	2009-11-06 21:20:56 UTC (rev 15938)
@@ -205,7 +205,7 @@
         Journal_Firewall(n!=0 && T!=0 && A!=0,errorStream,
                          "Error in NonNewtonian: T, n, and A must all be non-zero:\n\tT=%g\n\tn=%g\n\tA=%g\n",T,n,A);
 
-        viscosity=exp(T_0/(2*n*T))/(2*A);
+        viscosity=exp(T_0/(n*T))/(2*A);
 
 	if ( fabs( n - 1.0 ) > 1.0e-10 )
           {



More information about the CIG-COMMITS mailing list