[cig-commits] r12478 - in long/3D/Gale/trunk: . input/examples

[walter at geodynamics.org]

Author: walter
Date: 2008-07-25 07:15:02 -0700 (Fri, 25 Jul 2008)
New Revision: 12478

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/examples/rifting.xml
Log:
 r2284 at earth:  boo | 2008-07-25 07:12:58 -0700
 Make rifting use 2 layers and turn on gravity



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

Modified: long/3D/Gale/trunk/input/examples/rifting.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/rifting.xml	2008-07-25 07:51:33 UTC (rev 12477)
+++ long/3D/Gale/trunk/input/examples/rifting.xml	2008-07-25 14:15:02 UTC (rev 12478)
@@ -330,13 +330,13 @@
       <param name="nonLinearTolerance">1e-2</param>
       <param name="makeConvergenceFile">false</param>
     </struct>
-<!--     <struct name="buoyancyForceTerm"> -->
-<!--       <param name="Type">BuoyancyForceTerm</param> -->
-<!--       <param name="ForceVector">mom_force</param> -->
-<!--       <param name="TemperatureField">TemperatureField</param> -->
-<!--       <param name="Swarm">picIntegrationPoints</param> -->
-<!--       <param name="gravity">9.81</param> -->
-<!--     </struct> -->
+    <struct name="buoyancyForceTerm">
+      <param name="Type">BuoyancyForceTerm</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="TemperatureField">TemperatureField</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="gravity">9.81</param>
+    </struct>
     <struct name="background">
       <param name="Type">Everywhere</param>
     </struct>
@@ -382,6 +382,16 @@
       <param name="endZ">maxZ</param>
     </struct>
 
+    <struct name="crustShape">
+      <param name="Type">Box</param>
+      <param name="startX">minX</param>
+      <param name="endX">maxX</param>
+      <param name="startY">15000</param>
+      <param name="endY">maxY</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+
     <struct name="storeViscosity">
       <param name="Type">StoreVisc</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
@@ -422,9 +432,9 @@
       <param name="frictionCoefficientAfterSoftening">0.57735026919</param>
     </struct>
 
-    <struct name="crustViscous">
+    <struct name="crust">
       <param name="Type">RheologyMaterial</param>
-      <param name="Shape">boxShape</param>
+      <param name="Shape">crustShape</param>
       <param name="density">2800</param>
       <param name="alpha">3.0e-5</param>
       <param name="diffusivity">1.0e-6</param>
@@ -441,6 +451,7 @@
         <param>storeStress</param>
       </list>
     </struct>
+
     
     <struct name="divergenceForce">
       <param name="Type">DivergenceForce</param>
@@ -452,8 +463,43 @@
     </struct>
     
 
+    
+    <!-- Mantle -->
+    <struct name="mantleShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>boxShape</param>
+        <param>!crustShape</param>
+      </list>
+    </struct>
 
+    <struct name="mantleViscosity">
+      <param name="Type">FrankKamenetskii</param>
+      <param name="TemperatureField">TemperatureField</param>
+      <param name="eta0">2.31739465e26</param>
+      <param name="theta">0.011512925465</param>
+    </struct>	
 
+    <struct name="mantleViscous">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">mantleShape</param>
+      <param name="density">3300</param>
+      <param name="alpha">3.0e-5</param>
+      <param name="diffusivity">1.0e-6</param>
+      <list name="heatingElements">
+        <struct>
+          <param name="Q">8.48484848485e-13</param>
+          <param name="lambda">0.0</param>
+        </struct>
+      </list>
+      <list name="Rheology">
+        <param>mantleViscosity</param>
+        <param>storeViscosity</param>
+        <param>storeStress</param>
+      </list>
+    </struct>
+
+
     <!-- Temperature components -->
     <struct name="temperature">
       <param name="Type">MeshVariable</param>