[cig-commits] commit: Change NonNewtonian to use the convential scaling for A

[Mercurial]

changeset:   857:2612237febd1
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Oct 20 15:25:41 2010 -0700
files:       Rheology/src/NonNewtonian.c
description:
Change NonNewtonian to use the convential scaling for A


diff -r 20b6b32a9eef -r 2612237febd1 Rheology/src/NonNewtonian.c
--- a/Rheology/src/NonNewtonian.c	Wed Oct 20 10:43:46 2010 -0700
+++ b/Rheology/src/NonNewtonian.c	Wed Oct 20 15:25:41 2010 -0700
@@ -215,7 +215,7 @@ void _NonNewtonian_ModifyConstitutiveMat
         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/(n*T))/(2*A);
+        viscosity=exp(T_0/(n*T))*pow(A,-1/n)/2;
 
 	if ( fabs( n - 1.0 ) > 1.0e-10 )
           {