<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Walter - can you explain whether there was an actual bug in the code that could<div>affect other people's modeling?</div><div>- Magali</div><div><br><div><div>On Oct 14, 2010, at 11:29 AM, Walter Landry wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Walter Landry <<a href="mailto:walter@geodynamics.org">walter@geodynamics.org</a>> wrote:<br><blockquote type="cite">Karen Paczkowski <<a href="mailto:karen.paczkowski@yale.edu">karen.paczkowski@yale.edu</a>> wrote:<br></blockquote><blockquote type="cite"><blockquote type="cite">Hi,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Thank you for all the help. I am unfortunately still getting a very<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">unexpected flow field. The only velocities I have set are:<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Vx = 0 and Vy = -Vslab on the left side<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">Vy=0 across the top<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite"><br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">My input file is producing a flow field that flows in the opposite<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">direction as Vslab (up and to the right) and has a magnitude 2-3<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">orders of magnitude larger than Vslab. This was not occurring prior<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">to adding the thermal field, so I am assuming it has to do with that,<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">but I'm really not sure what is causing this exactly. I have attached<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">a picture of the velocity field and a copy of the input file. It is<br></blockquote></blockquote><blockquote type="cite"><blockquote type="cite">written for Gale 1.5.0. I would really appreciate any advice.<br></blockquote></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The problem is that there is no maximum temperature in the<br></blockquote><blockquote type="cite">HydrostaticTerm. I am on vacation until September, so I will not be<br></blockquote><blockquote type="cite">able to get to it until then.<br></blockquote><br>Sorry for the long delay, but I finally implemented this. You can get<br>the code from the repository. I am attaching a new input file and a<br>picture of the velocity and temperature. Please let me know if you<br>have any more questions.<br><br>Cheers,<br>Walter Landry<br><a href="mailto:walter@geodynamics.org">walter@geodynamics.org</a><br><?xml version="1.0"?><br><StGermainData xmlns="http://www.vpac.org/StGermain/XML_IO_Handler/Jun2003"><br> <import><br> <toolbox>Underworld</toolbox><br> </import><br> <plugins><br> <struct><br> <param name="Type">Underworld_EulerDeform</param><br> <param name="Context">context</param><br> </struct><br> <struct><br> <param name="Type">Underworld_VTKOutput</param><br> <param name="Context">context</param><br> </struct><br> <struct><br> <param name="Type">StgFEM_StandardConditionFunctions</param><br> <param name="Context">context</param><br> </struct><br> </plugins><br> <struct name="components"><br> <struct name="context"><br> <param name="Type">UnderworldContext</param><br> </struct><br> <struct name="mesh-linear"><br> <param name="Type">FeMesh</param><br> <param name="elementType">linear</param><br> </struct><br> <struct name="linearMesh-generator"><br> <param name="Type">CartesianGenerator</param><br> <param name="mesh">mesh-linear</param><br> <param name="dim">dim</param><br> <param name="shadowDepth">shadowDepth</param><br> <param name="regular">False</param><br> <list name="size"><br> <param>elementResI</param><br> <param>elementResJ</param><br> <param>elementResK</param><br> </list><br> <list name="minCoord"><br> <param>minX</param><br> <param>minY</param><br> <param>minZ</param><br> </list><br> <list name="maxCoord"><br> <param>maxX</param><br> <param>maxY</param><br> <param>maxZ</param><br> </list><br> </struct><br> <struct name="velocity"><br> <param name="Type">MeshVariable</param><br> <param name="mesh">mesh-linear</param><br> <param name="Rank">Vector</param><br> <param name="DataType">Double</param><br> <param name="VectorComponentCount">dim</param><br> <list name="names"><br> <param>vx</param><br> <param>vy</param><br> <param>vz</param><br> </list><br> </struct><br> <struct name="velocityBCs"><br> <param name="Type">CompositeVC</param><br> <param name="Data">mesh-linear</param><br> </struct><br> <struct name="velocityICs"><br> <param name="Type">CompositeVC</param><br> <param name="Data">mesh-linear</param><br> </struct><br> <struct name="velocityDofLayout"><br> <param name="Type">DofLayout</param><br> <param name="mesh">mesh-linear</param><br> <param name="BaseVariableCount">dim</param><br> <list name="BaseVariables"><br> <param>vx</param><br> <param>vy</param><br> <param>vz</param><br> </list><br> </struct><br> <struct name="VelocityField"><br> <param name="Type">FeVariable</param><br> <param name="FEMesh">mesh-linear</param><br> <param name="DofLayout">velocityDofLayout</param><br> <param name="BC">velocityBCs</param><br> <param name="IC">velocityICs</param><br> <param name="LinkedDofInfo">velocityLinkedDofs</param><br> </struct><br> <struct name="VelocityMagnitudeField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">Magnitude</param><br> <param name="FeVariable">VelocityField</param><br> </struct><br> <struct name="VelocityGradientsField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">Gradient</param><br> <param name="FeVariable">VelocityField</param><br> </struct><br> <struct name="VelocityGradientsInvariantField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">TensorInvariant</param><br> <param name="FeVariable">VelocityGradientsField</param><br> </struct><br> <struct name="StrainRateField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">TensorSymmetricPart</param><br> <param name="FeVariable">VelocityGradientsField</param><br> </struct><br> <struct name="VorticityField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">TensorAntisymmetricPart</param><br> <param name="FeVariable">VelocityGradientsField</param><br> </struct><br> <struct name="StrainRateInvariantField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">SymmetricTensor_Invariant</param><br> <param name="FeVariable">StrainRateField</param><br> </struct><br> <struct name="pressure"><br> <param name="Type">MeshVariable</param><br> <param name="mesh">mesh-linear</param><br> <param name="Rank">Scalar</param><br> <param name="DataType">Double</param><br> </struct><br> <struct name="pressureDofLayout"><br> <param name="Type">DofLayout</param><br> <param name="mesh">mesh-linear</param><br> <list name="BaseVariables"><br> <param>pressure</param><br> </list><br> </struct><br> <struct name="PressureField"><br> <param name="Type">FeVariable</param><br> <param name="FEMesh">mesh-linear</param><br> <param name="DofLayout">pressureDofLayout</param><br> <param name="LinkedDofInfo">pressureLinkedDofs</param><br> </struct><br> <struct name="StressField"><br> <param name="Type">StressField</param><br> <param name="StrainRateField">StrainRateField</param><br> <param name="Context">context</param><br> <param name="ConstitutiveMatrix">constitutiveMatrix</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="Mesh">mesh-linear</param><br> <param name="IC">stressICs</param><br> </struct> <br> <struct name="ViscosityField"><br> <param name="Type">ViscosityField</param><br> <param name="Context">context</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="Mesh">mesh-linear</param><br> <param name="ConstitutiveMatrix">constitutiveMatrix</param><br> </struct><br> <struct name="cellLayout"><br> <param name="Type">SingleCellLayout</param><br> </struct><br> <struct name="particleLayout"><br> <param name="Type">GaussParticleLayout</param><br> </struct><br> <struct name="gaussSwarm"><br> <param name="Type">IntegrationPointsSwarm</param><br> <param name="CellLayout">cellLayout</param><br> <param name="ParticleLayout">particleLayout</param><br> <param name="FeMesh">mesh-linear</param><br> <param name="TimeIntegrator">timeIntegrator</param><br> <param name="IntegrationPointMapper">gaussMapper</param><br> </struct><br> <struct name="gaussMapper"><br> <param name="Type">GaussMapper</param><br> <param name="IntegrationPointsSwarm">gaussSwarm</param><br> <param name="MaterialPointsSwarm">gaussMaterialSwarm</param><br> </struct><br> <struct name="backgroundLayout"><br> <param name="Type">BackgroundParticleLayout</param><br> </struct><br> <struct name="gaussMSwarmMovementHandler"><br> <param name="Type">ParticleMovementHandler</param><br> </struct><br> <struct name="gaussMaterialSwarm"><br> <param name="Type">MaterialPointsSwarm</param><br> <param name="CellLayout">cellLayout</param><br> <param name="ParticleLayout">backgroundLayout</param><br> <param name="FeMesh">mesh-linear</param><br> <list name="ParticleCommHandlers"><br> <param>gaussMSwarmMovementHandler</param><br> </list><br> </struct><br> <struct name="timeIntegrator"><br> <param name="Type">TimeIntegrator</param><br> <param name="order">1</param><br> <param name="simultaneous">t</param><br> <param name="Context">context</param><br> </struct><br> <struct name="elementCellLayout"><br> <param name="Type">ElementCellLayout</param><br> <param name="Mesh">mesh-linear</param><br> </struct><br> <struct name="weights"><br> <param name="Type">PCDVC</param><br> <param name="resolutionX">10</param><br> <param name="resolutionY">10</param><br> <param name="resolutionZ">10</param><br> <param name="lowerT">0.6</param><br> <param name="upperT">25</param><br> <param name="maxDeletions">3</param><br> <param name="maxSplits">3</param><br> <param name="MaterialPointsSwarm">materialSwarm</param><br> </struct><br> <struct name="localLayout"><br> <param name="Type">MappedParticleLayout</param><br> </struct><br> <struct name="picIntegrationPoints"><br> <param name="Type">IntegrationPointsSwarm</param><br> <param name="CellLayout">elementCellLayout</param><br> <param name="ParticleLayout">localLayout</param><br> <param name="FeMesh">mesh-linear</param><br> <param name="WeightsCalculator">weights</param><br> <param name="TimeIntegrator">timeIntegrator</param><br> <param name="IntegrationPointMapper">mapper</param><br> </struct><br> <struct name="mapper"><br> <param name="Type">CoincidentMapper</param><br> <param name="IntegrationPointsSwarm">picIntegrationPoints</param><br> <param name="MaterialPointsSwarm">materialSwarm</param><br> </struct><br> <struct name="materialSwarmParticleLayout"><br> <param name="Type">MeshParticleLayout</param><br> <param name="mesh">mesh-linear</param><br> <param name="cellParticleCount">particlesPerCell</param><br> </struct><br> <struct name="pMovementHandler"><br> <param name="Type">ParticleMovementHandler</param><br> </struct><br> <struct name="pShadowSync"><br> <param name="Type">ParticleShadowSync</param><br> </struct><br> <struct name="materialSwarm"><br> <param name="Type">MaterialPointsSwarm</param><br> <param name="CellLayout">elementCellLayout</param><br> <param name="ParticleLayout">materialSwarmParticleLayout</param><br> <param name="FeMesh">mesh-linear</param><br> <list name="ParticleCommHandlers"><br> <param>pMovementHandler</param><br> <param>pShadowSync</param><br> </list><br> <param name="SplittingRoutine">splittingRoutine</param><br> <param name="RemovalRoutine">removalRoutine</param><br> <param name="EscapedRoutine">escapedRoutine</param><br> </struct><br> <struct name="materialSwarmAdvector"><br> <param name="Type">SwarmAdvector</param><br> <param name="Swarm">materialSwarm</param><br> <param name="TimeIntegrator">timeIntegrator</param><br> <param name="VelocityField">VelocityField</param><br> <param name="PeriodicBCsManager">periodicBCsManager</param><br> <param name="allowFallbackToFirstOrder">True</param><br> </struct><br> <struct name="solutionVelocity"><br> <param name="Type">SolutionVector</param><br> <param name="FeVariable">VelocityField</param><br> </struct><br> <struct name="solutionPressure"><br> <param name="Type">SolutionVector</param><br> <param name="FeVariable">PressureField</param><br> </struct><br> <struct name="mom_force"><br> <param name="Type">ForceVector</param><br> <param name="FeVariable">VelocityField</param><br> <param name="ExtraInfo">context</param><br> </struct><br> <struct name="cont_force"><br> <param name="Type">ForceVector</param><br> <param name="FeVariable">PressureField</param><br> <param name="ExtraInfo">context</param><br> </struct><br> <struct name="k_matrix"><br> <param name="Type">StiffnessMatrix</param><br> <param name="RowVariable">VelocityField</param><br> <param name="ColumnVariable">VelocityField</param><br> <param name="RHS">mom_force</param><br> <param name="allowZeroElementContributions">False</param><br> </struct><br> <struct name="constitutiveMatrix"><br> <param name="Type">ConstitutiveMatrixCartesian</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="StiffnessMatrix">k_matrix</param><br> </struct><br> <struct name="g_matrix"><br> <param name="Type">StiffnessMatrix</param><br> <param name="RowVariable">VelocityField</param><br> <param name="ColumnVariable">PressureField</param><br> <param name="RHS">mom_force</param><br> <param name="transposeRHS">cont_force</param><br> <param name="allowZeroElementContributions">False</param><br> </struct><br> <struct name="gradientStiffnessMatrixTerm"><br> <param name="Type">GradientStiffnessMatrixTerm</param><br> <param name="Swarm">gaussSwarm</param><br> <param name="StiffnessMatrix">g_matrix</param><br> </struct><br> <struct name="preconditioner"><br> <param name="Type">StiffnessMatrix</param><br> <param name="RowVariable">PressureField</param><br> <param name="ColumnVariable">PressureField</param><br> <param name="RHS">cont_force</param><br> <param name="allowZeroElementContributions">True</param><br> </struct><br> <struct name="preconditionerTerm"><br> <param name="Type">UzawaPreconditionerTerm</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="StiffnessMatrix">preconditioner</param><br> </struct><br> <struct name="uzawa"><br> <param name="Type">Stokes_SLE_UzawaSolver</param><br> <param name="velocitySolver">matrixSolver</param><br> <param name="Preconditioner">preconditioner</param><br> <param name="tolerance">linearTolerance</param><br> <param name="monitor">false</param><br> <param name="maxIterations">5000</param><br> <param name="minIterations">1</param><br> </struct><br> <struct name="stokesEqn"><br> <param name="Type">Stokes_SLE</param><br> <param name="SLE_Solver">uzawa</param><br> <param name="Context">context</param><br> <param name="StressTensorMatrix">k_matrix</param><br> <param name="GradientMatrix">g_matrix</param><br> <param name="DivergenceMatrix"></param><br> <param name="CompressibilityMatrix">c_matrix</param><br> <param name="VelocityVector">solutionVelocity</param><br> <param name="PressureVector">solutionPressure</param><br> <param name="ForceVector">mom_force</param><br> <param name="ContinuityForceVector">cont_force</param><br> <param name="killNonConvergent">false</param><br> <param name="nonLinearMaxIterations">nonLinearMaxIterations</param><br> <param name="nonLinearTolerance">nonLinearTolerance</param><br> <param name="makeConvergenceFile">false</param><br> </struct><br> <struct name="c_matrix"><br> <param name="Type">StiffnessMatrix</param><br> <param name="RowVariable">PressureField</param><br> <param name="ColumnVariable">PressureField</param><br> <param name="RHS">cont_force</param><br> <param name="allowZeroElementContributions">True</param><br> </struct><br> <struct name="mixedStabiliser"><br> <param name="Type">MixedStabiliserTerm</param><br> <param name="Swarm">gaussSwarm</param><br> <param name="picSwarm">picIntegrationPoints</param><br> <param name="storeVisc">storeViscosity</param><br> <param name="StiffnessMatrix">c_matrix</param><br> </struct><br> <struct name="buoyancyForceTerm"><br> <param name="Type">BuoyancyForceTerm</param><br> <param name="ForceVector">mom_force</param><br> <param name="TemperatureField">TemperatureField</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="gravity">gravity</param><br> <param name="HydrostaticTerm">hydrostaticTerm</param><br> </struct><br> <struct name="background"><br> <param name="Type">Everywhere</param><br> </struct><br> <struct name="escapedRoutine"><br> <param name="Type">EscapedRoutine</param><br> <param name="idealParticleCount">0</param><br> </struct><br> <struct name="velocityRemesher"><br> <param name="Type">RegularRemesherCmpt</param><br> <param name="mesh">mesh-linear</param><br> <list name="remeshDims"><br> <param>0</param><br> <param>1</param><br> <param>2</param><br> </list><br> </struct><br><br> <struct name="bottomBoundaryShape"><br> <param name="Type">Box</param><br> <param name="startX">minX</param><br> <param name="endX">maxX</param><br> <param name="startY">minY</param><br> <param name="endY">10.00</param><br> <param name="startZ">minZ</param><br> <param name="endZ">maxZ</param><br> </struct><br><br> <!-- <struct name="fixedCornerShape"> --><br> <!-- <param name="Type">Box</param> --><br> <!-- <param name="startX">1999000</param> --><br> <!-- <param name="endX">maxX</param> --><br> <!-- <param name="startY">439000</param> --><br> <!-- <param name="endY">maxY</param> --><br> <!-- <param name="startZ">minZ</param> --><br> <!-- <param name="endZ">maxZ</param> --><br> <!-- </struct> --><br><br><br> <struct name="boxShape"><br> <param name="Type">Box</param><br> <param name="startX">minX</param><br> <param name="endX">maxX</param><br> <param name="startY">minY</param><br> <param name="endY">maxY</param><br> <param name="startZ">minZ</param><br> <param name="endZ">maxZ</param><br> </struct><br> <struct name="backgroundViscosity"><br> <param name="Type">MaterialViscosity</param><br> <param name="eta0">1.0</param><br> </struct><br> <struct name="storeViscosity"><br> <param name="Type">StoreVisc</param><br> <param name="MaterialPointsSwarm">materialSwarm</param><br> </struct><br> <struct name="storeStress"><br> <param name="Type">StoreStress</param><br> <param name="MaterialPointsSwarm">materialSwarm</param><br> </struct><br> <struct name="viscous"><br> <param name="Type">RheologyMaterial</param><br> <param name="Shape">boxShape</param><br> <param name="density">3000.01</param><br> <param name="alpha">3e-5</param><br> <param name="diffusivity">1e15</param><br> <list name="heatingElements"><br> <struct><br> <param name="Q">0</param><br> <param name="lambda">0</param><br> </struct><br> </list><br> <list name="Rheology"><br> <param>backgroundViscosity</param><br> <param>storeViscosity</param><br> <param>storeStress</param><br> </list><br> </struct><br><br> <struct name="hydrostaticTerm"><br> <param name="Type">HydrostaticTerm</param><br> <param name="upperDensity">3000</param><br> <param name="upperAlpha">3e-5</param><br> <param name="height">maxY</param><br> <param name="T_0">TemperatureProfileTop</param><br> <param name="T_max">TemperatureProfileMax</param><br> <param name="T_max_depth">100000</param><br> <param name="linearCoefficient">TemperatureProfileLinearCoefficient</param><br> <param name="TemperatureProfileMax">TemperatureProfileMax</param><br> <param name="gravity">gravity</param><br> </struct><br><br> <struct name="stressBCTop"><br> <param name="Type">StressBC</param><br> <param name="ForceVector">mom_force</param><br> <param name="Swarm">picIntegrationPoints</param><br> <param name="wall">top</param><br> <param name="y_type">HydrostaticTerm</param><br> <param name="y_value">hydrostaticTerm</param><br> </struct> <br><br> <!-- Temperature components --><br> <struct name="temperature"><br> <param name="Type">MeshVariable</param><br> <param name="Rank">Scalar</param><br> <param name="DataType">Double</param><br> <param name="mesh">mesh-linear</param><br> </struct><br> <struct name="temperatureBCs"><br> <param name="Type">CompositeVC</param><br> <param name="Data">mesh-linear</param><br> </struct><br> <struct name="temperatureICs"><br> <param name="Type">CompositeVC</param><br> <param name="Data">mesh-linear</param><br> </struct><br> <struct name="temperatureDofLayout"><br> <param name="Type">DofLayout</param><br> <param name="mesh">mesh-linear</param><br> <list name="BaseVariables"><br> <param>temperature</param><br> </list><br> </struct><br> <struct name="TemperatureField"><br> <param name="Type">FeVariable</param><br> <param name="FEMesh">mesh-linear</param><br> <param name="DofLayout">temperatureDofLayout</param><br> <param name="BC">temperatureBCs</param><br> <param name="IC">temperatureICs</param><br> <param name="LinkedDofInfo">temperatureLinkedDofs</param><br> </struct><br> <struct name="TemperatureGradientsField"><br> <param name="Type">OperatorFeVariable</param><br> <param name="Operator">Gradient</param><br> <param name="FeVariable">TemperatureField</param><br> </struct><br> <!-- Energy Equation --><br><br> <struct name="residual"><br> <param name="Type">ForceVector</param><br> <param name="FeVariable">TemperatureField</param><br> </struct><br> <struct name="massMatrix"><br> <param name="Type">ForceVector</param><br> <param name="FeVariable">TemperatureField</param><br> </struct><br> <struct name="predictorMulticorrector"><br> <param name="Type">AdvDiffMulticorrector</param><br> </struct><br> <struct name="EnergyEqn"><br> <param name="Type">AdvectionDiffusionSLE</param><br> <param name="SLE_Solver">predictorMulticorrector</param><br> <param name="Context">context</param><br> <param name="PhiField">TemperatureField</param><br> <param name="Residual">residual</param><br> <param name="MassMatrix">massMatrix</param><br> <param name="courantFactor">0.25</param><br> </struct><br> <struct name="lumpedMassMatrixForceTerm"><br> <param name="Type">LumpedMassMatrixForceTerm</param><br> <param name="Swarm">gaussSwarm</param><br> <param name="ForceVector">massMatrix</param><br> </struct><br> <struct name="defaultResidualForceTerm"><br> <param name="Type">AdvDiffResidualForceTerm</param><br> <param name="Swarm">gaussSwarm</param><br> <param name="ForceVector">residual</param><br> <param name="ExtraInfo">EnergyEqn</param><br> <param name="VelocityField">VelocityField</param><br> <param name="defaultDiffusivity">defaultDiffusivity</param><br> <param name="UpwindXiFunction">Exact</param><br> </struct><br> <struct name="internalHeatingTerm"><br> <param name="Type">RadiogenicHeatingTerm</param><br> <param name="ForceVector">residual</param><br> <param name="Swarm">picIntegrationPoints</param><br> </struct><span class="Apple-tab-span" style="white-space:pre">        </span><br><br><br> </struct><br> <list name="FieldVariablesToCheckpoint"><br> <param>VelocityField</param><br> <param>PressureField</param><br> </list><br> <param name="maxTimeSteps">60</param><br> <param name="outputEvery">10</param><br> <param name="dumpEvery">10</param><br> <param name="outputPath">output_3</param><br> <param name="dim">2</param><br> <param name="shadowDepth">1</param><br> <param name="minX">0.0f</param><br> <param name="minY">0.0f</param><br> <param name="minZ">0.0f</param><br> <param name="maxX">1320000.0f</param><br> <param name="maxY">440000.0f</param><br> <param name="maxZ">0.5f</param><br> <param name="elementResI">128</param><br> <param name="elementResJ">64</param><br> <param name="elementResK">4</param><br> <param name="allowUnbalancing">True</param><br> <param name="buildNodeNeighbourTbl">True</param><br> <param name="buildElementNeighbourTbl">True</param><br> <param name="particlesPerCell">30</param><br> <param name="seed">13</param><br> <param name="gaussParticlesX">2</param><br> <param name="gaussParticlesY">2</param><br> <param name="gaussParticlesZ">2</param><br> <param name="particlesPerCell">30</param><br> <struct name="EulerDeform"><br> <list name="systems"><br> <struct><br> <param name="mesh">mesh-linear</param><br> <param name="remesher">velocityRemesher</param><br> <param name="velocityField">VelocityField</param><br> <param name="wrapTop">False</param><br> <param name="staticTop">True</param> <br> <param name="staticBottom">True</param><br><span class="Apple-tab-span" style="white-space:pre">        </span><param name="staticRight">True</param><br><span class="Apple-tab-span" style="white-space:pre">        </span><param name="staticLeft">True</param><br> <list name="fields"><br> <struct><br> <param name="field">VelocityField</param><br> <param name="variable">velocity</param><br> </struct><br> <struct><br> <param name="field">PressureField</param><br> <param name="variable">pressure</param><br> </struct><br> </list><br> </struct><br> </list><br> </struct><br> <struct name="velocityBCs"><br> <param name="type">CompositeVC</param><br> <list name="vcList"><br><!-- <struct> --><br><!-- <param name="type">WallVC</param>--><br><!-- <param name="wall">front</param>--><br><!-- <list name="variables">--><br><!-- <struct>--><br><!-- <param name="name">vz</param>--><br><!-- <param name="type">double</param>--><br><!-- <param name="value">0</param>--><br><!-- </struct>--><br><!-- </list>--><br><!-- </struct>--><br><!-- <struct>--><br><!-- <param name="type">WallVC</param>--><br><!-- <param name="wall">back</param>--><br><!-- <list name="variables">--><br><!-- <struct>--><br><!-- <param name="name">vz</param>--><br><!-- <param name="type">double</param>--><br><!-- <param name="value">0</param>--><br><!-- </struct>--><br><!-- </list>--><br><!-- </struct>--><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">left</param><br> <list name="variables"><br> <struct><br> <param name="name">vx</param><br> <param name="type">double</param><br> <param name="value">0.0</param><br> </struct><br> <struct><br> <param name="name">vy</param><br> <param name="type">double</param><br> <param name="value">-1.0e12</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">top</param><br> <list name="variables"><br> <struct><br> <param name="name">vy</param><br> <param name="type">double</param><br> <param name="value">0</param><br> </struct><br> </list><br> </struct><br><!-- <struct>--><br><!-- <param name="type">WallVC</param> --><br><!-- <param name="wall">right</param>--><br><!-- <list name="variables">--><br><!-- <struct> --><br><!-- <param name="name">vy</param>--><br><!-- <param name="type">double</param>--><br><!-- <param name="value">0.0</param>--><br><!-- </struct>--><br><!--<span class="Apple-tab-span" style="white-space:pre">        </span></list>--><br><!-- </struct>--><br><!-- <struct> --><br><!-- <param name="name">vx</param> --><br><!-- <param name="type">double</param> --><br><!-- <param name="value">0.0</param> --><br><!-- </struct> --><br><!-- </list> --><br><!-- </struct> --><br> <!-- <struct> --><br> <!-- <param name="type">MeshShapeVC</param> --><br> <!-- <param name="Shape">fixedCornerShape</param> --><br> <!-- <list name="variables"> --><br> <!-- <struct> --><br> <!-- <param name="name">vy</param> --><br> <!-- <param name="type">double</param> --><br> <!-- <param name="value">0</param> --><br> <!-- </struct> --><br> <!-- </list> --><br> <!-- </struct> --><br><br> </list><br> </struct><br><br><!-- ============== Temperature BCs =============== --><br> <struct name="temperatureBCs"><br> <param name="type">CompositeVC</param><br> <list name="vcList"><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">top</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">double</param><br> <param name="value">273</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">bottom</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">func</param><br> <param name="value">TemperatureProfile</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">left</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">double</param><br> <param name="value">0.0</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">right</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">func</param><br> <param name="value">TemperatureProfile</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">front</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">func</param><br> <param name="value">TemperatureProfile</param><br> </struct><br> </list><br> </struct><br> <struct><br> <param name="type">WallVC</param><br> <param name="wall">back</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">func</param><br> <param name="value">TemperatureProfile</param><br> </struct><br> </list><br> </struct><br> </list><br> </struct><br><br><br><!-- ============== Temperature ICs ===================== --><br> <struct name="temperatureICs"><br> <param name="type">CompositeVC</param><br> <list name="vcList"><br> <struct><br> <param name="type">AllNodesVC</param><br> <list name="variables"><br> <struct><br> <param name="name">temperature</param><br> <param name="type">func</param><br> <param name="value">TemperatureProfile</param><br> </struct><br> </list><br> </struct><br> </list><br> </struct><br><br><br> <param name="TemperatureProfileTop">273.0</param><br> <param name="TemperatureProfileLinearCoefficient">0.0135</param><br> <param name="TemperatureProfileExponentialCoefficient1">0</param><br> <param name="TemperatureProfileExponentialCoefficient2">0</param><br> <param name="TemperatureProfileMax">1623.0</param><br> <param name="TemperatureProfileH0">100000.0</param><br><br><br><!-- <param name="checkpointEvery">1</param>--><br> <param name="dtFactor">0.1</param><br> <param name="checkpointEvery">10</param><br> <param name="gravity">10.0</param><br> <param name="journal.info">True</param><br> <param name="journal.debug">True</param><br> <param name="journal-level.info">2</param><br> <param name="journal-level.debug">2</param><br><br></StGermainData><br><span><vertical_slab.png></span>_______________________________________________<br>CIG-LONG mailing list<br>CIG-LONG@geodynamics.org<br>http://geodynamics.org/cgi-bin/mailman/listinfo/cig-long<br></blockquote></div><br><div apple-content-edited="true"> <span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: 'Lucida Grande'; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>-----------------------------</div><div>Associate Professor, U.C. Davis</div><div>Department of Geology/KeckCAVES</div><div>2129 Earth & Physical Sciences Bldg</div><div>Davis, CA 95616</div><div>-----------------</div><div><a href="mailto:mibillen@ucdavis.edu">mibillen@ucdavis.edu</a></div><div>(530) 754-5696</div><div><b>--------------------------</b></div></div></div></span></div></span></div></span></div></span> </div><br></div></body></html>