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

walter at geodynamics.org walter at geodynamics.org
Wed Jan 17 16:36:04 PST 2007


Author: walter
Date: 2007-01-17 16:36:03 -0800 (Wed, 17 Jan 2007)
New Revision: 5818

Added:
   long/3D/Gale/trunk/input/
   long/3D/Gale/trunk/input/extension.xml
   long/3D/Gale/trunk/input/fountain.xml
   long/3D/Gale/trunk/input/shortening.xml
   long/3D/Gale/trunk/input/sinusoid.xml
   long/3D/Gale/trunk/input/subduction.xml
Removed:
   long/3D/Gale/trunk/extension.xml
   long/3D/Gale/trunk/fountain.xml
   long/3D/Gale/trunk/shortening.xml
   long/3D/Gale/trunk/sinusoid.xml
   long/3D/Gale/trunk/subduction.xml
Modified:
   long/3D/Gale/trunk/
Log:
 r1437 at earth:  boo | 2007-01-17 16:34:08 -0800
 Move input files to input/



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

Deleted: long/3D/Gale/trunk/extension.xml
===================================================================
--- long/3D/Gale/trunk/extension.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/extension.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -1,535 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-  <struct name="components">
-    <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>
-    </struct>
-    <struct name="linearMesh-generator">
-      <param name="Type">CartesianGenerator</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="dims">dim</param>
-      <param name="shadowDepth">shadowDepth</param>
-      <list name="size">
-        <param>elementResI</param>
-        <param>elementResJ</param>
-        <param>elementResK</param>
-      </list>
-      <list name="minCoord">
-        <param>minX</param>
-        <param>minY</param>
-        <param>minZ</param>
-      </list>
-      <list name="maxCoord">
-        <param>maxX</param>
-        <param>maxY</param>
-        <param>maxZ</param>
-      </list>
-    </struct>
-    <struct name="velocity">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Vector</param>
-      <param name="DataType">Double</param>
-      <param name="VectorComponentCount">dim</param>
-      <list name="names">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="velocityBCs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityICs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="BaseVariableCount">dim</param>
-      <list name="BaseVariables">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="VelocityField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">velocityDofLayout</param>
-      <param name="BC">velocityBCs</param>
-      <param name="IC">velocityICs</param>
-      <param name="LinkedDofInfo">velocityLinkedDofs</param>
-    </struct>
-    <struct name="VelocityMagnitudeField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Magnitude</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Gradient</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorInvariant</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorSymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="VorticityField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorAntisymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">SymmetricTensor_Invariant</param>
-      <param name="FeVariable">StrainRateField</param>
-    </struct>
-    <struct name="pressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</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>
-      <list name="BaseVariables">
-        <param>pressure</param>
-      </list>
-    </struct>
-    <struct name="PressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
-      <param name="DofLayout">pressureDofLayout</param>
-      <param name="LinkedDofInfo">pressureLinkedDofs</param>
-    </struct>
-    <struct name="cellLayout">
-      <param name="Type">SingleCellLayout</param>
-    </struct>
-    <struct name="particleLayout">
-      <param name="Type">GaussParticleLayout</param>
-    </struct>
-    <struct name="gaussSwarm">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">particleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">gaussMapper</param>
-    </struct>
-    <struct name="gaussMapper">
-      <param name="Type">GaussMapper</param>
-      <param name="IntegrationPointsSwarm">gaussSwarm</param>
-      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
-    </struct>
-    <struct name="backgroundLayout">
-      <param name="Type">BackgroundParticleLayout</param>
-    </struct>
-    <struct name="gaussMaterialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">backgroundLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-    </struct>
-    <struct name="timeIntegrator">
-      <param name="Type">TimeIntegrator</param>
-      <param name="order">1</param>
-      <param name="simultaneous">t</param>
-      <param name="Context">context</param>
-    </struct>
-    <struct name="elementCellLayout">
-      <param name="Type">ElementCellLayout</param>
-      <param name="Mesh">mesh-linear</param>
-    </struct>
-    <struct name="discreteVoronoi">
-      <param name="Type">CellularAutomataVoronoi</param>
-      <param name="resolution">10</param>
-      <param name="diagonalNeighbours">true</param>
-    </struct>
-    <struct name="weights">
-      <param name="Type">DiscreteVoronoiWeights</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="localLayout">
-      <param name="Type">MappedParticleLayout</param>
-    </struct>
-    <struct name="picIntegrationPoints">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">localLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="WeightsCalculator">weights</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">mapper</param>
-    </struct>
-    <struct name="mapper">
-      <param name="Type">CoincidentMapper</param>
-      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="materialSwarmParticleLayout">
-      <param name="Type">RandomParticleLayout</param>
-      <param name="seed">13</param>
-      <param name="cellParticleCount">particlesPerCell</param>
-    </struct>
-    <struct name="materialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">materialSwarmParticleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="SplittingRoutine">splittingRoutine</param>
-      <param name="RemovalRoutine">removalRoutine</param>
-      <param name="EscapedRoutine">escapedRoutine</param>
-    </struct>
-    <struct name="materialSwarmAdvector">
-      <param name="Type">SwarmAdvector</param>
-      <param name="Swarm">materialSwarm</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="VelocityField">VelocityField</param>
-      <param name="PeriodicBCsManager">periodicBCsManager</param>
-      <param name="allowFallbackToFirstOrder">True</param>
-    </struct>
-    <struct name="removalRoutine">
-      <param name="Type">DiscreteVoronoiRemove</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="splittingRoutine">
-      <param name="Type">ReseedSplitting</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="solutionVelocity">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="solutionPressure">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">PressureField</param>
-    </struct>
-    <struct name="mom_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">VelocityField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="cont_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">PressureField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="k_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">VelocityField</param>
-      <param name="RHS">mom_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="constitutiveMatrix">
-      <param name="Type">ConstitutiveMatrixCartesian</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">k_matrix</param>
-    </struct>
-    <struct name="g_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">PressureField</param>
-      <param name="RHS">cont_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="gradientStiffnessMatrixTerm">
-      <param name="Type">GradientStiffnessMatrixTerm</param>
-      <param name="Swarm">gaussSwarm</param>
-      <param name="StiffnessMatrix">g_matrix</param>
-    </struct>
-    <struct name="preconditioner">
-      <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="preconditionerTerm">
-      <param name="Type">UzawaPreconditionerTerm</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">preconditioner</param>
-    </struct>
-    <struct name="uzawa">
-      <param name="Type">Stokes_SLE_UzawaSolver</param>
-      <param name="Preconditioner">preconditioner</param>
-      <param name="tolerance">1.0e-5</param>
-      <param name="maxIterations">5000</param>
-    </struct>
-    <struct name="stokesEqn">
-      <param name="Type">Stokes_SLE</param>
-      <param name="SLE_Solver">uzawa</param>
-      <param name="Context">context</param>
-      <param name="StressTensorMatrix">k_matrix</param>
-      <param name="GradientMatrix">g_matrix</param>
-      <param name="DivergenceMatrix"></param>
-      <param name="CompressibilityMatrix">c_matrix</param>
-      <param name="VelocityVector">solutionVelocity</param>
-      <param name="PressureVector">solutionPressure</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="ContinuityForceVector">cont_force</param>
-      <param name="killNonConvergent">false</param>
-      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">nonLinearTolerance</param>
-      <param name="makeConvergenceFile">false</param>
-    </struct>
-    <struct name="backgroundShape">
-      <param name="Type">Everywhere</param>
-    </struct>
-    <struct name="pdmsShape">
-      <param name="Type">Box</param>
-      <param name="startX">0.5</param>
-      <param name="endX">1.5</param>
-      <param name="startY">0.0</param>
-      <param name="endY">0.05</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-    <struct name="crustShape">
-      <param name="Type">Intersection</param>
-      <list name="shapes">
-        <param>backgroundShape</param>
-        <param>!pdmsShape</param>
-      </list>
-    </struct>
-    <struct name="crustViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">10.0</param>
-    </struct>
-    <struct name="strainWeakening">
-      <param name="Type">StrainWeakening</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="softeningStrain">0.1</param>
-      <param name="initialDamageFraction">0.0</param>
-      <param name="initialDamageWavenumber">0.5</param>
-      <param name="initialDamageFactor">0.5</param>
-      <param name="healingRate">0.0</param>
-    </struct>
-    <struct name="yielding">
-      <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">PressureField</param>
-      <param name="VelocityGradientsField">VelocityGradientsField</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="Context">context</param>
-      <param name="StrainWeakening">strainWeakening</param>
-      <param name="cohesion">10.0</param>
-      <param name="cohesionAfterSoftening">0.0001</param>
-      <param name="frictionCoefficient">0.0</param>
-      <param name="frictionCoefficientAfterSoftening">0.0</param>
-      <param name="minimumYieldStress">0.00001</param>
-    </struct>
-    <struct name="pdmsViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0</param>
-    </struct>
-    <struct name="storeViscosity">
-      <param name="Type">StoreVisc</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="crust">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">crustShape</param>
-      <list name="Rheology">
-        <param>crustViscosity</param>
-        <param>yielding</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="pdms">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">pdmsShape</param>
-      <list name="Rheology">
-        <param>pdmsViscosity</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="escapedRoutine">
-      <param name="Type">EscapedRoutine</param>
-      <param name="idealParticleCount">0</param>
-    </struct>
-    <struct name="velocityRemesher">
-      <param name="Type">StripRemesher</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="meshType">regular</param>
-      <list name="dims">
-        <param>true</param>
-        <param>true</param>
-        <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="dims">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
-    <struct name="conditionFunctions">
-      <param name="Type">StgFEM_StandardConditionFunctions</param>
-    </struct>
-  </struct>
-  <list name="plugins">
-    <param>Underworld_EulerDeform</param>
-    <param>Underworld_DumpSwarm</param>
-  </list>
-  <param name="maxTimeSteps">500</param>
-  <param name="dumpEvery">1</param>
-  <param name="outputPath">./output.extension</param>
-  <param name="dim">2</param>
-  <param name="shadowDepth">1</param>
-  <param name="minX">0.0f</param>
-  <param name="minY">0.0f</param>
-  <param name="minZ">0.0f</param>
-  <param name="maxX">2.0f</param>
-  <param name="maxY">0.35f</param>
-  <param name="maxZ">1.0f</param>
-  <param name="elementResI">16</param>
-  <param name="elementResJ">6</param>
-  <param name="elementResK">3</param>
-  <param name="allowUnbalancing">True</param>
-  <param name="buildNodeNeighbourTbl">True</param>
-  <param name="buildElementNeighbourTbl">True</param>
-  <param name="particlesPerCell">30</param>
-  <param name="seed">13</param>
-  <param name="gaussParticlesX">2</param>
-  <param name="gaussParticlesY">2</param>
-  <param name="gaussParticlesZ">2</param>
-  <param name="particlesPerCell">30</param>
-  <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>
-        <param name="wrapTop">True</param>
-        <list name="fields">
-          <struct>
-            <param name="field">VelocityField</param>
-            <param name="variable">velocity</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <struct name="velocityBCs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">front</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">back</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">left</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">right</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">1.0</param>
-          </struct>
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">bottom</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">func</param>
-            <param name="value">StepFunction</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <param name="StepFunctionLowerOffset">0.8</param>
-  <param name="StepFunctionUpperOffset">1.2</param>
-  <param name="StepFunctionValue">1.0</param>
-  <param name="StepFunctionDim">0</param>
-  <param name="StepFunctionLessThan">False</param>
-  <param name="checkpointEvery">1</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>

