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

walter at geodynamics.org walter at geodynamics.org
Wed Mar 25 06:39:32 PDT 2009


Author: walter
Date: 2009-03-25 06:39:31 -0700 (Wed, 25 Mar 2009)
New Revision: 14447

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/examples/tibet3D.xml
Log:
 r2588 at dante:  boo | 2009-03-25 06:39:12 -0700
 Make tibet3D work



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

Modified: long/3D/Gale/trunk/input/examples/tibet3D.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet3D.xml	2009-03-25 13:39:25 UTC (rev 14446)
+++ long/3D/Gale/trunk/input/examples/tibet3D.xml	2009-03-25 13:39:31 UTC (rev 14447)
@@ -4,15 +4,6 @@
     <struct name="conditionFunctions">
       <param name="Type">StgFEM_StandardConditionFunctions</param>
     </struct>
-    <struct name="mesh-constant">
-      <param name="Type">FeMesh</param>
-      <param name="elementType">constant</param>
-    </struct>
-    <struct name="constantMesh-generator">
-      <param name="Type">C0Generator</param>
-      <param name="mesh">mesh-constant</param>
-      <param name="elementMesh">mesh-linear</param>
-    </struct>
     <struct name="mesh-linear">
       <param name="Type">FeMesh</param>
       <param name="elementType">linear</param>
@@ -108,42 +99,23 @@
     </struct>
     <struct name="pressure">
       <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</param>
+      <param name="mesh">mesh-linear</param>
       <param name="Rank">Scalar</param>
       <param name="DataType">Double</param>
     </struct>
     <struct name="pressureDofLayout">
       <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-constant</param>
+      <param name="mesh">mesh-linear</param>
       <list name="BaseVariables">
         <param>pressure</param>
       </list>
     </struct>
     <struct name="PressureField">
       <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
+      <param name="FEMesh">mesh-linear</param>
       <param name="DofLayout">pressureDofLayout</param>
       <param name="LinkedDofInfo">pressureLinkedDofs</param>
     </struct>
-    <struct name="nodalPressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Scalar</param>
-      <param name="DataType">Double</param>
-    </struct>
-    <struct name="nodalPressureDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <list name="BaseVariables">
-        <param>nodalPressure</param>
-      </list>
-    </struct>
-    <struct name="NodalPressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">nodalPressureDofLayout</param>
-      <param name="LinkedDofInfo">nodalPressureLinkedDofs</param>
-    </struct>
     <struct name="StressField">
       <param name="Type">StressField</param>
       <param name="StrainRateField">StrainRateField</param>
@@ -330,12 +302,51 @@
       <param name="nonLinearTolerance">1e-2</param>
       <param name="makeConvergenceFile">false</param>
     </struct>
+
+    <struct name="c_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">PressureField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">True</param>
+    </struct>
+    <struct name="mixedStabiliser">
+      <param name="Type">MixedStabiliserTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="picSwarm">picIntegrationPoints</param>
+      <param name="storeVisc">storeViscosity</param>
+      <param name="StiffnessMatrix">c_matrix</param>
+    </struct>
+    <struct name="hydrostaticTerm">
+      <param name="Type">HydrostaticTerm</param>
+      <param name="upperDensity">2800</param>
+      <param name="upperAlpha">3.0e-5</param>
+      <param name="lowerDensity">3300</param>
+      <param name="lowerAlpha">3.0e-5</param>
+      <param name="height">maxY</param>
+      <param name="materialBoundary">68000</param>
+      <param name="T_0">TemperatureProfileTop</param>
+      <param name="linearCoefficient">TemperatureProfileLinearCoefficient</param>
+      <param name="exponentialCoefficient1">TemperatureProfileExponentialCoefficient1</param>
+      <param name="exponentialCoefficient2">TemperatureProfileExponentialCoefficient2</param>
+      <param name="gravity">gravity</param>
+    </struct>
+
+    <struct name="stressBC">
+      <param name="Type">StressBC</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="wall">top</param>
+      <param name="y_type">HydrostaticTerm</param>
+      <param name="y_value">hydrostaticTerm</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>
+      <param name="HydrostaticTerm">hydrostaticTerm</param>
     </struct>
     <struct name="background">
       <param name="Type">Everywhere</param>
