[cig-commits] r15971 - mc/3D/CitcomCU/trunk/src

[tan2 at geodynamics.org]

Author: tan2
Date: 2009-11-16 10:06:16 -0800 (Mon, 16 Nov 2009)
New Revision: 15971

Modified:
   mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
Log:
Added rheol=11, eta = eta0 * exp(E/(T+T0) - E/(0.5+T0))


Modified: mc/3D/CitcomCU/trunk/src/Viscosity_structures.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2009-11-15 17:46:27 UTC (rev 15970)
+++ mc/3D/CitcomCU/trunk/src/Viscosity_structures.c	2009-11-16 18:06:16 UTC (rev 15971)
@@ -499,6 +499,29 @@
 		  }
 	      }
 	    break;
+	  case 11:
+            /* 
+               eta = eta0 * exp(E/(T+T0) - E/(0.5+T0))
+            */
+            for(i = 1; i <= nel; i++)
+              {
+                l = E->mat[i] - 1 ;
+                tempa = E->viscosity.N0[l];
+
+                for(kk = 1; kk <= ends; kk++)
+                  TT[kk] = E->T[E->ien[i].node[kk]];
+
+                for(jj = 1; jj <= vpts; jj++)
+                  {
+                    temp = 1.0e-32;
+                    for(kk = 1; kk <= ends; kk++)
+                      {
+                        temp += max(zero, TT[kk]) * E->N.vpt[GNVINDEX(kk, jj)];;
+                      }
+                    EEta[(i - 1) * vpts + jj] = tempa * exp((E->viscosity.E[l]) / (temp + E->viscosity.T[l]) - (E->viscosity.E[l]) / (0.5 + E->viscosity.T[l]));
+                  }
+              }
+            break;
 	  default:
 	    myerror(E,"RHEOL option undefined");
 	    break;