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

walter at geodynamics.org walter at geodynamics.org
Sun Feb 22 23:39:57 PST 2009


Author: walter
Date: 2009-02-22 23:39:57 -0800 (Sun, 22 Feb 2009)
New Revision: 14120

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/input/examples/rifting.xml
Log:
 r2497 at dante:  boo | 2009-02-22 23:39:14 -0800
 Update rifting.xml



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

Modified: long/3D/Gale/trunk/input/examples/rifting.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/rifting.xml	2009-02-23 07:39:52 UTC (rev 14119)
+++ long/3D/Gale/trunk/input/examples/rifting.xml	2009-02-23 07:39:57 UTC (rev 14120)
@@ -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>
@@ -327,9 +299,23 @@
       <param name="ContinuityForceVector">cont_force</param>
       <param name="killNonConvergent">false</param>
       <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">1e-2</param>
+      <param name="nonLinearTolerance">1.0e-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="buoyancyForceTerm">
       <param name="Type">BuoyancyForceTerm</param>
       <param name="ForceVector">mom_force</param>
@@ -354,13 +340,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="boxShape">
       <param name="Type">Box</param>
@@ -405,7 +384,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>
 
@@ -421,7 +401,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>
@@ -479,7 +459,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>	
 
@@ -622,17 +603,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>
@@ -643,6 +613,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>
@@ -679,7 +653,8 @@
           <struct>
             <param name="name">vx</param>
             <param name="type">double</param>
-            <param name="value">-4e-10</param>
+            <param name="value">-4e16</param>
+<!--             <param name="value">-4e-10</param> -->
           </struct>
           <struct>
             <param name="name">vy</param>
@@ -696,7 +671,8 @@
             <param name="name">vx</param>
             <param name="type">double</param>
             <!-- 1 cm/yr -->
-            <param name="value">4e-10</param>
+            <param name="value">4e16</param>
+<!--             <param name="value">4e-10</param> -->
           </struct>
           <struct>
             <param name="name">vy</param>
@@ -820,7 +796,8 @@
 <!--   <param name="checkpointEvery">1</param> -->
   <param name="dtFactor">1.0</param>
 
-  <param name="ExtensionFactor">1.4e-13</param>
+  <param name="ExtensionFactor">1.4e13</param>
+<!--   <param name="ExtensionFactor">1.4e-13</param> -->
   <param name="ExtensionCentreX">20000</param>
 
   <param name="linearTolerance">1e-5</param>
@@ -831,8 +808,8 @@
   <param name="StepFunctionUpperOffset">-100000000</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