@@ -354,13 +365,6 @@
         <param>true</param>
       </list>
     </struct>
-    <struct name="pressureRemesher">
-      <param name="Type">CellRemesher</param>
-      <param name="mesh">mesh-constant</param>
-      <param name="meshType">regular</param>
-      <param name="dim">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
 
     <struct name="bottomBoundaryShape">
       <param name="Type">Box</param>
@@ -418,7 +422,8 @@
       <param name="Type">FrankKamenetskii</param>
       <param name="TemperatureField">TemperatureField</param>
       <!-- at T=273, viscosity=1e25, at T=1273, viscosity=1e20 -->
-      <param name="eta0">2.31739465e26</param>
+      <param name="eta0">2.31739465</param>
+<!--       <param name="eta0">2.31739465e26</param> -->
       <param name="theta">0.011512925465</param>
     </struct>	
 
@@ -434,7 +439,7 @@
     </struct>
     <struct name="yielding">
       <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">NodalPressureField</param>
+      <param name="PressureField">PressureField</param>
       <param name="VelocityGradientsField">VelocityGradientsField</param>
       <param name="MaterialPointsSwarm">materialSwarm</param>
       <param name="Context">context</param>
@@ -477,7 +482,8 @@
     <struct name="mantleViscosity">
       <param name="Type">FrankKamenetskii</param>
       <param name="TemperatureField">TemperatureField</param>
-      <param name="eta0">2.31739465e26</param>
+      <param name="eta0">2.31739465</param>
+<!--       <param name="eta0">2.31739465e26</param> -->
       <param name="theta">0.011512925465</param>
     </struct>	
 
@@ -582,7 +588,6 @@
   </struct>
   <list name="plugins">
     <param>Underworld_EulerDeform</param>
-<!--     <param>Underworld_DumpSwarm</param> -->
     <param>Underworld_VTKOutput</param>
   </list>
   <param name="maxTimeSteps">10</param>
@@ -612,17 +617,6 @@
   <struct name="EulerDeform">
     <list name="systems">
       <struct>
-        <param name="mesh">mesh-constant</param>
-        <param name="remesher">pressureRemesher</param>
-        <param name="velocityField">VelocityField</param>
-        <list name="fields">
-          <struct>
-            <param name="field">PressureField</param>
-            <param name="variable">pressure</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
         <param name="mesh">mesh-linear</param>
         <param name="remesher">velocityRemesher</param>
         <param name="velocityField">VelocityField</param>
@@ -633,6 +627,10 @@
             <param name="field">VelocityField</param>
             <param name="variable">velocity</param>
           </struct>
+          <struct>
+            <param name="field">PressureField</param>
+            <param name="variable">pressure</param>
+          </struct>
         </list>
       </struct>
     </list>
@@ -659,7 +657,8 @@
             <param name="name">vz</param>
             <param name="type">double</param>
             <!-- 1 cm/yr -->
-            <param name="value">-3.17097919838e-10</param>
+            <param name="value">-3.17097919838e16</param>
+<!--             <param name="value">-3.17097919838e-10</param> -->
           </struct>
         </list>
       </struct>
@@ -807,8 +806,8 @@
 <!--   <param name="checkpointEvery">1</param> -->
   <param name="dtFactor">1.0</param>
 
-<!--   <param name="journal.info">True</param> -->
-<!--   <param name="journal.debug">True</param> -->
-<!--   <param name="journal-level.info">2</param> -->
-<!--   <param name="journal-level.debug">2</param> -->
+  <param name="journal.info">True</param>
+  <param name="journal.debug">True</param>
+  <param name="journal-level.info">2</param>
+  <param name="journal-level.debug">2</param>
 </StGermainData>



More information about the CIG-COMMITS mailing list