Deleted: long/3D/Gale/trunk/fountain.xml
===================================================================
--- long/3D/Gale/trunk/fountain.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/fountain.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -1,493 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-  <struct name="components">
-    <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>
-    </struct>
-    <struct name="linearMesh-generator">
-      <param name="Type">CartesianGenerator</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="dims">dim</param>
-      <param name="shadowDepth">shadowDepth</param>
-      <list name="size">
-        <param>elementResI</param>
-        <param>elementResJ</param>
-        <param>elementResK</param>
-      </list>
-      <list name="minCoord">
-        <param>minX</param>
-        <param>minY</param>
-        <param>minZ</param>
-      </list>
-      <list name="maxCoord">
-        <param>maxX</param>
-        <param>maxY</param>
-        <param>maxZ</param>
-      </list>
-    </struct>
-    <struct name="velocity">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Vector</param>
-      <param name="DataType">Double</param>
-      <param name="VectorComponentCount">dim</param>
-      <list name="names">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="velocityBCs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityICs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="BaseVariableCount">dim</param>
-      <list name="BaseVariables">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="VelocityField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">velocityDofLayout</param>
-      <param name="BC">velocityBCs</param>
-      <param name="IC">velocityICs</param>
-      <param name="LinkedDofInfo">velocityLinkedDofs</param>
-    </struct>
-    <struct name="VelocityMagnitudeField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Magnitude</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Gradient</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorInvariant</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorSymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="VorticityField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorAntisymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">SymmetricTensor_Invariant</param>
-      <param name="FeVariable">StrainRateField</param>
-    </struct>
-    <struct name="pressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</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>
-      <list name="BaseVariables">
-        <param>pressure</param>
-      </list>
-    </struct>
-    <struct name="PressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
-      <param name="DofLayout">pressureDofLayout</param>
-      <param name="LinkedDofInfo">pressureLinkedDofs</param>
-    </struct>
-    <struct name="cellLayout">
-      <param name="Type">SingleCellLayout</param>
-    </struct>
-    <struct name="particleLayout">
-      <param name="Type">GaussParticleLayout</param>
-    </struct>
-    <struct name="gaussSwarm">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">particleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">gaussMapper</param>
-    </struct>
-    <struct name="gaussMapper">
-      <param name="Type">GaussMapper</param>
-      <param name="IntegrationPointsSwarm">gaussSwarm</param>
-      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
-    </struct>
-    <struct name="backgroundLayout">
-      <param name="Type">BackgroundParticleLayout</param>
-    </struct>
-    <struct name="gaussMaterialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">backgroundLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-    </struct>
-    <struct name="timeIntegrator">
-      <param name="Type">TimeIntegrator</param>
-      <param name="order">1</param>
-      <param name="simultaneous">t</param>
-      <param name="Context">context</param>
-    </struct>
-    <struct name="elementCellLayout">
-      <param name="Type">ElementCellLayout</param>
-      <param name="Mesh">mesh-linear</param>
-    </struct>
-    <struct name="discreteVoronoi">
-      <param name="Type">CellularAutomataVoronoi</param>
-      <param name="resolution">10</param>
-      <param name="diagonalNeighbours">true</param>
-    </struct>
-    <struct name="weights">
-      <param name="Type">DiscreteVoronoiWeights</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="localLayout">
-      <param name="Type">MappedParticleLayout</param>
-    </struct>
-    <struct name="picIntegrationPoints">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">localLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="WeightsCalculator">weights</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">mapper</param>
-    </struct>
-    <struct name="mapper">
-      <param name="Type">CoincidentMapper</param>
-      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="materialSwarmParticleLayout">
-      <param name="Type">RandomParticleLayout</param>
-      <param name="seed">13</param>
-      <param name="cellParticleCount">particlesPerCell</param>
-    </struct>
-    <struct name="materialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">materialSwarmParticleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="SplittingRoutine">splittingRoutine</param>
-      <param name="RemovalRoutine">removalRoutine</param>
-      <param name="EscapedRoutine">escapedRoutine</param>
-    </struct>
-    <struct name="materialSwarmAdvector">
-      <param name="Type">SwarmAdvector</param>
-      <param name="Swarm">materialSwarm</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="VelocityField">VelocityField</param>
-      <param name="PeriodicBCsManager">periodicBCsManager</param>
-      <param name="allowFallbackToFirstOrder">True</param>
-    </struct>
-    <struct name="removalRoutine">
-      <param name="Type">DiscreteVoronoiRemove</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="splittingRoutine">
-      <param name="Type">ReseedSplitting</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="solutionVelocity">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="solutionPressure">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">PressureField</param>
-    </struct>
-    <struct name="mom_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">VelocityField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="cont_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">PressureField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="k_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">VelocityField</param>
-      <param name="RHS">mom_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="constitutiveMatrix">
-      <param name="Type">ConstitutiveMatrixCartesian</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">k_matrix</param>
-    </struct>
-    <struct name="g_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">PressureField</param>
-      <param name="RHS">cont_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="gradientStiffnessMatrixTerm">
-      <param name="Type">GradientStiffnessMatrixTerm</param>
-      <param name="Swarm">gaussSwarm</param>
-      <param name="StiffnessMatrix">g_matrix</param>
-    </struct>
-    <struct name="preconditioner">
-      <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="preconditionerTerm">
-      <param name="Type">UzawaPreconditionerTerm</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">preconditioner</param>
-    </struct>
-    <struct name="uzawa">
-      <param name="Type">Stokes_SLE_UzawaSolver</param>
-      <param name="Preconditioner">preconditioner</param>
-      <param name="tolerance">1.0e-5</param>
-      <param name="maxIterations">5000</param>
-    </struct>
-    <struct name="stokesEqn">
-      <param name="Type">Stokes_SLE</param>
-      <param name="SLE_Solver">uzawa</param>
-      <param name="Context">context</param>
-      <param name="StressTensorMatrix">k_matrix</param>
-      <param name="GradientMatrix">g_matrix</param>
-      <param name="DivergenceMatrix"></param>
-      <param name="CompressibilityMatrix">c_matrix</param>
-      <param name="VelocityVector">solutionVelocity</param>
-      <param name="PressureVector">solutionPressure</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="ContinuityForceVector">cont_force</param>
-      <param name="killNonConvergent">false</param>
-      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">nonLinearTolerance</param>
-      <param name="makeConvergenceFile">false</param>
-    </struct>
-<!--     <struct name="buoyancyForceTerm"> -->
-<!--       <param name="Type">BuoyancyForceTerm</param> -->
-<!--       <param name="ForceVector">mom_force</param> -->
-<!--       <param name="Swarm">picIntegrationPoints</param> -->
-<!--       <param name="gravity">gravity</param> -->
-<!--     </struct> -->
-    <struct name="conditionFunctions">
-      <param name="Type">StgFEM_StandardConditionFunctions</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="x_type">double</param> -->
-<!--       <param name="x_value">1.0</param> -->
-      <param name="y_type">func</param>
-      <param name="y_value">Gaussian</param>
-    </struct>
-    <struct name="background">
-      <param name="Type">Everywhere</param>
-    </struct>
-    <struct name="pdmsShape">
-      <param name="Type">Box</param>
-      <param name="startX">minX</param>
-      <param name="endX">maxX</param>
-      <param name="startY">minY</param>
-      <param name="endY">maxY</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-    <struct name="pdmsViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0</param>
-    </struct>
-    <struct name="storeViscosity">
-      <param name="Type">StoreVisc</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="pdms">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">pdmsShape</param>
-      <param name="density">1.0</param>
-      <list name="Rheology">
-        <param>pdmsViscosity</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="escapedRoutine">
-      <param name="Type">EscapedRoutine</param>
-      <param name="idealParticleCount">0</param>
-    </struct>
-    <struct name="velocityRemesher">
-      <param name="Type">StripRemesher</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="meshType">regular</param>
-      <list name="dims">
-        <param>true</param>
-        <param>true</param>
-        <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="dims">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
-  </struct>
-  <list name="plugins">
-    <param>Underworld_EulerDeform</param>
-  </list>
-  <param name="maxTimeSteps">500</param>
-  <param name="dumpEvery">1</param>
-  <param name="outputPath">./output.fountain</param>
-  <param name="dim">2</param>
-  <param name="shadowDepth">1</param>
-  <param name="minX">0.0f</param>
-  <param name="minY">0.0f</param>
-  <param name="minZ">0.0f</param>
-  <param name="maxX">1.0f</param>
-  <param name="maxY">0.5f</param>
-  <param name="maxZ">0.3f</param>
-  <param name="elementResI">32</param>
-  <param name="elementResJ">16</param>
-  <param name="elementResK">3</param>
-  <param name="allowUnbalancing">True</param>
-  <param name="buildNodeNeighbourTbl">True</param>
-  <param name="buildElementNeighbourTbl">True</param>
-  <param name="particlesPerCell">30</param>
-  <param name="seed">13</param>
-  <param name="gaussParticlesX">2</param>
-  <param name="gaussParticlesY">2</param>
-  <param name="gaussParticlesZ">2</param>
-  <param name="particlesPerCell">30</param>
-  <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>
-        <param name="wrapTop">True</param>
-        <param name="staticBottom">True</param>
-        <list name="fields">
-          <struct>
-            <param name="field">VelocityField</param>
-            <param name="variable">velocity</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <struct name="velocityBCs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">front</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">back</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">left</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">right</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-<!--       <struct> -->
-<!--         <param name="type">WallVC</param> -->
-<!--         <param name="wall">bottom</param> -->
-<!--         <list name="variables"> -->
-<!--           <struct> -->
-<!--             <param name="name">vy</param> -->
-<!--             <param name="type">double</param> -->
-<!--             <param name="value">0.0</param> -->
-<!--           </struct> -->
-<!--         </list> -->
-<!--       </struct> -->
-    </list>
-  </struct>
-  <param name="checkpointEvery">1</param>
-  <param name="gravity">1.</param>
-  <param name="GaussianHeight">1.0</param>
-  <param name="GaussianWidth">0.1</param>
-  <param name="GaussianCenter">0.5</param>
-  <param name="GaussianDim">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> -->
-</StGermainData>

Added: long/3D/Gale/trunk/input/extension.xml
===================================================================
--- long/3D/Gale/trunk/input/extension.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/input/extension.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -0,0 +1,535 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+  <struct name="components">
+    <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>
+    </struct>
+    <struct name="linearMesh-generator">
+      <param name="Type">CartesianGenerator</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="dims">dim</param>
+      <param name="shadowDepth">shadowDepth</param>
+      <list name="size">
+        <param>elementResI</param>
+        <param>elementResJ</param>
+        <param>elementResK</param>
+      </list>
+      <list name="minCoord">
+        <param>minX</param>
+        <param>minY</param>
+        <param>minZ</param>
+      </list>
+      <list name="maxCoord">
+        <param>maxX</param>
+        <param>maxY</param>
+        <param>maxZ</param>
+      </list>
+    </struct>
+    <struct name="velocity">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="Rank">Vector</param>
+      <param name="DataType">Double</param>
+      <param name="VectorComponentCount">dim</param>
+      <list name="names">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="velocityBCs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityICs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityDofLayout">
+      <param name="Type">DofLayout</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="BaseVariableCount">dim</param>
+      <list name="BaseVariables">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="VelocityField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-linear</param>
+      <param name="DofLayout">velocityDofLayout</param>
+      <param name="BC">velocityBCs</param>
+      <param name="IC">velocityICs</param>
+      <param name="LinkedDofInfo">velocityLinkedDofs</param>
+    </struct>
+    <struct name="VelocityMagnitudeField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Magnitude</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Gradient</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorInvariant</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorSymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="VorticityField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorAntisymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">SymmetricTensor_Invariant</param>
+      <param name="FeVariable">StrainRateField</param>
+    </struct>
+    <struct name="pressure">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-constant</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>
+      <list name="BaseVariables">
+        <param>pressure</param>
+      </list>
+    </struct>
+    <struct name="PressureField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-constant</param>
+      <param name="DofLayout">pressureDofLayout</param>
+      <param name="LinkedDofInfo">pressureLinkedDofs</param>
+    </struct>
+    <struct name="cellLayout">
+      <param name="Type">SingleCellLayout</param>
+    </struct>
+    <struct name="particleLayout">
+      <param name="Type">GaussParticleLayout</param>
+    </struct>
+    <struct name="gaussSwarm">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">particleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">gaussMapper</param>
+    </struct>
+    <struct name="gaussMapper">
+      <param name="Type">GaussMapper</param>
+      <param name="IntegrationPointsSwarm">gaussSwarm</param>
+      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+    </struct>
+    <struct name="backgroundLayout">
+      <param name="Type">BackgroundParticleLayout</param>
+    </struct>
+    <struct name="gaussMaterialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">backgroundLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+    </struct>
+    <struct name="timeIntegrator">
+      <param name="Type">TimeIntegrator</param>
+      <param name="order">1</param>
+      <param name="simultaneous">t</param>
+      <param name="Context">context</param>
+    </struct>
+    <struct name="elementCellLayout">
+      <param name="Type">ElementCellLayout</param>
+      <param name="Mesh">mesh-linear</param>
+    </struct>
+    <struct name="discreteVoronoi">
+      <param name="Type">CellularAutomataVoronoi</param>
+      <param name="resolution">10</param>
+      <param name="diagonalNeighbours">true</param>
+    </struct>
+    <struct name="weights">
+      <param name="Type">DiscreteVoronoiWeights</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="localLayout">
+      <param name="Type">MappedParticleLayout</param>
+    </struct>
+    <struct name="picIntegrationPoints">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">localLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="WeightsCalculator">weights</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">mapper</param>
+    </struct>
+    <struct name="mapper">
+      <param name="Type">CoincidentMapper</param>
+      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="materialSwarmParticleLayout">
+      <param name="Type">RandomParticleLayout</param>
+      <param name="seed">13</param>
+      <param name="cellParticleCount">particlesPerCell</param>
+    </struct>
+    <struct name="materialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">materialSwarmParticleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="SplittingRoutine">splittingRoutine</param>
+      <param name="RemovalRoutine">removalRoutine</param>
+      <param name="EscapedRoutine">escapedRoutine</param>
+    </struct>
+    <struct name="materialSwarmAdvector">
+      <param name="Type">SwarmAdvector</param>
+      <param name="Swarm">materialSwarm</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="VelocityField">VelocityField</param>
+      <param name="PeriodicBCsManager">periodicBCsManager</param>
+      <param name="allowFallbackToFirstOrder">True</param>
+    </struct>
+    <struct name="removalRoutine">
+      <param name="Type">DiscreteVoronoiRemove</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="splittingRoutine">
+      <param name="Type">ReseedSplitting</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="solutionVelocity">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="solutionPressure">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">PressureField</param>
+    </struct>
+    <struct name="mom_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">VelocityField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="cont_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">PressureField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="k_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">VelocityField</param>
+      <param name="RHS">mom_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="constitutiveMatrix">
+      <param name="Type">ConstitutiveMatrixCartesian</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">k_matrix</param>
+    </struct>
+    <struct name="g_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="gradientStiffnessMatrixTerm">
+      <param name="Type">GradientStiffnessMatrixTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="StiffnessMatrix">g_matrix</param>
+    </struct>
+    <struct name="preconditioner">
+      <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="preconditionerTerm">
+      <param name="Type">UzawaPreconditionerTerm</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">preconditioner</param>
+    </struct>
+    <struct name="uzawa">
+      <param name="Type">Stokes_SLE_UzawaSolver</param>
+      <param name="Preconditioner">preconditioner</param>
+      <param name="tolerance">1.0e-5</param>
+      <param name="maxIterations">5000</param>
+    </struct>
+    <struct name="stokesEqn">
+      <param name="Type">Stokes_SLE</param>
+      <param name="SLE_Solver">uzawa</param>
+      <param name="Context">context</param>
+      <param name="StressTensorMatrix">k_matrix</param>
+      <param name="GradientMatrix">g_matrix</param>
+      <param name="DivergenceMatrix"></param>
+      <param name="CompressibilityMatrix">c_matrix</param>
+      <param name="VelocityVector">solutionVelocity</param>
+      <param name="PressureVector">solutionPressure</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="ContinuityForceVector">cont_force</param>
+      <param name="killNonConvergent">false</param>
+      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+      <param name="nonLinearTolerance">nonLinearTolerance</param>
+      <param name="makeConvergenceFile">false</param>
+    </struct>
+    <struct name="backgroundShape">
+      <param name="Type">Everywhere</param>
+    </struct>
+    <struct name="pdmsShape">
+      <param name="Type">Box</param>
+      <param name="startX">0.5</param>
+      <param name="endX">1.5</param>
+      <param name="startY">0.0</param>
+      <param name="endY">0.05</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+    <struct name="crustShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>backgroundShape</param>
+        <param>!pdmsShape</param>
+      </list>
+    </struct>
+    <struct name="crustViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">10.0</param>
+    </struct>
+    <struct name="strainWeakening">
+      <param name="Type">StrainWeakening</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="softeningStrain">0.1</param>
+      <param name="initialDamageFraction">0.0</param>
+      <param name="initialDamageWavenumber">0.5</param>
+      <param name="initialDamageFactor">0.5</param>
+      <param name="healingRate">0.0</param>
+    </struct>
+    <struct name="yielding">
+      <param name="Type">MohrCoulomb</param>
+      <param name="PressureField">PressureField</param>
+      <param name="VelocityGradientsField">VelocityGradientsField</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="Context">context</param>
+      <param name="StrainWeakening">strainWeakening</param>
+      <param name="cohesion">10.0</param>
+      <param name="cohesionAfterSoftening">0.0001</param>
+      <param name="frictionCoefficient">0.0</param>
+      <param name="frictionCoefficientAfterSoftening">0.0</param>
+      <param name="minimumYieldStress">0.00001</param>
+    </struct>
+    <struct name="pdmsViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">1.0</param>
+    </struct>
+    <struct name="storeViscosity">
+      <param name="Type">StoreVisc</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="crust">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">crustShape</param>
+      <list name="Rheology">
+        <param>crustViscosity</param>
+        <param>yielding</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="pdms">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">pdmsShape</param>
+      <list name="Rheology">
+        <param>pdmsViscosity</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="escapedRoutine">
+      <param name="Type">EscapedRoutine</param>
+      <param name="idealParticleCount">0</param>
+    </struct>
+    <struct name="velocityRemesher">
+      <param name="Type">StripRemesher</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="meshType">regular</param>
+      <list name="dims">
+        <param>true</param>
+        <param>true</param>
+        <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="dims">3</param>
+      <param name="cellMesh">mesh-linear</param>
+    </struct>
+    <struct name="conditionFunctions">
+      <param name="Type">StgFEM_StandardConditionFunctions</param>
+    </struct>
+  </struct>
+  <list name="plugins">
+    <param>Underworld_EulerDeform</param>
+    <param>Underworld_DumpSwarm</param>
+  </list>
+  <param name="maxTimeSteps">500</param>
+  <param name="dumpEvery">1</param>
+  <param name="outputPath">./output.extension</param>
+  <param name="dim">2</param>
+  <param name="shadowDepth">1</param>
+  <param name="minX">0.0f</param>
+  <param name="minY">0.0f</param>
+  <param name="minZ">0.0f</param>
+  <param name="maxX">2.0f</param>
+  <param name="maxY">0.35f</param>
+  <param name="maxZ">1.0f</param>
+  <param name="elementResI">16</param>
+  <param name="elementResJ">6</param>
+  <param name="elementResK">3</param>
+  <param name="allowUnbalancing">True</param>
+  <param name="buildNodeNeighbourTbl">True</param>
+  <param name="buildElementNeighbourTbl">True</param>
+  <param name="particlesPerCell">30</param>
+  <param name="seed">13</param>
+  <param name="gaussParticlesX">2</param>
+  <param name="gaussParticlesY">2</param>
+  <param name="gaussParticlesZ">2</param>
+  <param name="particlesPerCell">30</param>
+  <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>
+        <param name="wrapTop">True</param>
+        <list name="fields">
+          <struct>
+            <param name="field">VelocityField</param>
+            <param name="variable">velocity</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <struct name="velocityBCs">
+    <param name="type">CompositeVC</param>
+    <list name="vcList">
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">front</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">back</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">left</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">right</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">1.0</param>
+          </struct>
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">bottom</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">func</param>
+            <param name="value">StepFunction</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <param name="StepFunctionLowerOffset">0.8</param>
+  <param name="StepFunctionUpperOffset">1.2</param>
+  <param name="StepFunctionValue">1.0</param>
+  <param name="StepFunctionDim">0</param>
+  <param name="StepFunctionLessThan">False</param>
+  <param name="checkpointEvery">1</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>

