[cig-commits] r13429 - in long/3D/Gale/trunk: . input/examples
walter at geodynamics.org
walter at geodynamics.org
Mon Dec 1 10:56:59 PST 2008
Author: walter
Date: 2008-12-01 10:56:59 -0800 (Mon, 01 Dec 2008)
New Revision: 13429
Modified:
long/3D/Gale/trunk/
long/3D/Gale/trunk/input/examples/tibet.xml
Log:
r2382 at dante: boo | 2008-12-01 10:55:12 -0800
Make tibet.xml work
Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
- 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2381
+ 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2382
Modified: long/3D/Gale/trunk/input/examples/tibet.xml
===================================================================
--- long/3D/Gale/trunk/input/examples/tibet.xml 2008-12-01 18:56:57 UTC (rev 13428)
+++ long/3D/Gale/trunk/input/examples/tibet.xml 2008-12-01 18:56:59 UTC (rev 13429)
@@ -115,6 +115,8 @@
<param name="FEMesh">mesh-linear</param>
<param name="DofLayout">pressureDofLayout</param>
<param name="LinkedDofInfo">pressureLinkedDofs</param>
+ <param name="BC">pressureBCs</param>
+ <param name="IC">pressureICs</param>
</struct>
<struct name="StressField">
<param name="Type">StressField</param>
@@ -282,7 +284,7 @@
<struct name="uzawa">
<param name="Type">Stokes_SLE_UzawaSolver</param>
<param name="Preconditioner">preconditioner</param>
- <param name="tolerance">1.0e-5</param>
+ <param name="tolerance">linearTolerance</param>
<param name="maxIterations">5000</param>
</struct>
<struct name="stokesEqn">
@@ -302,13 +304,29 @@
<param name="nonLinearTolerance">1e-2</param>
<param name="makeConvergenceFile">false</param>
</struct>
-<!-- <struct name="buoyancyForceTerm"> -->
-<!-- <param name="Type">BuoyancyForceTerm</param> -->
-<!-- <param name="ForceVector">mom_force</param> -->
-<!-- <param name="TemperatureField">TemperatureField</param> -->
-<!-- <param name="Swarm">picIntegrationPoints</param> -->
-<!-- <param name="gravity">9.81</param> -->
-<!-- </struct> -->
+
+ <struct name="c_matrix">
+ <param name="Type">StiffnessMatrix</param>
+ <param name="RowVariable">PressureField</param>
+ <param name="ColumnVariable">PressureField</param>
+ <param name="RHS">cont_force</param>
+ <param name="allowZeroElementContributions">True</param>
+ </struct>
+ <struct name="mixedStabiliser">
+ <param name="Type">MixedStabiliserTerm</param>
+ <param name="Swarm">gaussSwarm</param>
+ <param name="picSwarm">picIntegrationPoints</param>
+ <param name="storeVisc">storeViscosity</param>
+ <param name="StiffnessMatrix">c_matrix</param>
+ </struct>
+
+ <struct name="buoyancyForceTerm">
+ <param name="Type">BuoyancyForceTerm</param>
+ <param name="ForceVector">mom_force</param>
+ <param name="TemperatureField">TemperatureField</param>
+ <param name="Swarm">picIntegrationPoints</param>
+ <param name="gravity">9.81</param>
+ </struct>
<struct name="background">
<param name="Type">Everywhere</param>
</struct>
@@ -326,21 +344,6 @@
<param>true</param>
</list>
</struct>
-
- <element type="struct" name="c_matrix">
- <element type="param" name="Type">StiffnessMatrix</element>
- <element type="param" name="RowVariable">PressureField</element>
- <element type="param" name="ColumnVariable">PressureField</element>
- <element type="param" name="RHS">cont_force</element>
- <element type="param" name="allowZeroElementContributions">True</element>
- </element>
- <element type="struct" name="mixedStabiliser">
- <element type="param" name="Type">MixedStabiliserTerm</element>
- <element type="param" name="Swarm">gaussSwarm</element>
- <element type="param" name="materialSwarm">materialSwarm</element>
- <element type="param" name="StiffnessMatrix">c_matrix</element>
- </element>
-
<struct name="bottomBoundaryShape">
<param name="Type">Box</param>
<param name="startX">minX</param>
@@ -361,21 +364,21 @@
<!-- <param name="y_value">3080340000</param> -->
<!-- </struct> -->
-<!-- <struct name="surfaceAdaptor"> -->
-<!-- <param name="Type">SurfaceAdaptor</param> -->
-<!-- <param name="mesh">mesh-linear</param> -->
-<!-- <param name="sourceGenerator">linearMesh-generator</param> -->
-<!-- <param name="surfaceType">plateau</param> -->
-<!-- <param name="x1">300000</param> -->
-<!-- <param name="x2">350000</param> -->
-<!-- <param name="x3">650000</param> -->
-<!-- <param name="x4">700000</param> -->
-<!-- <param name="z1">200000</param> -->
-<!-- <param name="z2">250000</param> -->
-<!-- <param name="z3">750000</param> -->
-<!-- <param name="z4">800000</param> -->
-<!-- <param name="height">3030.3030303</param> -->
-<!-- </struct> -->
+ <struct name="surfaceAdaptor">
+ <param name="Type">SurfaceAdaptor</param>
+ <param name="mesh">mesh-linear</param>
+ <param name="sourceGenerator">linearMesh-generator</param>
+ <param name="surfaceType">plateau</param>
+ <param name="x1">300000</param>
+ <param name="x2">350000</param>
+ <param name="x3">650000</param>
+ <param name="x4">700000</param>
+ <param name="z1">200000</param>
+ <param name="z2">250000</param>
+ <param name="z3">750000</param>
+ <param name="z4">800000</param>
+ <param name="height">3030.3030303</param>
+ </struct>
<struct name="boxShape">
<param name="Type">Box</param>
@@ -400,91 +403,91 @@
</struct>
<!-- Crust -->
-<!-- <struct name="crustShape"> -->
-<!-- <param name="Type">PolygonShape</param> -->
-<!-- <list name="verticies"> -->
-<!-- <asciidata> -->
-<!-- <columnDefinition name = "x" type="double"/> -->
-<!-- <columnDefinition name = "y" type="double"/> -->
-<!-- 0 68000 -->
-<!-- 300000 68000 -->
-<!-- 350000 51030.3030303 -->
-<!-- 650000 51030.3030303 -->
-<!-- 700000 68000 -->
-<!-- 1000000 68000 -->
-<!-- 1000000 100000 -->
-<!-- 0 100000 -->
-<!-- </asciidata> -->
-<!-- </list> -->
-<!-- </struct> -->
+ <struct name="crustShape">
+ <param name="Type">PolygonShape</param>
+ <list name="verticies">
+ <asciidata>
+ <columnDefinition name = "x" type="double"/>
+ <columnDefinition name = "y" type="double"/>
+ 0 68000
+ 300000 68000
+ 350000 51030.3030303
+ 650000 51030.3030303
+ 700000 68000
+ 1000000 68000
+ 1000000 100000
+ 0 100000
+ </asciidata>
+ </list>
+ </struct>
-<!-- <struct name="crustViscosity"> -->
-<!-- <param name="Type">FrankKamenetskii</param> -->
-<!-- <param name="TemperatureField">TemperatureField</param> -->
+ <struct name="crustViscosity">
+ <param name="Type">FrankKamenetskii</param>
+ <param name="TemperatureField">TemperatureField</param>
<!-- at T=273, viscosity=1e25, at T=1273, viscosity=1e20 -->
+ <param name="eta0">2.31739465</param>
<!-- <param name="eta0">2.31739465e26</param> -->
-<!-- <param name="theta">0.011512925465</param> -->
-<!-- </struct> -->
+ <param name="theta">0.011512925465</param>
+ </struct>
-<!-- <struct name="strainWeakening"> -->
-<!-- <param name="Type">StrainWeakening</param> -->
-<!-- <param name="TimeIntegrator">timeIntegrator</param> -->
-<!-- <param name="MaterialPointsSwarm">materialSwarm</param> -->
-<!-- <param name="finalSofteningStrain">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">44e6</param> -->
-<!-- <param name="cohesionAfterSoftening">4e6</param> -->
-<!-- <param name="frictionCoefficient">0.57735026919</param> -->
-<!-- <param name="frictionCoefficientAfterSoftening">0.57735026919</param> -->
-<!-- </struct> -->
+ <struct name="strainWeakening">
+ <param name="Type">StrainWeakening</param>
+ <param name="TimeIntegrator">timeIntegrator</param>
+ <param name="MaterialPointsSwarm">materialSwarm</param>
+ <param name="finalSofteningStrain">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">44e6</param>
+ <param name="cohesionAfterSoftening">4e6</param>
+ <param name="frictionCoefficient">0.57735026919</param>
+ <param name="frictionCoefficientAfterSoftening">0.57735026919</param>
+ </struct>
-<!-- <struct name="crustViscous"> -->
-<!-- <param name="Type">RheologyMaterial</param> -->
-<!-- <param name="Shape">crustShape</param> -->
-<!-- <param name="density">2800</param> -->
-<!-- <param name="alpha">3.0e-5</param> -->
-<!-- <param name="diffusivity">1.0e-6</param> -->
-<!-- <list name="heatingElements"> -->
-<!-- <struct> -->
-<!-- <param name="Q">1.0e-12</param> -->
-<!-- <param name="lambda">0.0</param> -->
-<!-- </struct> -->
-<!-- </list> -->
-<!-- <list name="Rheology"> -->
-<!-- <param>crustViscosity</param> -->
-<!-- <param>yielding</param> -->
-<!-- <param>storeViscosity</param> -->
-<!-- <param>storeStress</param> -->
-<!-- </list> -->
-<!-- </struct> -->
+ <struct name="crustViscous">
+ <param name="Type">RheologyMaterial</param>
+ <param name="Shape">crustShape</param>
+ <param name="density">2800</param>
+ <param name="alpha">3.0e-5</param>
+ <param name="diffusivity">1.0e-6</param>
+ <list name="heatingElements">
+ <struct>
+ <param name="Q">1.0e-12</param>
+ <param name="lambda">0.0</param>
+ </struct>
+ </list>
+ <list name="Rheology">
+ <param>crustViscosity</param>
+ <param>yielding</param>
+ <param>storeViscosity</param>
+ <param>storeStress</param>
+ </list>
+ </struct>
<!-- Mantle -->
<struct name="mantleShape">
<param name="Type">Intersection</param>
<list name="shapes">
<param>boxShape</param>
-<!-- <param>!crustShape</param> -->
+ <param>!crustShape</param>
</list>
</struct>
<struct name="mantleViscosity">
- <param name="Type">MaterialViscosity</param>
- <param name="eta0">1.0</param>
-<!-- <param name="Type">FrankKamenetskii</param> -->
-<!-- <param name="TemperatureField">TemperatureField</param> -->
+ <param name="Type">FrankKamenetskii</param>
+ <param name="TemperatureField">TemperatureField</param>
+ <param name="eta0">2.31739465</param>
<!-- <param name="eta0">2.31739465e26</param> -->
-<!-- <param name="theta">0.011512925465</param> -->
+ <param name="theta">0.011512925465</param>
</struct>
<struct name="mantleViscous">
@@ -588,10 +591,9 @@
</struct>
<list name="plugins">
<param>Underworld_EulerDeform</param>
-<!-- <param>Underworld_DumpSwarm</param> -->
<param>Underworld_VTKOutput</param>
</list>
- <param name="maxTimeSteps">10</param>
+ <param name="maxTimeSteps">1000</param>
<param name="outputEvery">1</param>
<param name="dumpEvery">1</param>
<param name="outputPath">./tibet</param>
@@ -680,7 +682,9 @@
<param name="name">vx</param>
<param name="type">double</param>
<!-- 1 cm/yr -->
- <param name="value">3.17097919838e-10</param>
+ <param name="value">3.17097919838e17</param>
+<!-- <param name="value">3.17097919838e16</param> -->
+<!-- <param name="value">3.17097919838e-10</param> -->
</struct>
</list>
</struct>
More information about the CIG-COMMITS
mailing list