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

[tan2 at geodynamics.org]

Author: tan2
Date: 2008-11-07 15:32:00 -0800 (Fri, 07 Nov 2008)
New Revision: 13272

Modified:
   mc/3D/CitcomS/trunk/lib/Viscosity_structures.c
Log:
Reserved 'rheol=100' for user-defined viscosity law.


Modified: mc/3D/CitcomS/trunk/lib/Viscosity_structures.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Viscosity_structures.c	2008-11-07 23:31:29 UTC (rev 13271)
+++ mc/3D/CitcomS/trunk/lib/Viscosity_structures.c	2008-11-07 23:32:00 UTC (rev 13272)
@@ -382,7 +382,7 @@
         break;
 
     case 4:
-
+        /* eta = N_0 exp( (E + (1-z)Z_0) / (T+T_0) ) */
         for(m=1;m<=E->sphere.caps_per_proc;m++)
             for(i=1;i<=nel;i++)   {
                 l = E->mat[m][i] - 1;
@@ -618,8 +618,22 @@
             }
         break;
 
+    case 100:
+        /* user-defined viscosity law goes here */
+        fprintf(stderr, "Need user definition for viscosity law: 'rheol=%d'\n",
+                E->viscosity.RHEOL);
+        parallel_process_termination();
+        break;
+
+    default:
+        /* unknown option */
+        fprintf(stderr, "Invalid value of 'rheol=%d'\n", E->viscosity.RHEOL);
+
+        parallel_process_termination();
+        break;
     }
 
+
     return;
 }