Added: long/3D/Gale/trunk/input/fountain.xml
===================================================================
--- long/3D/Gale/trunk/input/fountain.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/input/fountain.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -0,0 +1,493 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+  <struct name="components">
+    <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>
+    </struct>
+    <struct name="linearMesh-generator">
+      <param name="Type">CartesianGenerator</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="dims">dim</param>
+      <param name="shadowDepth">shadowDepth</param>
+      <list name="size">
+        <param>elementResI</param>
+        <param>elementResJ</param>
+        <param>elementResK</param>
+      </list>
+      <list name="minCoord">
+        <param>minX</param>
+        <param>minY</param>
+        <param>minZ</param>
+      </list>
+      <list name="maxCoord">
+        <param>maxX</param>
+        <param>maxY</param>
+        <param>maxZ</param>
+      </list>
+    </struct>
+    <struct name="velocity">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="Rank">Vector</param>
+      <param name="DataType">Double</param>
+      <param name="VectorComponentCount">dim</param>
+      <list name="names">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="velocityBCs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityICs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityDofLayout">
+      <param name="Type">DofLayout</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="BaseVariableCount">dim</param>
+      <list name="BaseVariables">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="VelocityField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-linear</param>
+      <param name="DofLayout">velocityDofLayout</param>
+      <param name="BC">velocityBCs</param>
+      <param name="IC">velocityICs</param>
+      <param name="LinkedDofInfo">velocityLinkedDofs</param>
+    </struct>
+    <struct name="VelocityMagnitudeField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Magnitude</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Gradient</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorInvariant</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorSymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="VorticityField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorAntisymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">SymmetricTensor_Invariant</param>
+      <param name="FeVariable">StrainRateField</param>
+    </struct>
+    <struct name="pressure">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-constant</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>
+      <list name="BaseVariables">
+        <param>pressure</param>
+      </list>
+    </struct>
+    <struct name="PressureField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-constant</param>
+      <param name="DofLayout">pressureDofLayout</param>
+      <param name="LinkedDofInfo">pressureLinkedDofs</param>
+    </struct>
+    <struct name="cellLayout">
+      <param name="Type">SingleCellLayout</param>
+    </struct>
+    <struct name="particleLayout">
+      <param name="Type">GaussParticleLayout</param>
+    </struct>
+    <struct name="gaussSwarm">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">particleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">gaussMapper</param>
+    </struct>
+    <struct name="gaussMapper">
+      <param name="Type">GaussMapper</param>
+      <param name="IntegrationPointsSwarm">gaussSwarm</param>
+      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+    </struct>
+    <struct name="backgroundLayout">
+      <param name="Type">BackgroundParticleLayout</param>
+    </struct>
+    <struct name="gaussMaterialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">backgroundLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+    </struct>
+    <struct name="timeIntegrator">
+      <param name="Type">TimeIntegrator</param>
+      <param name="order">1</param>
+      <param name="simultaneous">t</param>
+      <param name="Context">context</param>
+    </struct>
+    <struct name="elementCellLayout">
+      <param name="Type">ElementCellLayout</param>
+      <param name="Mesh">mesh-linear</param>
+    </struct>
+    <struct name="discreteVoronoi">
+      <param name="Type">CellularAutomataVoronoi</param>
+      <param name="resolution">10</param>
+      <param name="diagonalNeighbours">true</param>
+    </struct>
+    <struct name="weights">
+      <param name="Type">DiscreteVoronoiWeights</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="localLayout">
+      <param name="Type">MappedParticleLayout</param>
+    </struct>
+    <struct name="picIntegrationPoints">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">localLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="WeightsCalculator">weights</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">mapper</param>
+    </struct>
+    <struct name="mapper">
+      <param name="Type">CoincidentMapper</param>
+      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="materialSwarmParticleLayout">
+      <param name="Type">RandomParticleLayout</param>
+      <param name="seed">13</param>
+      <param name="cellParticleCount">particlesPerCell</param>
+    </struct>
+    <struct name="materialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">materialSwarmParticleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="SplittingRoutine">splittingRoutine</param>
+      <param name="RemovalRoutine">removalRoutine</param>
+      <param name="EscapedRoutine">escapedRoutine</param>
+    </struct>
+    <struct name="materialSwarmAdvector">
+      <param name="Type">SwarmAdvector</param>
+      <param name="Swarm">materialSwarm</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="VelocityField">VelocityField</param>
+      <param name="PeriodicBCsManager">periodicBCsManager</param>
+      <param name="allowFallbackToFirstOrder">True</param>
+    </struct>
+    <struct name="removalRoutine">
+      <param name="Type">DiscreteVoronoiRemove</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="splittingRoutine">
+      <param name="Type">ReseedSplitting</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="solutionVelocity">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="solutionPressure">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">PressureField</param>
+    </struct>
+    <struct name="mom_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">VelocityField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="cont_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">PressureField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="k_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">VelocityField</param>
+      <param name="RHS">mom_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="constitutiveMatrix">
+      <param name="Type">ConstitutiveMatrixCartesian</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">k_matrix</param>
+    </struct>
+    <struct name="g_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="gradientStiffnessMatrixTerm">
+      <param name="Type">GradientStiffnessMatrixTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="StiffnessMatrix">g_matrix</param>
+    </struct>
+    <struct name="preconditioner">
+      <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="preconditionerTerm">
+      <param name="Type">UzawaPreconditionerTerm</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">preconditioner</param>
+    </struct>
+    <struct name="uzawa">
+      <param name="Type">Stokes_SLE_UzawaSolver</param>
+      <param name="Preconditioner">preconditioner</param>
+      <param name="tolerance">1.0e-5</param>
+      <param name="maxIterations">5000</param>
+    </struct>
+    <struct name="stokesEqn">
+      <param name="Type">Stokes_SLE</param>
+      <param name="SLE_Solver">uzawa</param>
+      <param name="Context">context</param>
+      <param name="StressTensorMatrix">k_matrix</param>
+      <param name="GradientMatrix">g_matrix</param>
+      <param name="DivergenceMatrix"></param>
+      <param name="CompressibilityMatrix">c_matrix</param>
+      <param name="VelocityVector">solutionVelocity</param>
+      <param name="PressureVector">solutionPressure</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="ContinuityForceVector">cont_force</param>
+      <param name="killNonConvergent">false</param>
+      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+      <param name="nonLinearTolerance">nonLinearTolerance</param>
+      <param name="makeConvergenceFile">false</param>
+    </struct>
+<!--     <struct name="buoyancyForceTerm"> -->
+<!--       <param name="Type">BuoyancyForceTerm</param> -->
+<!--       <param name="ForceVector">mom_force</param> -->
+<!--       <param name="Swarm">picIntegrationPoints</param> -->
+<!--       <param name="gravity">gravity</param> -->
+<!--     </struct> -->
+    <struct name="conditionFunctions">
+      <param name="Type">StgFEM_StandardConditionFunctions</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="x_type">double</param> -->
+<!--       <param name="x_value">1.0</param> -->
+      <param name="y_type">func</param>
+      <param name="y_value">Gaussian</param>
+    </struct>
+    <struct name="background">
+      <param name="Type">Everywhere</param>
+    </struct>
+    <struct name="pdmsShape">
+      <param name="Type">Box</param>
+      <param name="startX">minX</param>
+      <param name="endX">maxX</param>
+      <param name="startY">minY</param>
+      <param name="endY">maxY</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+    <struct name="pdmsViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">1.0</param>
+    </struct>
+    <struct name="storeViscosity">
+      <param name="Type">StoreVisc</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="pdms">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">pdmsShape</param>
+      <param name="density">1.0</param>
+      <list name="Rheology">
+        <param>pdmsViscosity</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="escapedRoutine">
+      <param name="Type">EscapedRoutine</param>
+      <param name="idealParticleCount">0</param>
+    </struct>
+    <struct name="velocityRemesher">
+      <param name="Type">StripRemesher</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="meshType">regular</param>
+      <list name="dims">
+        <param>true</param>
+        <param>true</param>
+        <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="dims">3</param>
+      <param name="cellMesh">mesh-linear</param>
+    </struct>
+  </struct>
+  <list name="plugins">
+    <param>Underworld_EulerDeform</param>
+  </list>
+  <param name="maxTimeSteps">500</param>
+  <param name="dumpEvery">1</param>
+  <param name="outputPath">./output.fountain</param>
+  <param name="dim">2</param>
+  <param name="shadowDepth">1</param>
+  <param name="minX">0.0f</param>
+  <param name="minY">0.0f</param>
+  <param name="minZ">0.0f</param>
+  <param name="maxX">1.0f</param>
+  <param name="maxY">0.5f</param>
+  <param name="maxZ">0.3f</param>
+  <param name="elementResI">32</param>
+  <param name="elementResJ">16</param>
+  <param name="elementResK">3</param>
+  <param name="allowUnbalancing">True</param>
+  <param name="buildNodeNeighbourTbl">True</param>
+  <param name="buildElementNeighbourTbl">True</param>
+  <param name="particlesPerCell">30</param>
+  <param name="seed">13</param>
+  <param name="gaussParticlesX">2</param>
+  <param name="gaussParticlesY">2</param>
+  <param name="gaussParticlesZ">2</param>
+  <param name="particlesPerCell">30</param>
+  <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>
+        <param name="wrapTop">True</param>
+        <param name="staticBottom">True</param>
+        <list name="fields">
+          <struct>
+            <param name="field">VelocityField</param>
+            <param name="variable">velocity</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <struct name="velocityBCs">
+    <param name="type">CompositeVC</param>
+    <list name="vcList">
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">front</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">back</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">left</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">right</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+<!--       <struct> -->
+<!--         <param name="type">WallVC</param> -->
+<!--         <param name="wall">bottom</param> -->
+<!--         <list name="variables"> -->
+<!--           <struct> -->
+<!--             <param name="name">vy</param> -->
+<!--             <param name="type">double</param> -->
+<!--             <param name="value">0.0</param> -->
+<!--           </struct> -->
+<!--         </list> -->
+<!--       </struct> -->
+    </list>
+  </struct>
+  <param name="checkpointEvery">1</param>
+  <param name="gravity">1.</param>
+  <param name="GaussianHeight">1.0</param>
+  <param name="GaussianWidth">0.1</param>
+  <param name="GaussianCenter">0.5</param>
+  <param name="GaussianDim">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> -->
+</StGermainData>

