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

[walter at geodynamics.org]

Author: walter
Date: 2007-09-24 12:55:35 -0700 (Mon, 24 Sep 2007)
New Revision: 8014

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/examples/tibet.xml
Log:
 r1934 at earth:  boo | 2007-09-24 12:42:49 -0700
 Make tibet have the right size and normal stress bc



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

Modified: long/3D/Gale/trunk/input/examples/tibet.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet.xml	2007-09-24 19:55:31 UTC (rev 8013)
+++ long/3D/Gale/trunk/input/examples/tibet.xml	2007-09-24 19:55:35 UTC (rev 8014)
@@ -299,7 +299,7 @@
       <param name="Type">BuoyancyForceTerm</param>
       <param name="ForceVector">mom_force</param>
       <param name="Swarm">picIntegrationPoints</param>
-      <param name="gravity">1.0</param>
+      <param name="gravity">9.8</param>
     </struct>
     <struct name="background">
       <param name="Type">Everywhere</param>
@@ -331,11 +331,32 @@
       <param name="startX">minX</param>
       <param name="endX">maxX</param>
       <param name="startY">minY</param>
-      <param name="endY">0.05</param>
+      <param name="endY">10000</param>
       <param name="startZ">minZ</param>
       <param name="endZ">maxZ</param>
     </struct>
 
+    <struct name="stressBC">
+      <param name="Type">StressBC</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="wall">bottom</param>
+      <param name="y_type">double</param>
+      <!-- This number is just adding up rho*g*h at the edge -->
+      <param name="y_value">3077200000</param>
+    </struct>
+    <struct name="surfaceAdaptor">
+      <param name="Type">SurfaceAdaptor</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="sourceGenerator">linearMesh-generator</param>
+      <param name="surfaceType">plateau</param>
+      <param name="x1">300000</param>
+      <param name="x2">350000</param>
+      <param name="x3">650000</param>
+      <param name="x4">700000</param>
+      <param name="height">5000</param>
+    </struct>
+
     <struct name="boxShape">
       <param name="Type">Box</param>
       <param name="startX">minX</param>
@@ -357,36 +378,53 @@
       <param name="Type">StoreStress</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
     </struct>
-    <struct name="viscous">
+
+    <!-- Crust -->
+    <struct name="crustShape">
+      <param name="Type">PolygonShape</param>
+      <list name="verticies">
+        <asciidata>
+          <columnDefinition name = "x" type="double"/>
+          <columnDefinition name = "y" type="double"/>
+          0         68000
+          300000    68000
+          350000    53000
+          650000    53000
+          700000    68000
+          1000000   68000
+          1000000   100000
+          0         100000
+        </asciidata>
+      </list>
+    </struct>
+    <struct name="crustViscous">
       <param name="Type">RheologyMaterial</param>
-      <param name="Shape">boxShape</param>
-      <param name="density">1.0</param>
+      <param name="Shape">crustShape</param>
+      <param name="density">2800</param>
       <list name="Rheology">
         <param>backgroundViscosity</param>
         <param>storeViscosity</param>
         <param>storeStress</param>
       </list>
     </struct>
-
-    <struct name="stressBC">
-      <param name="Type">StressBC</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="wall">bottom</param>
-      <param name="y_type">double</param>
-      <param name="y_value">0.35</param>
+    
+    <!-- Mantle -->
+    <struct name="mantleShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>boxShape</param>
+        <param>!crustShape</param>
+      </list>
     </struct>
-
-    <struct name="surfaceAdaptor">
-      <param name="Type">SurfaceAdaptor</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="sourceGenerator">linearMesh-generator</param>
-      <param name="surfaceType">plateau</param>
-      <param name="x1">300</param>
-      <param name="x2">350</param>
-      <param name="x3">650</param>
-      <param name="x4">700</param>
-      <param name="height">10</param>
+    <struct name="mantleViscous">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">mantleShape</param>
+      <param name="density">3300</param>
+      <list name="Rheology">
+        <param>backgroundViscosity</param>
+        <param>storeViscosity</param>
+        <param>storeStress</param>
+      </list>
     </struct>
 
   </struct>
@@ -404,9 +442,9 @@
   <param name="minX">0.0f</param>
   <param name="minY">0.0f</param>
   <param name="minZ">0.0f</param>
-  <param name="maxX">1000.0f</param>
-  <param name="maxY">100f</param>
-  <param name="maxZ">1000f</param>
+  <param name="maxX">1000000.0f</param>
+  <param name="maxY">100000f</param>
+  <param name="maxZ">1000000f</param>
   <param name="elementResI">32</param>
   <param name="elementResJ">16</param>
   <param name="elementResK">4</param>