[cig-commits] r13786 - mc/3D/CitcomS/trunk/lib

[tan2 at geodynamics.org]

Author: tan2
Date: 2008-12-22 10:48:50 -0800 (Mon, 22 Dec 2008)
New Revision: 13786

Modified:
   mc/3D/CitcomS/trunk/lib/Viscosity_structures.c
Log:
Fixed a bug in viscosity rheol=7. Bug reported by Laura Alisic.


Modified: mc/3D/CitcomS/trunk/lib/Viscosity_structures.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Viscosity_structures.c	2008-12-17 10:34:12 UTC (rev 13785)
+++ mc/3D/CitcomS/trunk/lib/Viscosity_structures.c	2008-12-22 18:48:50 UTC (rev 13786)
@@ -561,10 +561,10 @@
                     */
 
                     EEta[m][ (i-1)*vpts + jj ] = tempa*
-                        exp( (E->viscosity.E[l] +  E->viscosity.Z[l-1]*zzz )
+                        exp( (E->viscosity.E[l] +  E->viscosity.Z[l]*zzz )
                              / (E->viscosity.T[l] + temp)
                              - (E->viscosity.E[l] +
-                                E->viscosity.Z[l]*(one-E->sphere.ri) )
+                                E->viscosity.Z[l]*(E->sphere.ro-E->sphere.ri) )
                              / (E->viscosity.T[l] + one) );
                 }
             }