Added: long/3D/Gale/trunk/input/shortening.xml
===================================================================
--- long/3D/Gale/trunk/input/shortening.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/input/shortening.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -0,0 +1,521 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+  <struct name="components">
+    <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>
+    </struct>
+    <struct name="linearMesh-generator">
+      <param name="Type">CartesianGenerator</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="dims">dim</param>
+      <param name="shadowDepth">shadowDepth</param>
+      <list name="size">
+        <param>elementResI</param>
+        <param>elementResJ</param>
+        <param>elementResK</param>
+      </list>
+      <list name="minCoord">
+        <param>minX</param>
+        <param>minY</param>
+        <param>minZ</param>
+      </list>
+      <list name="maxCoord">
+        <param>maxX</param>
+        <param>maxY</param>
+        <param>maxZ</param>
+      </list>
+    </struct>
+    <struct name="velocity">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="Rank">Vector</param>
+      <param name="DataType">Double</param>
+      <param name="VectorComponentCount">dim</param>
+      <list name="names">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="velocityBCs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityICs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityDofLayout">
+      <param name="Type">DofLayout</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="BaseVariableCount">dim</param>
+      <list name="BaseVariables">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="VelocityField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-linear</param>
+      <param name="DofLayout">velocityDofLayout</param>
+      <param name="BC">velocityBCs</param>
+      <param name="IC">velocityICs</param>
+      <param name="LinkedDofInfo">velocityLinkedDofs</param>
+    </struct>
+    <struct name="VelocityMagnitudeField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Magnitude</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Gradient</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorInvariant</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorSymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="VorticityField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorAntisymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">SymmetricTensor_Invariant</param>
+      <param name="FeVariable">StrainRateField</param>
+    </struct>
+    <struct name="pressure">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-constant</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>
+      <list name="BaseVariables">
+        <param>pressure</param>
+      </list>
+    </struct>
+    <struct name="PressureField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-constant</param>
+      <param name="DofLayout">pressureDofLayout</param>
+      <param name="LinkedDofInfo">pressureLinkedDofs</param>
+    </struct>
+    <struct name="cellLayout">
+      <param name="Type">SingleCellLayout</param>
+    </struct>
+    <struct name="particleLayout">
+      <param name="Type">GaussParticleLayout</param>
+    </struct>
+    <struct name="gaussSwarm">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">particleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">gaussMapper</param>
+    </struct>
+    <struct name="gaussMapper">
+      <param name="Type">GaussMapper</param>
+      <param name="IntegrationPointsSwarm">gaussSwarm</param>
+      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+    </struct>
+    <struct name="backgroundLayout">
+      <param name="Type">BackgroundParticleLayout</param>
+    </struct>
+    <struct name="gaussMaterialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">backgroundLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+    </struct>
+    <struct name="timeIntegrator">
+      <param name="Type">TimeIntegrator</param>
+      <param name="order">1</param>
+      <param name="simultaneous">t</param>
+      <param name="Context">context</param>
+    </struct>
+    <struct name="elementCellLayout">
+      <param name="Type">ElementCellLayout</param>
+      <param name="Mesh">mesh-linear</param>
+    </struct>
+    <struct name="discreteVoronoi">
+      <param name="Type">CellularAutomataVoronoi</param>
+      <param name="resolution">10</param>
+      <param name="diagonalNeighbours">true</param>
+    </struct>
+    <struct name="weights">
+      <param name="Type">DiscreteVoronoiWeights</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="localLayout">
+      <param name="Type">MappedParticleLayout</param>
+    </struct>
+    <struct name="picIntegrationPoints">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">localLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="WeightsCalculator">weights</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">mapper</param>
+    </struct>
+    <struct name="mapper">
+      <param name="Type">CoincidentMapper</param>
+      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="materialSwarmParticleLayout">
+      <param name="Type">RandomParticleLayout</param>
+      <param name="seed">13</param>
+      <param name="cellParticleCount">particlesPerCell</param>
+    </struct>
+    <struct name="materialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">materialSwarmParticleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="SplittingRoutine">splittingRoutine</param>
+      <param name="RemovalRoutine">removalRoutine</param>
+      <param name="EscapedRoutine">escapedRoutine</param>
+    </struct>
+    <struct name="materialSwarmAdvector">
+      <param name="Type">SwarmAdvector</param>
+      <param name="Swarm">materialSwarm</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="VelocityField">VelocityField</param>
+      <param name="PeriodicBCsManager">periodicBCsManager</param>
+      <param name="allowFallbackToFirstOrder">True</param>
+    </struct>
+    <struct name="removalRoutine">
+      <param name="Type">DiscreteVoronoiRemove</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="splittingRoutine">
+      <param name="Type">ReseedSplitting</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="solutionVelocity">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="solutionPressure">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">PressureField</param>
+    </struct>
+    <struct name="mom_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">VelocityField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="cont_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">PressureField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="k_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">VelocityField</param>
+      <param name="RHS">mom_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="constitutiveMatrix">
+      <param name="Type">ConstitutiveMatrixCartesian</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">k_matrix</param>
+    </struct>
+    <struct name="g_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="gradientStiffnessMatrixTerm">
+      <param name="Type">GradientStiffnessMatrixTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="StiffnessMatrix">g_matrix</param>
+    </struct>
+    <struct name="preconditioner">
+      <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="preconditionerTerm">
+      <param name="Type">UzawaPreconditionerTerm</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">preconditioner</param>
+    </struct>
+    <struct name="uzawa">
+      <param name="Type">Stokes_SLE_UzawaSolver</param>
+      <param name="Preconditioner">preconditioner</param>
+      <param name="tolerance">1.0e-5</param>
+      <param name="maxIterations">5000</param>
+    </struct>
+    <struct name="stokesEqn">
+      <param name="Type">Stokes_SLE</param>
+      <param name="SLE_Solver">uzawa</param>
+      <param name="Context">context</param>
+      <param name="StressTensorMatrix">k_matrix</param>
+      <param name="GradientMatrix">g_matrix</param>
+      <param name="DivergenceMatrix"></param>
+      <param name="CompressibilityMatrix">c_matrix</param>
+      <param name="VelocityVector">solutionVelocity</param>
+      <param name="PressureVector">solutionPressure</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="ContinuityForceVector">cont_force</param>
+      <param name="killNonConvergent">false</param>
+      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+      <param name="nonLinearTolerance">nonLinearTolerance</param>
+      <param name="makeConvergenceFile">false</param>
+    </struct>
+    <struct name="backgroundShape">
+      <param name="Type">Everywhere</param>
+    </struct>
+    <struct name="beadsShape">
+      <param name="Type">Box</param>
+      <param name="startX">minX</param>
+      <param name="endX">maxX</param>
+      <param name="startY">0.05</param>
+      <param name="endY">0.1</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+    <struct name="crustShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>backgroundShape</param>
+        <param>!beadsShape</param>
+      </list>
+    </struct>
+    <struct name="crustViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">10.0</param>
+    </struct>
+    <struct name="strainWeakening">
+      <param name="Type">StrainWeakening</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="softeningStrain">0.1</param>
+      <param name="initialDamageFraction">0.0</param>
+      <param name="initialDamageWavenumber">0.5</param>
+      <param name="initialDamageFactor">0.5</param>
+      <param name="healingRate">0.0</param>
+    </struct>
+    <struct name="yielding">
+      <param name="Type">MohrCoulomb</param>
+      <param name="PressureField">PressureField</param>
+      <param name="VelocityGradientsField">VelocityGradientsField</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="Context">context</param>
+      <param name="StrainWeakening">strainWeakening</param>
+      <param name="cohesion">10.0</param>
+      <param name="cohesionAfterSoftening">0.0001</param>
+      <param name="frictionCoefficient">0.0</param>
+      <param name="frictionCoefficientAfterSoftening">0.0</param>
+      <param name="minimumYieldStress">0.00001</param>
+    </struct>
+    <struct name="beadsViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">1.0</param>
+    </struct>
+    <struct name="storeViscosity">
+      <param name="Type">StoreVisc</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="crust">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">crustShape</param>
+      <list name="Rheology">
+        <param>crustViscosity</param>
+        <param>yielding</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="beads">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">beadsShape</param>
+      <list name="Rheology">
+        <param>beadsViscosity</param>
+        <param>storeViscosity</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">wedge</param>
+      <list name="origin">
+        <param>0.0</param>
+      </list>
+      <param name="offset">1.5</param>
+      <param name="gradient">0.1763269807</param>
+    </struct>
+    <struct name="escapedRoutine">
+      <param name="Type">EscapedRoutine</param>
+      <param name="idealParticleCount">0</param>
+    </struct>
+    <struct name="velocityRemesher">
+      <param name="Type">StripRemesher</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="meshType">regular</param>
+      <list name="dims">
+        <param>true</param>
+        <param>true</param>
+        <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="dims">3</param>
+      <param name="cellMesh">mesh-linear</param>
+    </struct>
+    <struct name="conditionFunctions">
+      <param name="Type">StgFEM_StandardConditionFunctions</param>
+    </struct>
+  </struct>
+  <list name="plugins">
+    <param>Underworld_EulerDeform</param>
+    <param>Underworld_DumpSwarm</param>
+  </list>
+  <param name="maxTimeSteps">500</param>
+  <param name="outputEvery">1</param>
+  <param name="dumpEvery">1</param>
+  <param name="outputPath">./output.shortening</param>
+  <param name="dim">2</param>
+  <param name="shadowDepth">1</param>
+  <param name="minX">0.0f</param>
+  <param name="minY">0.0f</param>
+  <param name="minZ">0.0f</param>
+  <param name="maxX">2.0f</param>
+  <param name="maxY">0.35f</param>
+  <param name="maxZ">1.0f</param>
+  <param name="elementResI">32</param>
+  <param name="elementResJ">16</param>
+  <param name="elementResK">4</param>
+  <param name="allowUnbalancing">True</param>
+  <param name="buildNodeNeighbourTbl">True</param>
+  <param name="buildElementNeighbourTbl">True</param>
+  <param name="particlesPerCell">30</param>
+  <param name="seed">13</param>
+  <param name="gaussParticlesX">2</param>
+  <param name="gaussParticlesY">2</param>
+  <param name="gaussParticlesZ">2</param>
+  <param name="particlesPerCell">30</param>
+  <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>
+        <param name="wrapTop">True</param>
+        <list name="fields">
+          <struct>
+            <param name="field">VelocityField</param>
+            <param name="variable">velocity</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <struct name="velocityBCs">
+    <param name="type">CompositeVC</param>
+    <list name="vcList">
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">front</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">back</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">left</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">right</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">-1.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">bottom</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <param name="checkpointEvery">1</param>
+</StGermainData>

Added: long/3D/Gale/trunk/input/sinusoid.xml
===================================================================
--- long/3D/Gale/trunk/input/sinusoid.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/input/sinusoid.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -0,0 +1,486 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+  <struct name="components">
+    <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>
+    </struct>
+    <struct name="linearMesh-generator">
+      <param name="Type">CartesianGenerator</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="dims">dim</param>
+      <param name="shadowDepth">shadowDepth</param>
+      <list name="size">
+        <param>elementResI</param>
+        <param>elementResJ</param>
+        <param>elementResK</param>
+      </list>
+      <list name="minCoord">
+        <param>minX</param>
+        <param>minY</param>
+        <param>minZ</param>
+      </list>
+      <list name="maxCoord">
+        <param>maxX</param>
+        <param>maxY</param>
+        <param>maxZ</param>
+      </list>
+    </struct>
+    <struct name="velocity">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="Rank">Vector</param>
+      <param name="DataType">Double</param>
+      <param name="VectorComponentCount">dim</param>
+      <list name="names">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="velocityBCs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityICs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityDofLayout">
+      <param name="Type">DofLayout</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="BaseVariableCount">dim</param>
+      <list name="BaseVariables">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="VelocityField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-linear</param>
+      <param name="DofLayout">velocityDofLayout</param>
+      <param name="BC">velocityBCs</param>
+      <param name="IC">velocityICs</param>
+      <param name="LinkedDofInfo">velocityLinkedDofs</param>
+    </struct>
+    <struct name="VelocityMagnitudeField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Magnitude</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Gradient</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorInvariant</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorSymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="VorticityField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorAntisymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">SymmetricTensor_Invariant</param>
+      <param name="FeVariable">StrainRateField</param>
+    </struct>
+    <struct name="pressure">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-constant</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>
+      <list name="BaseVariables">
+        <param>pressure</param>
+      </list>
+    </struct>
+    <struct name="PressureField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-constant</param>
+      <param name="DofLayout">pressureDofLayout</param>
+      <param name="LinkedDofInfo">pressureLinkedDofs</param>
+    </struct>
+    <struct name="cellLayout">
+      <param name="Type">SingleCellLayout</param>
+    </struct>
+    <struct name="particleLayout">
+      <param name="Type">GaussParticleLayout</param>
+    </struct>
+    <struct name="gaussSwarm">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">particleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">gaussMapper</param>
+    </struct>
+    <struct name="gaussMapper">
+      <param name="Type">GaussMapper</param>
+      <param name="IntegrationPointsSwarm">gaussSwarm</param>
+      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+    </struct>
+    <struct name="backgroundLayout">
+      <param name="Type">BackgroundParticleLayout</param>
+    </struct>
+    <struct name="gaussMaterialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">backgroundLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+    </struct>
+    <struct name="timeIntegrator">
+      <param name="Type">TimeIntegrator</param>
+      <param name="order">1</param>
+      <param name="simultaneous">t</param>
+      <param name="Context">context</param>
+    </struct>
+    <struct name="elementCellLayout">
+      <param name="Type">ElementCellLayout</param>
+      <param name="Mesh">mesh-linear</param>
+    </struct>
+    <struct name="discreteVoronoi">
+      <param name="Type">CellularAutomataVoronoi</param>
+      <param name="resolution">10</param>
+      <param name="diagonalNeighbours">true</param>
+    </struct>
+    <struct name="weights">
+      <param name="Type">DiscreteVoronoiWeights</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="localLayout">
+      <param name="Type">MappedParticleLayout</param>
+    </struct>
+    <struct name="picIntegrationPoints">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">localLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="WeightsCalculator">weights</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">mapper</param>
+    </struct>
+    <struct name="mapper">
+      <param name="Type">CoincidentMapper</param>
+      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="materialSwarmParticleLayout">
+      <param name="Type">RandomParticleLayout</param>
+      <param name="seed">13</param>
+      <param name="cellParticleCount">particlesPerCell</param>
+    </struct>
+    <struct name="materialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">materialSwarmParticleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="SplittingRoutine">splittingRoutine</param>
+      <param name="RemovalRoutine">removalRoutine</param>
+      <param name="EscapedRoutine">escapedRoutine</param>
+    </struct>
+    <struct name="materialSwarmAdvector">
+      <param name="Type">SwarmAdvector</param>
+      <param name="Swarm">materialSwarm</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="VelocityField">VelocityField</param>
+      <param name="PeriodicBCsManager">periodicBCsManager</param>
+      <param name="allowFallbackToFirstOrder">True</param>
+    </struct>
+    <struct name="removalRoutine">
+      <param name="Type">DiscreteVoronoiRemove</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="splittingRoutine">
+      <param name="Type">ReseedSplitting</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="solutionVelocity">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="solutionPressure">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">PressureField</param>
+    </struct>
+    <struct name="mom_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">VelocityField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="cont_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">PressureField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="k_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">VelocityField</param>
+      <param name="RHS">mom_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="constitutiveMatrix">
+      <param name="Type">ConstitutiveMatrixCartesian</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">k_matrix</param>
+    </struct>
+    <struct name="g_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="gradientStiffnessMatrixTerm">
+      <param name="Type">GradientStiffnessMatrixTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="StiffnessMatrix">g_matrix</param>
+    </struct>
+    <struct name="preconditioner">
+      <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="preconditionerTerm">
+      <param name="Type">UzawaPreconditionerTerm</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">preconditioner</param>
+    </struct>
+    <struct name="uzawa">
+      <param name="Type">Stokes_SLE_UzawaSolver</param>
+      <param name="Preconditioner">preconditioner</param>
+      <param name="tolerance">1.0e-5</param>
+      <param name="maxIterations">5000</param>
+    </struct>
+    <struct name="stokesEqn">
+      <param name="Type">Stokes_SLE</param>
+      <param name="SLE_Solver">uzawa</param>
+      <param name="Context">context</param>
+      <param name="StressTensorMatrix">k_matrix</param>
+      <param name="GradientMatrix">g_matrix</param>
+      <param name="DivergenceMatrix"></param>
+      <param name="CompressibilityMatrix">c_matrix</param>
+      <param name="VelocityVector">solutionVelocity</param>
+      <param name="PressureVector">solutionPressure</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="ContinuityForceVector">cont_force</param>
+      <param name="killNonConvergent">false</param>
+      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+      <param name="nonLinearTolerance">nonLinearTolerance</param>
+      <param name="makeConvergenceFile">false</param>
+    </struct>
+    <struct name="buoyancyForceTerm">
+      <param name="Type">BuoyancyForceTerm</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="gravity">gravity</param>
+    </struct>
+    <struct name="background">
+      <param name="Type">Everywhere</param>
+    </struct>
+    <struct name="pdmsShape">
+      <param name="Type">Box</param>
+      <param name="startX">minX</param>
+      <param name="endX">maxX</param>
+      <param name="startY">minY</param>
+      <param name="endY">maxY</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+    <struct name="pdmsViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">1.0</param>
+    </struct>
+    <struct name="storeViscosity">
+      <param name="Type">StoreVisc</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="pdms">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">pdmsShape</param>
+      <param name="density">1.0</param>
+      <list name="Rheology">
+        <param>pdmsViscosity</param>
+        <param>storeViscosity</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">cosine</param>
+      <list name="origin">
+        <param>0.0</param>
+      </list>
+      <param name="amplitude">0.1</param>
+      <param name="frequency">6.28318530718</param>
+    </struct>
+    <struct name="escapedRoutine">
+      <param name="Type">EscapedRoutine</param>
+      <param name="idealParticleCount">0</param>
+    </struct>
+    <struct name="velocityRemesher">
+      <param name="Type">StripRemesher</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="meshType">regular</param>
+      <list name="dims">
+        <param>true</param>
+        <param>true</param>
+        <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="dims">3</param>
+      <param name="cellMesh">mesh-linear</param>
+    </struct>
+    <struct name="conditionFunctions">
+      <param name="Type">StgFEM_StandardConditionFunctions</param>
+    </struct>
+  </struct>
+  <list name="plugins">
+    <param>Underworld_EulerDeform</param>
+    <param>Underworld_DumpSwarm</param>
+  </list>
+  <param name="maxTimeSteps">500</param>
+  <param name="outputEvery">1</param>
+  <param name="dumpEvery">1</param>
+  <param name="outputPath">./output.sinusoid</param>
+  <param name="dim">2</param>
+  <param name="shadowDepth">1</param>
+  <param name="minX">0.0f</param>
+  <param name="minY">0.0f</param>
+  <param name="minZ">0.0f</param>
+  <param name="maxX">1.0f</param>
+  <param name="maxY">0.5f</param>
+  <param name="maxZ">0.3f</param>
+  <param name="elementResI">32</param>
+  <param name="elementResJ">16</param>
+  <param name="elementResK">3</param>
+  <param name="allowUnbalancing">True</param>
+  <param name="buildNodeNeighbourTbl">True</param>
+  <param name="buildElementNeighbourTbl">True</param>
+  <param name="particlesPerCell">30</param>
+  <param name="seed">13</param>
+  <param name="gaussParticlesX">2</param>
+  <param name="gaussParticlesY">2</param>
+  <param name="gaussParticlesZ">2</param>
+  <param name="particlesPerCell">30</param>
+  <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>
+        <param name="wrapTop">True</param>
+        <list name="fields">
+          <struct>
+            <param name="field">VelocityField</param>
+            <param name="variable">velocity</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <struct name="velocityBCs">
+    <param name="type">CompositeVC</param>
+    <list name="vcList">
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">front</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">back</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">left</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">right</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">bottom</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <param name="checkpointEvery">1</param>
+  <param name="gravity">1.</param>
+</StGermainData>

