[cig-commits] r14542 - long/3D/Gale/trunk

[walter at geodynamics.org]

Author: walter
Date: 2009-03-30 02:44:16 -0700 (Mon, 30 Mar 2009)
New Revision: 14542

Added:
   long/3D/Gale/trunk/UPGRADE
Modified:
   long/3D/Gale/trunk/
Log:
 r2636 at dante:  boo | 2009-03-30 02:42:51 -0700
 Add UPGRADE



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

Added: long/3D/Gale/trunk/UPGRADE
===================================================================
--- long/3D/Gale/trunk/UPGRADE	                        (rev 0)
+++ long/3D/Gale/trunk/UPGRADE	2009-03-30 09:44:16 UTC (rev 14542)
@@ -0,0 +1,105 @@
+To upgrade your files from Gale 1.3 or 1.3.1 to 1.4, you need to take
+out the mesh-contant structs at the beginning.
+
+    <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>
+
+Next, replace the references "mesh-constant" in "pressure",
+"pressureDofLayout", and "PressureField" with "mesh-linear".
+
+Remove the nodal pressures
+
+    <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>
+
+Replace the reference to NodalPressureField in any of your MohrCoulomb
+rheologies to PressureField.  At the same time, you need to rename all
+of the MohrCoulomb structs to DruckerPrager.
+
+Remove the "pressureRemesher" 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>
+
+Remove any references to UnderworldDumpSwarm.
+
+In EulerDeform, remove the constant mesh remesher
+
+      <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>
+
+and add PressureField to the linear mesh remesher.  So it would look like
+
+      <struct>
+        <param name="mesh">mesh-linear</param>
+        <param name="remesher">velocityRemesher</param>
+        <param name="velocityField">VelocityField</param>
+        <param name="wrapTop">True</param>
+        <list name="fields">
+          <struct>
+            <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>
+
+Add the stabilization term right after "stokesEqn"
+
+    <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>
+
+You may also want to add a HydrostaticTerm.  See the manual for details.