Added: long/3D/Gale/trunk/input/subduction.xml
===================================================================
--- long/3D/Gale/trunk/input/subduction.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/input/subduction.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -0,0 +1,595 @@
+<?xml version="1.0"?>
+<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
+  <struct name="components">
+    <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>
+    </struct>
+    <struct name="linearMesh-generator">
+      <param name="Type">CartesianGenerator</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="dims">dim</param>
+      <param name="shadowDepth">shadowDepth</param>
+      <list name="size">
+        <param>elementResI</param>
+        <param>elementResJ</param>
+        <param>elementResK</param>
+      </list>
+      <list name="minCoord">
+        <param>minX</param>
+        <param>minY</param>
+        <param>minZ</param>
+      </list>
+      <list name="maxCoord">
+        <param>maxX</param>
+        <param>maxY</param>
+        <param>maxZ</param>
+      </list>
+    </struct>
+    <struct name="velocity">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="Rank">Vector</param>
+      <param name="DataType">Double</param>
+      <param name="VectorComponentCount">dim</param>
+      <list name="names">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="velocityBCs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityICs">
+      <param name="Type">CompositeVC</param>
+      <param name="Data">mesh-linear</param>
+    </struct>
+    <struct name="velocityDofLayout">
+      <param name="Type">DofLayout</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="BaseVariableCount">dim</param>
+      <list name="BaseVariables">
+        <param>vx</param>
+        <param>vy</param>
+        <param>vz</param>
+      </list>
+    </struct>
+    <struct name="VelocityField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-linear</param>
+      <param name="DofLayout">velocityDofLayout</param>
+      <param name="BC">velocityBCs</param>
+      <param name="IC">velocityICs</param>
+      <param name="LinkedDofInfo">velocityLinkedDofs</param>
+    </struct>
+    <struct name="VelocityMagnitudeField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Magnitude</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">Gradient</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="VelocityGradientsInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorInvariant</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorSymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="VorticityField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">TensorAntisymmetricPart</param>
+      <param name="FeVariable">VelocityGradientsField</param>
+    </struct>
+    <struct name="StrainRateInvariantField">
+      <param name="Type">OperatorFeVariable</param>
+      <param name="Operator">SymmetricTensor_Invariant</param>
+      <param name="FeVariable">StrainRateField</param>
+    </struct>
+    <struct name="pressure">
+      <param name="Type">MeshVariable</param>
+      <param name="mesh">mesh-constant</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>
+      <list name="BaseVariables">
+        <param>pressure</param>
+      </list>
+    </struct>
+    <struct name="PressureField">
+      <param name="Type">FeVariable</param>
+      <param name="FEMesh">mesh-constant</param>
+      <param name="DofLayout">pressureDofLayout</param>
+      <param name="LinkedDofInfo">pressureLinkedDofs</param>
+    </struct>
+    <struct name="cellLayout">
+      <param name="Type">SingleCellLayout</param>
+    </struct>
+    <struct name="particleLayout">
+      <param name="Type">GaussParticleLayout</param>
+    </struct>
+    <struct name="gaussSwarm">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">particleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">gaussMapper</param>
+    </struct>
+    <struct name="gaussMapper">
+      <param name="Type">GaussMapper</param>
+      <param name="IntegrationPointsSwarm">gaussSwarm</param>
+      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
+    </struct>
+    <struct name="backgroundLayout">
+      <param name="Type">BackgroundParticleLayout</param>
+    </struct>
+    <struct name="gaussMaterialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">cellLayout</param>
+      <param name="ParticleLayout">backgroundLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+    </struct>
+    <struct name="timeIntegrator">
+      <param name="Type">TimeIntegrator</param>
+      <param name="order">1</param>
+      <param name="simultaneous">t</param>
+      <param name="Context">context</param>
+    </struct>
+    <struct name="elementCellLayout">
+      <param name="Type">ElementCellLayout</param>
+      <param name="Mesh">mesh-linear</param>
+    </struct>
+    <struct name="discreteVoronoi">
+      <param name="Type">CellularAutomataVoronoi</param>
+      <param name="resolution">10</param>
+      <param name="diagonalNeighbours">true</param>
+    </struct>
+    <struct name="weights">
+      <param name="Type">DiscreteVoronoiWeights</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="localLayout">
+      <param name="Type">MappedParticleLayout</param>
+    </struct>
+    <struct name="picIntegrationPoints">
+      <param name="Type">IntegrationPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">localLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="WeightsCalculator">weights</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="IntegrationPointMapper">mapper</param>
+    </struct>
+    <struct name="mapper">
+      <param name="Type">CoincidentMapper</param>
+      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="materialSwarmParticleLayout">
+      <param name="Type">RandomParticleLayout</param>
+      <param name="seed">13</param>
+      <param name="cellParticleCount">particlesPerCell</param>
+    </struct>
+    <struct name="materialSwarm">
+      <param name="Type">MaterialPointsSwarm</param>
+      <param name="CellLayout">elementCellLayout</param>
+      <param name="ParticleLayout">materialSwarmParticleLayout</param>
+      <param name="FeMesh">mesh-linear</param>
+      <param name="SplittingRoutine">splittingRoutine</param>
+      <param name="RemovalRoutine">removalRoutine</param>
+      <param name="EscapedRoutine">escapedRoutine</param>
+    </struct>
+    <struct name="materialSwarmAdvector">
+      <param name="Type">SwarmAdvector</param>
+      <param name="Swarm">materialSwarm</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="VelocityField">VelocityField</param>
+      <param name="PeriodicBCsManager">periodicBCsManager</param>
+      <param name="allowFallbackToFirstOrder">True</param>
+    </struct>
+    <struct name="removalRoutine">
+      <param name="Type">DiscreteVoronoiRemove</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="splittingRoutine">
+      <param name="Type">ReseedSplitting</param>
+      <param name="idealParticleCount">particlesPerCell</param>
+      <param name="originalParticleShape">leftBoundaryShape</param>
+      <param name="DiscreteVoronoi">discreteVoronoi</param>
+    </struct>
+    <struct name="solutionVelocity">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">VelocityField</param>
+    </struct>
+    <struct name="solutionPressure">
+      <param name="Type">SolutionVector</param>
+      <param name="FeVariable">PressureField</param>
+    </struct>
+    <struct name="mom_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">VelocityField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="cont_force">
+      <param name="Type">ForceVector</param>
+      <param name="FeVariable">PressureField</param>
+      <param name="ExtraInfo">context</param>
+    </struct>
+    <struct name="k_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">VelocityField</param>
+      <param name="RHS">mom_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="constitutiveMatrix">
+      <param name="Type">ConstitutiveMatrixCartesian</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">k_matrix</param>
+    </struct>
+    <struct name="g_matrix">
+      <param name="Type">StiffnessMatrix</param>
+      <param name="RowVariable">VelocityField</param>
+      <param name="ColumnVariable">PressureField</param>
+      <param name="RHS">cont_force</param>
+      <param name="allowZeroElementContributions">False</param>
+    </struct>
+    <struct name="gradientStiffnessMatrixTerm">
+      <param name="Type">GradientStiffnessMatrixTerm</param>
+      <param name="Swarm">gaussSwarm</param>
+      <param name="StiffnessMatrix">g_matrix</param>
+    </struct>
+    <struct name="preconditioner">
+      <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="preconditionerTerm">
+      <param name="Type">UzawaPreconditionerTerm</param>
+      <param name="Swarm">picIntegrationPoints</param>
+      <param name="StiffnessMatrix">preconditioner</param>
+    </struct>
+    <struct name="uzawa">
+      <param name="Type">Stokes_SLE_UzawaSolver</param>
+      <param name="Preconditioner">preconditioner</param>
+      <param name="tolerance">1.0e-5</param>
+      <param name="maxIterations">5000</param>
+    </struct>
+    <struct name="stokesEqn">
+      <param name="Type">Stokes_SLE</param>
+      <param name="SLE_Solver">uzawa</param>
+      <param name="Context">context</param>
+      <param name="StressTensorMatrix">k_matrix</param>
+      <param name="GradientMatrix">g_matrix</param>
+      <param name="DivergenceMatrix"></param>
+      <param name="CompressibilityMatrix">c_matrix</param>
+      <param name="VelocityVector">solutionVelocity</param>
+      <param name="PressureVector">solutionPressure</param>
+      <param name="ForceVector">mom_force</param>
+      <param name="ContinuityForceVector">cont_force</param>
+      <param name="killNonConvergent">false</param>
+      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
+      <param name="nonLinearTolerance">nonLinearTolerance</param>
+      <param name="makeConvergenceFile">false</param>
+    </struct>
+    <struct name="backgroundShape">
+      <param name="Type">Everywhere</param>
+    </struct>
+    <struct name="leftBoundaryShape">
+      <param name="Type">Box</param>
+      <param name="startX">minX</param>
+      <param name="endX">0.1</param>
+      <param name="startY">minY</param>
+      <param name="endY">maxY</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+    <struct name="rightBoundaryShape">
+      <param name="Type">Box</param>
+      <param name="startX">1.9</param>
+      <param name="endX">maxX</param>
+      <param name="startY">minY</param>
+      <param name="endY">maxY</param>
+      <param name="startZ">minZ</param>
+      <param name="endZ">maxZ</param>
+    </struct>
+
+    <struct name="leftSlabOuterSphere">
+      <param name="Type">Sphere</param>
+      <param name="CentreX">0.0</param>
+      <param name="CentreY">-7.15</param>
+      <param name="radius">7.25</param>
+    </struct>
+    <struct name="leftSlabInnerSphere">
+      <param name="Type">Sphere</param>
+      <param name="CentreX">0.0</param>
+      <param name="CentreY">-6.15</param>
+      <param name="radius">6.2</param>
+    </struct>
+    <struct name="leftSlabShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>leftSlabOuterSphere</param>
+        <param>!leftSlabInnerSphere</param>
+      </list>
+    </struct>
+
+
+
+    <struct name="rightSlabOuterSphere">
+      <param name="Type">Sphere</param>
+      <param name="CentreX">maxX</param>
+      <param name="CentreY">-7.15</param>
+      <param name="radius">7.25</param>
+    </struct>
+    <struct name="rightSlabInnerSphere">
+      <param name="Type">Sphere</param>
+      <param name="CentreX">maxX</param>
+      <param name="CentreY">-6.15</param>
+      <param name="radius">6.2</param>
+    </struct>
+    <struct name="rightSlabShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+              <param>rightSlabOuterSphere</param>
+        <param>!rightSlabInnerSphere</param>
+        <param>!leftSlabShape</param>
+      </list>
+    </struct>
+    
+    <struct name="beadsShape">
+      <param name="Type">Intersection</param>
+      <list name="shapes">
+        <param>backgroundShape</param>
+        <param>!leftSlabShape</param>
+        <param>!rightSlabShape</param>
+      </list>
+    </struct>
+    <struct name="crustShape">
+      <param name="Type">Union</param>
+      <list name="shapes">
+        <param>leftSlabShape</param>
+        <param>rightSlabShape</param>
+      </list>
+    </struct>
+    <struct name="crustViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">10.0</param>
+    </struct>
+    <struct name="strainWeakening">
+      <param name="Type">StrainWeakening</param>
+      <param name="TimeIntegrator">timeIntegrator</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="softeningStrain">0.1</param>
+      <param name="initialDamageFraction">0.0</param>
+      <param name="initialDamageWavenumber">0.5</param>
+      <param name="initialDamageFactor">0.5</param>
+      <param name="healingRate">0.0</param>
+    </struct>
+    <struct name="yielding">
+      <param name="Type">MohrCoulomb</param>
+      <param name="PressureField">PressureField</param>
+      <param name="VelocityGradientsField">VelocityGradientsField</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+      <param name="Context">context</param>
+      <param name="StrainWeakening">strainWeakening</param>
+      <param name="cohesion">100.0</param>
+      <param name="cohesionAfterSoftening">0.0001</param>
+      <param name="frictionCoefficient">0.0</param>
+      <param name="frictionCoefficientAfterSoftening">0.0</param>
+      <param name="minimumYieldStress">0.00001</param>
+    </struct>
+    <struct name="beadsViscosity">
+      <param name="Type">MaterialViscosity</param>
+      <param name="eta0">0.1</param>
+    </struct>
+    <struct name="storeViscosity">
+      <param name="Type">StoreVisc</param>
+      <param name="MaterialPointsSwarm">materialSwarm</param>
+    </struct>
+    <struct name="crust">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">crustShape</param>
+      <list name="Rheology">
+        <param>crustViscosity</param>
+        <param>yielding</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="beads">
+      <param name="Type">RheologyMaterial</param>
+      <param name="Shape">beadsShape</param>
+      <list name="Rheology">
+        <param>beadsViscosity</param>
+        <param>storeViscosity</param>
+      </list>
+    </struct>
+    <struct name="escapedRoutine">
+      <param name="Type">EscapedRoutine</param>
+      <param name="idealParticleCount">0</param>
+    </struct>
+    <struct name="velocityRemesher">
+      <param name="Type">StripRemesher</param>
+      <param name="mesh">mesh-linear</param>
+      <param name="meshType">regular</param>
+      <list name="dims">
+        <param>true</param>
+        <param>true</param>
+        <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="dims">3</param>
+      <param name="cellMesh">mesh-linear</param>
+    </struct>
+    <struct name="conditionFunctions">
+      <param name="Type">StgFEM_StandardConditionFunctions</param>
+    </struct>
+  </struct>
+  <list name="plugins">
+    <param>Underworld_EulerDeform</param>
+    <param>Underworld_DumpSwarm</param>
+    <param>Gale_SurfaceProcess</param>
+  </list>
+  <param name="maxTimeSteps">500</param>
+  <param name="outputEvery">1</param>
+  <param name="dumpEvery">1</param>
+  <param name="outputPath">./output.subduction</param>
+  <param name="dim">2</param>
+  <param name="shadowDepth">1</param>
+  <param name="minX">0.0f</param>
+  <param name="minY">0.0f</param>
+  <param name="minZ">0.0f</param>
+  <param name="maxX">2.0f</param>
+  <param name="maxY">0.2f</param>
+  <param name="maxZ">1.0f</param>
+  <param name="elementResI">32</param>
+  <param name="elementResJ">16</param>
+  <param name="elementResK">4</param>
+  <param name="allowUnbalancing">True</param>
+  <param name="buildNodeNeighbourTbl">True</param>
+  <param name="buildElementNeighbourTbl">True</param>
+  <param name="particlesPerCell">30</param>
+  <param name="seed">13</param>
+  <param name="gaussParticlesX">2</param>
+  <param name="gaussParticlesY">2</param>
+  <param name="gaussParticlesZ">2</param>
+  <param name="particlesPerCell">30</param>
+  <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>
+        <param name="wrapTop">True</param>
+        <param name="staticLeft">True</param>
+        <param name="staticBottom">True</param>
+        <list name="fields">
+          <struct>
+            <param name="field">VelocityField</param>
+            <param name="variable">velocity</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <struct name="velocityBCs">
+    <param name="type">CompositeVC</param>
+    <list name="vcList">
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">front</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">back</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vz</param>
+            <param name="type">double</param>
+            <param name="value">0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">left</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">func</param>
+            <param name="value">StepFunctionProduct3</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">right</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">double</param>
+            <param name="value">0.0</param>
+          </struct>
+        </list>
+      </struct>
+      <struct>
+        <param name="type">WallVC</param>
+        <param name="wall">bottom</param>
+        <list name="variables">
+          <struct>
+            <param name="name">vy</param>
+            <param name="type">func</param>
+            <param name="value">StepFunctionProduct2</param>
+          </struct>
+          <struct>
+            <param name="name">vx</param>
+            <param name="type">func</param>
+            <param name="value">StepFunctionProduct1</param>
+          </struct>
+        </list>
+      </struct>
+    </list>
+  </struct>
+  <param name="StepFunctionProduct1Start">0.8</param>
+  <param name="StepFunctionProduct1End">1.2</param>
+  <param name="StepFunctionProduct1Value">0.5</param>
+  <param name="StepFunctionProduct2Start">0.8</param>
+  <param name="StepFunctionProduct2End">1.2</param>
+  <param name="StepFunctionProduct2Value">-0.1</param>
+
+  <param name="StepFunctionProduct3Start">0.05</param>
+  <param name="StepFunctionProduct3End">0.1</param>
+  <param name="StepFunctionProduct3Value">0.5</param>
+
+  <struct name="SurfaceProcess">
+    <param name="mesh">mesh-linear</param>
+  </struct>
+
+  <param name="checkpointEvery">1</param>
+</StGermainData>

Deleted: long/3D/Gale/trunk/shortening.xml
===================================================================
--- long/3D/Gale/trunk/shortening.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/shortening.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -1,521 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-  <struct name="components">
-    <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>
-    </struct>
-    <struct name="linearMesh-generator">
-      <param name="Type">CartesianGenerator</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="dims">dim</param>
-      <param name="shadowDepth">shadowDepth</param>
-      <list name="size">
-        <param>elementResI</param>
-        <param>elementResJ</param>
-        <param>elementResK</param>
-      </list>
-      <list name="minCoord">
-        <param>minX</param>
-        <param>minY</param>
-        <param>minZ</param>
-      </list>
-      <list name="maxCoord">
-        <param>maxX</param>
-        <param>maxY</param>
-        <param>maxZ</param>
-      </list>
-    </struct>
-    <struct name="velocity">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Vector</param>
-      <param name="DataType">Double</param>
-      <param name="VectorComponentCount">dim</param>
-      <list name="names">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="velocityBCs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityICs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="BaseVariableCount">dim</param>
-      <list name="BaseVariables">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="VelocityField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">velocityDofLayout</param>
-      <param name="BC">velocityBCs</param>
-      <param name="IC">velocityICs</param>
-      <param name="LinkedDofInfo">velocityLinkedDofs</param>
-    </struct>
-    <struct name="VelocityMagnitudeField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Magnitude</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Gradient</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorInvariant</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorSymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="VorticityField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorAntisymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">SymmetricTensor_Invariant</param>
-      <param name="FeVariable">StrainRateField</param>
-    </struct>
-    <struct name="pressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</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>
-      <list name="BaseVariables">
-        <param>pressure</param>
-      </list>
-    </struct>
-    <struct name="PressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
-      <param name="DofLayout">pressureDofLayout</param>
-      <param name="LinkedDofInfo">pressureLinkedDofs</param>
-    </struct>
-    <struct name="cellLayout">
-      <param name="Type">SingleCellLayout</param>
-    </struct>
-    <struct name="particleLayout">
-      <param name="Type">GaussParticleLayout</param>
-    </struct>
-    <struct name="gaussSwarm">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">particleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">gaussMapper</param>
-    </struct>
-    <struct name="gaussMapper">
-      <param name="Type">GaussMapper</param>
-      <param name="IntegrationPointsSwarm">gaussSwarm</param>
-      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
-    </struct>
-    <struct name="backgroundLayout">
-      <param name="Type">BackgroundParticleLayout</param>
-    </struct>
-    <struct name="gaussMaterialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">backgroundLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-    </struct>
-    <struct name="timeIntegrator">
-      <param name="Type">TimeIntegrator</param>
-      <param name="order">1</param>
-      <param name="simultaneous">t</param>
-      <param name="Context">context</param>
-    </struct>
-    <struct name="elementCellLayout">
-      <param name="Type">ElementCellLayout</param>
-      <param name="Mesh">mesh-linear</param>
-    </struct>
-    <struct name="discreteVoronoi">
-      <param name="Type">CellularAutomataVoronoi</param>
-      <param name="resolution">10</param>
-      <param name="diagonalNeighbours">true</param>
-    </struct>
-    <struct name="weights">
-      <param name="Type">DiscreteVoronoiWeights</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="localLayout">
-      <param name="Type">MappedParticleLayout</param>
-    </struct>
-    <struct name="picIntegrationPoints">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">localLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="WeightsCalculator">weights</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">mapper</param>
-    </struct>
-    <struct name="mapper">
-      <param name="Type">CoincidentMapper</param>
-      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="materialSwarmParticleLayout">
-      <param name="Type">RandomParticleLayout</param>
-      <param name="seed">13</param>
-      <param name="cellParticleCount">particlesPerCell</param>
-    </struct>
-    <struct name="materialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">materialSwarmParticleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="SplittingRoutine">splittingRoutine</param>
-      <param name="RemovalRoutine">removalRoutine</param>
-      <param name="EscapedRoutine">escapedRoutine</param>
-    </struct>
-    <struct name="materialSwarmAdvector">
-      <param name="Type">SwarmAdvector</param>
-      <param name="Swarm">materialSwarm</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="VelocityField">VelocityField</param>
-      <param name="PeriodicBCsManager">periodicBCsManager</param>
-      <param name="allowFallbackToFirstOrder">True</param>
-    </struct>
-    <struct name="removalRoutine">
-      <param name="Type">DiscreteVoronoiRemove</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="splittingRoutine">
-      <param name="Type">ReseedSplitting</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="solutionVelocity">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="solutionPressure">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">PressureField</param>
-    </struct>
-    <struct name="mom_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">VelocityField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="cont_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">PressureField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="k_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">VelocityField</param>
-      <param name="RHS">mom_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="constitutiveMatrix">
-      <param name="Type">ConstitutiveMatrixCartesian</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">k_matrix</param>
-    </struct>
-    <struct name="g_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">PressureField</param>
-      <param name="RHS">cont_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="gradientStiffnessMatrixTerm">
-      <param name="Type">GradientStiffnessMatrixTerm</param>
-      <param name="Swarm">gaussSwarm</param>
-      <param name="StiffnessMatrix">g_matrix</param>
-    </struct>
-    <struct name="preconditioner">
-      <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="preconditionerTerm">
-      <param name="Type">UzawaPreconditionerTerm</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">preconditioner</param>
-    </struct>
-    <struct name="uzawa">
-      <param name="Type">Stokes_SLE_UzawaSolver</param>
-      <param name="Preconditioner">preconditioner</param>
-      <param name="tolerance">1.0e-5</param>
-      <param name="maxIterations">5000</param>
-    </struct>
-    <struct name="stokesEqn">
-      <param name="Type">Stokes_SLE</param>
-      <param name="SLE_Solver">uzawa</param>
-      <param name="Context">context</param>
-      <param name="StressTensorMatrix">k_matrix</param>
-      <param name="GradientMatrix">g_matrix</param>
-      <param name="DivergenceMatrix"></param>
-      <param name="CompressibilityMatrix">c_matrix</param>
-      <param name="VelocityVector">solutionVelocity</param>
-      <param name="PressureVector">solutionPressure</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="ContinuityForceVector">cont_force</param>
-      <param name="killNonConvergent">false</param>
-      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">nonLinearTolerance</param>
-      <param name="makeConvergenceFile">false</param>
-    </struct>
-    <struct name="backgroundShape">
-      <param name="Type">Everywhere</param>
-    </struct>
-    <struct name="beadsShape">
-      <param name="Type">Box</param>
-      <param name="startX">minX</param>
-      <param name="endX">maxX</param>
-      <param name="startY">0.05</param>
-      <param name="endY">0.1</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-    <struct name="crustShape">
-      <param name="Type">Intersection</param>
-      <list name="shapes">
-        <param>backgroundShape</param>
-        <param>!beadsShape</param>
-      </list>
-    </struct>
-    <struct name="crustViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">10.0</param>
-    </struct>
-    <struct name="strainWeakening">
-      <param name="Type">StrainWeakening</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="softeningStrain">0.1</param>
-      <param name="initialDamageFraction">0.0</param>
-      <param name="initialDamageWavenumber">0.5</param>
-      <param name="initialDamageFactor">0.5</param>
-      <param name="healingRate">0.0</param>
-    </struct>
-    <struct name="yielding">
-      <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">PressureField</param>
-      <param name="VelocityGradientsField">VelocityGradientsField</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="Context">context</param>
-      <param name="StrainWeakening">strainWeakening</param>
-      <param name="cohesion">10.0</param>
-      <param name="cohesionAfterSoftening">0.0001</param>
-      <param name="frictionCoefficient">0.0</param>
-      <param name="frictionCoefficientAfterSoftening">0.0</param>
-      <param name="minimumYieldStress">0.00001</param>
-    </struct>
-    <struct name="beadsViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0</param>
-    </struct>
-    <struct name="storeViscosity">
-      <param name="Type">StoreVisc</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="crust">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">crustShape</param>
-      <list name="Rheology">
-        <param>crustViscosity</param>
-        <param>yielding</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="beads">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">beadsShape</param>
-      <list name="Rheology">
-        <param>beadsViscosity</param>
-        <param>storeViscosity</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">wedge</param>
-      <list name="origin">
-        <param>0.0</param>
-      </list>
-      <param name="offset">1.5</param>
-      <param name="gradient">0.1763269807</param>
-    </struct>
-    <struct name="escapedRoutine">
-      <param name="Type">EscapedRoutine</param>
-      <param name="idealParticleCount">0</param>
-    </struct>
-    <struct name="velocityRemesher">
-      <param name="Type">StripRemesher</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="meshType">regular</param>
-      <list name="dims">
-        <param>true</param>
-        <param>true</param>
-        <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="dims">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
-    <struct name="conditionFunctions">
-      <param name="Type">StgFEM_StandardConditionFunctions</param>
-    </struct>
-  </struct>
-  <list name="plugins">
-    <param>Underworld_EulerDeform</param>
-    <param>Underworld_DumpSwarm</param>
-  </list>
-  <param name="maxTimeSteps">500</param>
-  <param name="outputEvery">1</param>
-  <param name="dumpEvery">1</param>
-  <param name="outputPath">./output.shortening</param>
-  <param name="dim">2</param>
-  <param name="shadowDepth">1</param>
-  <param name="minX">0.0f</param>
-  <param name="minY">0.0f</param>
-  <param name="minZ">0.0f</param>
-  <param name="maxX">2.0f</param>
-  <param name="maxY">0.35f</param>
-  <param name="maxZ">1.0f</param>
-  <param name="elementResI">32</param>
-  <param name="elementResJ">16</param>
-  <param name="elementResK">4</param>
-  <param name="allowUnbalancing">True</param>
-  <param name="buildNodeNeighbourTbl">True</param>
-  <param name="buildElementNeighbourTbl">True</param>
-  <param name="particlesPerCell">30</param>
-  <param name="seed">13</param>
-  <param name="gaussParticlesX">2</param>
-  <param name="gaussParticlesY">2</param>
-  <param name="gaussParticlesZ">2</param>
-  <param name="particlesPerCell">30</param>
-  <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>
-        <param name="wrapTop">True</param>
-        <list name="fields">
-          <struct>
-            <param name="field">VelocityField</param>
-            <param name="variable">velocity</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <struct name="velocityBCs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">front</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">back</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">left</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">right</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">-1.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">bottom</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <param name="checkpointEvery">1</param>
-</StGermainData>

Deleted: long/3D/Gale/trunk/sinusoid.xml
===================================================================
--- long/3D/Gale/trunk/sinusoid.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/sinusoid.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -1,486 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-  <struct name="components">
-    <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>
-    </struct>
-    <struct name="linearMesh-generator">
-      <param name="Type">CartesianGenerator</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="dims">dim</param>
-      <param name="shadowDepth">shadowDepth</param>
-      <list name="size">
-        <param>elementResI</param>
-        <param>elementResJ</param>
-        <param>elementResK</param>
-      </list>
-      <list name="minCoord">
-        <param>minX</param>
-        <param>minY</param>
-        <param>minZ</param>
-      </list>
-      <list name="maxCoord">
-        <param>maxX</param>
-        <param>maxY</param>
-        <param>maxZ</param>
-      </list>
-    </struct>
-    <struct name="velocity">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Vector</param>
-      <param name="DataType">Double</param>
-      <param name="VectorComponentCount">dim</param>
-      <list name="names">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="velocityBCs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityICs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="BaseVariableCount">dim</param>
-      <list name="BaseVariables">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="VelocityField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">velocityDofLayout</param>
-      <param name="BC">velocityBCs</param>
-      <param name="IC">velocityICs</param>
-      <param name="LinkedDofInfo">velocityLinkedDofs</param>
-    </struct>
-    <struct name="VelocityMagnitudeField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Magnitude</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Gradient</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorInvariant</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorSymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="VorticityField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorAntisymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">SymmetricTensor_Invariant</param>
-      <param name="FeVariable">StrainRateField</param>
-    </struct>
-    <struct name="pressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</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>
-      <list name="BaseVariables">
-        <param>pressure</param>
-      </list>
-    </struct>
-    <struct name="PressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
-      <param name="DofLayout">pressureDofLayout</param>
-      <param name="LinkedDofInfo">pressureLinkedDofs</param>
-    </struct>
-    <struct name="cellLayout">
-      <param name="Type">SingleCellLayout</param>
-    </struct>
-    <struct name="particleLayout">
-      <param name="Type">GaussParticleLayout</param>
-    </struct>
-    <struct name="gaussSwarm">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">particleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">gaussMapper</param>
-    </struct>
-    <struct name="gaussMapper">
-      <param name="Type">GaussMapper</param>
-      <param name="IntegrationPointsSwarm">gaussSwarm</param>
-      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
-    </struct>
-    <struct name="backgroundLayout">
-      <param name="Type">BackgroundParticleLayout</param>
-    </struct>
-    <struct name="gaussMaterialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">backgroundLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-    </struct>
-    <struct name="timeIntegrator">
-      <param name="Type">TimeIntegrator</param>
-      <param name="order">1</param>
-      <param name="simultaneous">t</param>
-      <param name="Context">context</param>
-    </struct>
-    <struct name="elementCellLayout">
-      <param name="Type">ElementCellLayout</param>
-      <param name="Mesh">mesh-linear</param>
-    </struct>
-    <struct name="discreteVoronoi">
-      <param name="Type">CellularAutomataVoronoi</param>
-      <param name="resolution">10</param>
-      <param name="diagonalNeighbours">true</param>
-    </struct>
-    <struct name="weights">
-      <param name="Type">DiscreteVoronoiWeights</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="localLayout">
-      <param name="Type">MappedParticleLayout</param>
-    </struct>
-    <struct name="picIntegrationPoints">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">localLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="WeightsCalculator">weights</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">mapper</param>
-    </struct>
-    <struct name="mapper">
-      <param name="Type">CoincidentMapper</param>
-      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="materialSwarmParticleLayout">
-      <param name="Type">RandomParticleLayout</param>
-      <param name="seed">13</param>
-      <param name="cellParticleCount">particlesPerCell</param>
-    </struct>
-    <struct name="materialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">materialSwarmParticleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="SplittingRoutine">splittingRoutine</param>
-      <param name="RemovalRoutine">removalRoutine</param>
-      <param name="EscapedRoutine">escapedRoutine</param>
-    </struct>
-    <struct name="materialSwarmAdvector">
-      <param name="Type">SwarmAdvector</param>
-      <param name="Swarm">materialSwarm</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="VelocityField">VelocityField</param>
-      <param name="PeriodicBCsManager">periodicBCsManager</param>
-      <param name="allowFallbackToFirstOrder">True</param>
-    </struct>
-    <struct name="removalRoutine">
-      <param name="Type">DiscreteVoronoiRemove</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="splittingRoutine">
-      <param name="Type">ReseedSplitting</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="solutionVelocity">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="solutionPressure">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">PressureField</param>
-    </struct>
-    <struct name="mom_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">VelocityField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="cont_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">PressureField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="k_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">VelocityField</param>
-      <param name="RHS">mom_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="constitutiveMatrix">
-      <param name="Type">ConstitutiveMatrixCartesian</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">k_matrix</param>
-    </struct>
-    <struct name="g_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">PressureField</param>
-      <param name="RHS">cont_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="gradientStiffnessMatrixTerm">
-      <param name="Type">GradientStiffnessMatrixTerm</param>
-      <param name="Swarm">gaussSwarm</param>
-      <param name="StiffnessMatrix">g_matrix</param>
-    </struct>
-    <struct name="preconditioner">
-      <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="preconditionerTerm">
-      <param name="Type">UzawaPreconditionerTerm</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">preconditioner</param>
-    </struct>
-    <struct name="uzawa">
-      <param name="Type">Stokes_SLE_UzawaSolver</param>
-      <param name="Preconditioner">preconditioner</param>
-      <param name="tolerance">1.0e-5</param>
-      <param name="maxIterations">5000</param>
-    </struct>
-    <struct name="stokesEqn">
-      <param name="Type">Stokes_SLE</param>
-      <param name="SLE_Solver">uzawa</param>
-      <param name="Context">context</param>
-      <param name="StressTensorMatrix">k_matrix</param>
-      <param name="GradientMatrix">g_matrix</param>
-      <param name="DivergenceMatrix"></param>
-      <param name="CompressibilityMatrix">c_matrix</param>
-      <param name="VelocityVector">solutionVelocity</param>
-      <param name="PressureVector">solutionPressure</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="ContinuityForceVector">cont_force</param>
-      <param name="killNonConvergent">false</param>
-      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">nonLinearTolerance</param>
-      <param name="makeConvergenceFile">false</param>
-    </struct>
-    <struct name="buoyancyForceTerm">
-      <param name="Type">BuoyancyForceTerm</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="gravity">gravity</param>
-    </struct>
-    <struct name="background">
-      <param name="Type">Everywhere</param>
-    </struct>
-    <struct name="pdmsShape">
-      <param name="Type">Box</param>
-      <param name="startX">minX</param>
-      <param name="endX">maxX</param>
-      <param name="startY">minY</param>
-      <param name="endY">maxY</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-    <struct name="pdmsViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">1.0</param>
-    </struct>
-    <struct name="storeViscosity">
-      <param name="Type">StoreVisc</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="pdms">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">pdmsShape</param>
-      <param name="density">1.0</param>
-      <list name="Rheology">
-        <param>pdmsViscosity</param>
-        <param>storeViscosity</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">cosine</param>
-      <list name="origin">
-        <param>0.0</param>
-      </list>
-      <param name="amplitude">0.1</param>
-      <param name="frequency">6.28318530718</param>
-    </struct>
-    <struct name="escapedRoutine">
-      <param name="Type">EscapedRoutine</param>
-      <param name="idealParticleCount">0</param>
-    </struct>
-    <struct name="velocityRemesher">
-      <param name="Type">StripRemesher</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="meshType">regular</param>
-      <list name="dims">
-        <param>true</param>
-        <param>true</param>
-        <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="dims">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
-    <struct name="conditionFunctions">
-      <param name="Type">StgFEM_StandardConditionFunctions</param>
-    </struct>
-  </struct>
-  <list name="plugins">
-    <param>Underworld_EulerDeform</param>
-    <param>Underworld_DumpSwarm</param>
-  </list>
-  <param name="maxTimeSteps">500</param>
-  <param name="outputEvery">1</param>
-  <param name="dumpEvery">1</param>
-  <param name="outputPath">./output.sinusoid</param>
-  <param name="dim">2</param>
-  <param name="shadowDepth">1</param>
-  <param name="minX">0.0f</param>
-  <param name="minY">0.0f</param>
-  <param name="minZ">0.0f</param>
-  <param name="maxX">1.0f</param>
-  <param name="maxY">0.5f</param>
-  <param name="maxZ">0.3f</param>
-  <param name="elementResI">32</param>
-  <param name="elementResJ">16</param>
-  <param name="elementResK">3</param>
-  <param name="allowUnbalancing">True</param>
-  <param name="buildNodeNeighbourTbl">True</param>
-  <param name="buildElementNeighbourTbl">True</param>
-  <param name="particlesPerCell">30</param>
-  <param name="seed">13</param>
-  <param name="gaussParticlesX">2</param>
-  <param name="gaussParticlesY">2</param>
-  <param name="gaussParticlesZ">2</param>
-  <param name="particlesPerCell">30</param>
-  <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>
-        <param name="wrapTop">True</param>
-        <list name="fields">
-          <struct>
-            <param name="field">VelocityField</param>
-            <param name="variable">velocity</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <struct name="velocityBCs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">front</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">back</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">left</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">right</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">bottom</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <param name="checkpointEvery">1</param>
-  <param name="gravity">1.</param>
-</StGermainData>

Deleted: long/3D/Gale/trunk/subduction.xml
===================================================================
--- long/3D/Gale/trunk/subduction.xml	2007-01-18 00:35:59 UTC (rev 5817)
+++ long/3D/Gale/trunk/subduction.xml	2007-01-18 00:36:03 UTC (rev 5818)
@@ -1,595 +0,0 @@
-<?xml version="1.0"?>
-<StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003">
-  <struct name="components">
-    <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>
-    </struct>
-    <struct name="linearMesh-generator">
-      <param name="Type">CartesianGenerator</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="dims">dim</param>
-      <param name="shadowDepth">shadowDepth</param>
-      <list name="size">
-        <param>elementResI</param>
-        <param>elementResJ</param>
-        <param>elementResK</param>
-      </list>
-      <list name="minCoord">
-        <param>minX</param>
-        <param>minY</param>
-        <param>minZ</param>
-      </list>
-      <list name="maxCoord">
-        <param>maxX</param>
-        <param>maxY</param>
-        <param>maxZ</param>
-      </list>
-    </struct>
-    <struct name="velocity">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="Rank">Vector</param>
-      <param name="DataType">Double</param>
-      <param name="VectorComponentCount">dim</param>
-      <list name="names">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="velocityBCs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityICs">
-      <param name="Type">CompositeVC</param>
-      <param name="Data">mesh-linear</param>
-    </struct>
-    <struct name="velocityDofLayout">
-      <param name="Type">DofLayout</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="BaseVariableCount">dim</param>
-      <list name="BaseVariables">
-        <param>vx</param>
-        <param>vy</param>
-        <param>vz</param>
-      </list>
-    </struct>
-    <struct name="VelocityField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-linear</param>
-      <param name="DofLayout">velocityDofLayout</param>
-      <param name="BC">velocityBCs</param>
-      <param name="IC">velocityICs</param>
-      <param name="LinkedDofInfo">velocityLinkedDofs</param>
-    </struct>
-    <struct name="VelocityMagnitudeField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Magnitude</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">Gradient</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="VelocityGradientsInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorInvariant</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorSymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="VorticityField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">TensorAntisymmetricPart</param>
-      <param name="FeVariable">VelocityGradientsField</param>
-    </struct>
-    <struct name="StrainRateInvariantField">
-      <param name="Type">OperatorFeVariable</param>
-      <param name="Operator">SymmetricTensor_Invariant</param>
-      <param name="FeVariable">StrainRateField</param>
-    </struct>
-    <struct name="pressure">
-      <param name="Type">MeshVariable</param>
-      <param name="mesh">mesh-constant</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>
-      <list name="BaseVariables">
-        <param>pressure</param>
-      </list>
-    </struct>
-    <struct name="PressureField">
-      <param name="Type">FeVariable</param>
-      <param name="FEMesh">mesh-constant</param>
-      <param name="DofLayout">pressureDofLayout</param>
-      <param name="LinkedDofInfo">pressureLinkedDofs</param>
-    </struct>
-    <struct name="cellLayout">
-      <param name="Type">SingleCellLayout</param>
-    </struct>
-    <struct name="particleLayout">
-      <param name="Type">GaussParticleLayout</param>
-    </struct>
-    <struct name="gaussSwarm">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">particleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">gaussMapper</param>
-    </struct>
-    <struct name="gaussMapper">
-      <param name="Type">GaussMapper</param>
-      <param name="IntegrationPointsSwarm">gaussSwarm</param>
-      <param name="MaterialPointsSwarm">gaussMaterialSwarm</param>
-    </struct>
-    <struct name="backgroundLayout">
-      <param name="Type">BackgroundParticleLayout</param>
-    </struct>
-    <struct name="gaussMaterialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">cellLayout</param>
-      <param name="ParticleLayout">backgroundLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-    </struct>
-    <struct name="timeIntegrator">
-      <param name="Type">TimeIntegrator</param>
-      <param name="order">1</param>
-      <param name="simultaneous">t</param>
-      <param name="Context">context</param>
-    </struct>
-    <struct name="elementCellLayout">
-      <param name="Type">ElementCellLayout</param>
-      <param name="Mesh">mesh-linear</param>
-    </struct>
-    <struct name="discreteVoronoi">
-      <param name="Type">CellularAutomataVoronoi</param>
-      <param name="resolution">10</param>
-      <param name="diagonalNeighbours">true</param>
-    </struct>
-    <struct name="weights">
-      <param name="Type">DiscreteVoronoiWeights</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="localLayout">
-      <param name="Type">MappedParticleLayout</param>
-    </struct>
-    <struct name="picIntegrationPoints">
-      <param name="Type">IntegrationPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">localLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="WeightsCalculator">weights</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="IntegrationPointMapper">mapper</param>
-    </struct>
-    <struct name="mapper">
-      <param name="Type">CoincidentMapper</param>
-      <param name="IntegrationPointsSwarm">picIntegrationPoints</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="materialSwarmParticleLayout">
-      <param name="Type">RandomParticleLayout</param>
-      <param name="seed">13</param>
-      <param name="cellParticleCount">particlesPerCell</param>
-    </struct>
-    <struct name="materialSwarm">
-      <param name="Type">MaterialPointsSwarm</param>
-      <param name="CellLayout">elementCellLayout</param>
-      <param name="ParticleLayout">materialSwarmParticleLayout</param>
-      <param name="FeMesh">mesh-linear</param>
-      <param name="SplittingRoutine">splittingRoutine</param>
-      <param name="RemovalRoutine">removalRoutine</param>
-      <param name="EscapedRoutine">escapedRoutine</param>
-    </struct>
-    <struct name="materialSwarmAdvector">
-      <param name="Type">SwarmAdvector</param>
-      <param name="Swarm">materialSwarm</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="VelocityField">VelocityField</param>
-      <param name="PeriodicBCsManager">periodicBCsManager</param>
-      <param name="allowFallbackToFirstOrder">True</param>
-    </struct>
-    <struct name="removalRoutine">
-      <param name="Type">DiscreteVoronoiRemove</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="splittingRoutine">
-      <param name="Type">ReseedSplitting</param>
-      <param name="idealParticleCount">particlesPerCell</param>
-      <param name="originalParticleShape">leftBoundaryShape</param>
-      <param name="DiscreteVoronoi">discreteVoronoi</param>
-    </struct>
-    <struct name="solutionVelocity">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">VelocityField</param>
-    </struct>
-    <struct name="solutionPressure">
-      <param name="Type">SolutionVector</param>
-      <param name="FeVariable">PressureField</param>
-    </struct>
-    <struct name="mom_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">VelocityField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="cont_force">
-      <param name="Type">ForceVector</param>
-      <param name="FeVariable">PressureField</param>
-      <param name="ExtraInfo">context</param>
-    </struct>
-    <struct name="k_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">VelocityField</param>
-      <param name="RHS">mom_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="constitutiveMatrix">
-      <param name="Type">ConstitutiveMatrixCartesian</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">k_matrix</param>
-    </struct>
-    <struct name="g_matrix">
-      <param name="Type">StiffnessMatrix</param>
-      <param name="RowVariable">VelocityField</param>
-      <param name="ColumnVariable">PressureField</param>
-      <param name="RHS">cont_force</param>
-      <param name="allowZeroElementContributions">False</param>
-    </struct>
-    <struct name="gradientStiffnessMatrixTerm">
-      <param name="Type">GradientStiffnessMatrixTerm</param>
-      <param name="Swarm">gaussSwarm</param>
-      <param name="StiffnessMatrix">g_matrix</param>
-    </struct>
-    <struct name="preconditioner">
-      <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="preconditionerTerm">
-      <param name="Type">UzawaPreconditionerTerm</param>
-      <param name="Swarm">picIntegrationPoints</param>
-      <param name="StiffnessMatrix">preconditioner</param>
-    </struct>
-    <struct name="uzawa">
-      <param name="Type">Stokes_SLE_UzawaSolver</param>
-      <param name="Preconditioner">preconditioner</param>
-      <param name="tolerance">1.0e-5</param>
-      <param name="maxIterations">5000</param>
-    </struct>
-    <struct name="stokesEqn">
-      <param name="Type">Stokes_SLE</param>
-      <param name="SLE_Solver">uzawa</param>
-      <param name="Context">context</param>
-      <param name="StressTensorMatrix">k_matrix</param>
-      <param name="GradientMatrix">g_matrix</param>
-      <param name="DivergenceMatrix"></param>
-      <param name="CompressibilityMatrix">c_matrix</param>
-      <param name="VelocityVector">solutionVelocity</param>
-      <param name="PressureVector">solutionPressure</param>
-      <param name="ForceVector">mom_force</param>
-      <param name="ContinuityForceVector">cont_force</param>
-      <param name="killNonConvergent">false</param>
-      <param name="nonLinearMaxIterations">nonLinearMaxIterations</param>
-      <param name="nonLinearTolerance">nonLinearTolerance</param>
-      <param name="makeConvergenceFile">false</param>
-    </struct>
-    <struct name="backgroundShape">
-      <param name="Type">Everywhere</param>
-    </struct>
-    <struct name="leftBoundaryShape">
-      <param name="Type">Box</param>
-      <param name="startX">minX</param>
-      <param name="endX">0.1</param>
-      <param name="startY">minY</param>
-      <param name="endY">maxY</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-    <struct name="rightBoundaryShape">
-      <param name="Type">Box</param>
-      <param name="startX">1.9</param>
-      <param name="endX">maxX</param>
-      <param name="startY">minY</param>
-      <param name="endY">maxY</param>
-      <param name="startZ">minZ</param>
-      <param name="endZ">maxZ</param>
-    </struct>
-
-    <struct name="leftSlabOuterSphere">
-      <param name="Type">Sphere</param>
-      <param name="CentreX">0.0</param>
-      <param name="CentreY">-7.15</param>
-      <param name="radius">7.25</param>
-    </struct>
-    <struct name="leftSlabInnerSphere">
-      <param name="Type">Sphere</param>
-      <param name="CentreX">0.0</param>
-      <param name="CentreY">-6.15</param>
-      <param name="radius">6.2</param>
-    </struct>
-    <struct name="leftSlabShape">
-      <param name="Type">Intersection</param>
-      <list name="shapes">
-        <param>leftSlabOuterSphere</param>
-        <param>!leftSlabInnerSphere</param>
-      </list>
-    </struct>
-
-
-
-    <struct name="rightSlabOuterSphere">
-      <param name="Type">Sphere</param>
-      <param name="CentreX">maxX</param>
-      <param name="CentreY">-7.15</param>
-      <param name="radius">7.25</param>
-    </struct>
-    <struct name="rightSlabInnerSphere">
-      <param name="Type">Sphere</param>
-      <param name="CentreX">maxX</param>
-      <param name="CentreY">-6.15</param>
-      <param name="radius">6.2</param>
-    </struct>
-    <struct name="rightSlabShape">
-      <param name="Type">Intersection</param>
-      <list name="shapes">
-              <param>rightSlabOuterSphere</param>
-        <param>!rightSlabInnerSphere</param>
-        <param>!leftSlabShape</param>
-      </list>
-    </struct>
-    
-    <struct name="beadsShape">
-      <param name="Type">Intersection</param>
-      <list name="shapes">
-        <param>backgroundShape</param>
-        <param>!leftSlabShape</param>
-        <param>!rightSlabShape</param>
-      </list>
-    </struct>
-    <struct name="crustShape">
-      <param name="Type">Union</param>
-      <list name="shapes">
-        <param>leftSlabShape</param>
-        <param>rightSlabShape</param>
-      </list>
-    </struct>
-    <struct name="crustViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">10.0</param>
-    </struct>
-    <struct name="strainWeakening">
-      <param name="Type">StrainWeakening</param>
-      <param name="TimeIntegrator">timeIntegrator</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="softeningStrain">0.1</param>
-      <param name="initialDamageFraction">0.0</param>
-      <param name="initialDamageWavenumber">0.5</param>
-      <param name="initialDamageFactor">0.5</param>
-      <param name="healingRate">0.0</param>
-    </struct>
-    <struct name="yielding">
-      <param name="Type">MohrCoulomb</param>
-      <param name="PressureField">PressureField</param>
-      <param name="VelocityGradientsField">VelocityGradientsField</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-      <param name="Context">context</param>
-      <param name="StrainWeakening">strainWeakening</param>
-      <param name="cohesion">100.0</param>
-      <param name="cohesionAfterSoftening">0.0001</param>
-      <param name="frictionCoefficient">0.0</param>
-      <param name="frictionCoefficientAfterSoftening">0.0</param>
-      <param name="minimumYieldStress">0.00001</param>
-    </struct>
-    <struct name="beadsViscosity">
-      <param name="Type">MaterialViscosity</param>
-      <param name="eta0">0.1</param>
-    </struct>
-    <struct name="storeViscosity">
-      <param name="Type">StoreVisc</param>
-      <param name="MaterialPointsSwarm">materialSwarm</param>
-    </struct>
-    <struct name="crust">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">crustShape</param>
-      <list name="Rheology">
-        <param>crustViscosity</param>
-        <param>yielding</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="beads">
-      <param name="Type">RheologyMaterial</param>
-      <param name="Shape">beadsShape</param>
-      <list name="Rheology">
-        <param>beadsViscosity</param>
-        <param>storeViscosity</param>
-      </list>
-    </struct>
-    <struct name="escapedRoutine">
-      <param name="Type">EscapedRoutine</param>
-      <param name="idealParticleCount">0</param>
-    </struct>
-    <struct name="velocityRemesher">
-      <param name="Type">StripRemesher</param>
-      <param name="mesh">mesh-linear</param>
-      <param name="meshType">regular</param>
-      <list name="dims">
-        <param>true</param>
-        <param>true</param>
-        <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="dims">3</param>
-      <param name="cellMesh">mesh-linear</param>
-    </struct>
-    <struct name="conditionFunctions">
-      <param name="Type">StgFEM_StandardConditionFunctions</param>
-    </struct>
-  </struct>
-  <list name="plugins">
-    <param>Underworld_EulerDeform</param>
-    <param>Underworld_DumpSwarm</param>
-    <param>Gale_SurfaceProcess</param>
-  </list>
-  <param name="maxTimeSteps">500</param>
-  <param name="outputEvery">1</param>
-  <param name="dumpEvery">1</param>
-  <param name="outputPath">./output.subduction</param>
-  <param name="dim">2</param>
-  <param name="shadowDepth">1</param>
-  <param name="minX">0.0f</param>
-  <param name="minY">0.0f</param>
-  <param name="minZ">0.0f</param>
-  <param name="maxX">2.0f</param>
-  <param name="maxY">0.2f</param>
-  <param name="maxZ">1.0f</param>
-  <param name="elementResI">32</param>
-  <param name="elementResJ">16</param>
-  <param name="elementResK">4</param>
-  <param name="allowUnbalancing">True</param>
-  <param name="buildNodeNeighbourTbl">True</param>
-  <param name="buildElementNeighbourTbl">True</param>
-  <param name="particlesPerCell">30</param>
-  <param name="seed">13</param>
-  <param name="gaussParticlesX">2</param>
-  <param name="gaussParticlesY">2</param>
-  <param name="gaussParticlesZ">2</param>
-  <param name="particlesPerCell">30</param>
-  <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>
-        <param name="wrapTop">True</param>
-        <param name="staticLeft">True</param>
-        <param name="staticBottom">True</param>
-        <list name="fields">
-          <struct>
-            <param name="field">VelocityField</param>
-            <param name="variable">velocity</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <struct name="velocityBCs">
-    <param name="type">CompositeVC</param>
-    <list name="vcList">
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">front</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">back</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vz</param>
-            <param name="type">double</param>
-            <param name="value">0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">left</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">func</param>
-            <param name="value">StepFunctionProduct3</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">right</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">double</param>
-            <param name="value">0.0</param>
-          </struct>
-        </list>
-      </struct>
-      <struct>
-        <param name="type">WallVC</param>
-        <param name="wall">bottom</param>
-        <list name="variables">
-          <struct>
-            <param name="name">vy</param>
-            <param name="type">func</param>
-            <param name="value">StepFunctionProduct2</param>
-          </struct>
-          <struct>
-            <param name="name">vx</param>
-            <param name="type">func</param>
-            <param name="value">StepFunctionProduct1</param>
-          </struct>
-        </list>
-      </struct>
-    </list>
-  </struct>
-  <param name="StepFunctionProduct1Start">0.8</param>
-  <param name="StepFunctionProduct1End">1.2</param>
-  <param name="StepFunctionProduct1Value">0.5</param>
-  <param name="StepFunctionProduct2Start">0.8</param>
-  <param name="StepFunctionProduct2End">1.2</param>
-  <param name="StepFunctionProduct2Value">-0.1</param>
-
-  <param name="StepFunctionProduct3Start">0.05</param>
-  <param name="StepFunctionProduct3End">0.1</param>
-  <param name="StepFunctionProduct3Value">0.5</param>
-
-  <struct name="SurfaceProcess">
-    <param name="mesh">mesh-linear</param>
-  </struct>
-
-  <param name="checkpointEvery">1</param>
-</StGermainData>



More information about the cig-commits mailing list