!!!!!!Modified from Stokes_diagnostic.sif !echo on Header !CHECK KEYWORDS Warn Mesh DB "." "testglacier2" Include Path "" Results Directory "" End Simulation Max Output Level = 4 Coordinate System = "Cartesian 2D" Coordinate Mapping(3) = 1 2 3 Simulation Type = "Transient" Steady State Max Iterations = 1 Timestepping Method = "BDF" BDF Order = 1 Timestep Sizes = 1.0 Timestep Intervals = 5 Output Intervals = 1 Output File = "scarp_test.result" Post File = "scarp_test.vtu" ! use .ep suffix for leagcy format Initialize Dirichlet Conditions = Logical False End Constants Gravity(4) = 0 -1 0 3.711 Stefan Boltzmann = 5.67e-08 End Body 1 Name = "Glacier" Body Force = 1 Equation = 1 Material = 1 Initial Condition = 1 End Body 2 Name = "surface" Equation = 2 Material = 2 Body Force = 2 Initial Condition = 2 End Equation 1 Name = "Equation1" Convection = "computed" Flow Solution Name = String "Flow Solution" Active Solvers(5) = 1 2 3 4 6 End Equation 2 Active Solvers(1) = 5 Flow Solution Name = String "Flow Solution" Convection = Computed End Initial Condition 1 Velocity 1 = 0.0 Velocity 2 = 0.0 Pressure = 0.0 Depth = Real 0.0 End Initial Condition 2 Zs = Equals Coordinate 2 RefZs = Equals Coordinate 2 Mesh Update 3 = Real 0.0 Mesh Update 2 = Real 0.0 Mesh Update 1 = Real 0.0 End ! maps DEM's to originally rectangular mesh ! see Top and Bottom Surface in BC's Solver 1 Exec Solver = "before Simulation" Equation = "MapCoordinate" Procedure = "StructuredMeshMapper" "StructuredMeshMapper" Active Coordinate = Integer 2 ! the mesh-update is y-direction ! For time being this is currently externally allocated Mesh Velocity Variable = String "Mesh Velocity 2" ! The 1st value is special as the mesh velocity could be unrelistically high Mesh Velocity First Zero = Logical True Dot Product Tolerance = Real 0.01 End ! FlowDepth on generally unstructured mesh ! (could be replaced by strucuterd version) Solver 2 Equation = "Flowdepth" Exec Solver = "Before Simulation" Procedure = File "ElmerIceSolvers" "FlowDepthSolver" Variable = String "Depth" Variable DOFs = 1 Linear System Solver = "Direct" Linear System Direct Method = "UMFPACK" Linear System Max Iterations = 200 Linear System Convergence Tolerance = 1.0E-09 Linear System Abort Not Converged = False Linear System Preconditioning = "ILU0" Linear System Residual Output = 1 Gradient = Real -1.0E00 Exported Variable 1 = -dofs 3 "Mesh Velocity" Calc Free Surface = Logical True ! the name for the exported (if not existing) added variable ! the gradients will be stored in variables with the base ! name given and "Grad1" and (in 3 dimensions) "Grad2" added, ! so in our case "FreeSurfGrad1" and "FreeSurfGrad2" ! again, if those variables did not exist, they will be ! automatically created !----------------------------------------------------------- Freesurf Name = String "FreeSurf" End Solver 3 ! Exec Solver = "Never" # uncommenting would switch this off Equation = "Navier-Stokes" Optimize Bandwidth = Logical True ! direct solver Linear System Solver = Direct Linear System Direct Method = "UMFPACK" Linear System Max Iterations = 5000 Linear System Convergence Tolerance = 1.0E-06 Linear System Abort Not Converged = False Linear System Preconditioning = "ILU1" Linear System Residual Output = 1 Steady State Convergence Tolerance = 1.0E-05 Stabilization Method = Stabilized Nonlinear System Convergence Tolerance = 1.0E-04 Nonlinear System Convergence Measure = Solution Nonlinear System Max Iterations = 50 Nonlinear System Newton After Iterations = 3 Nonlinear System Newton After Tolerance = 1.0E-01 ! Nonlinear System Relaxation Factor = 0.75 End Solver 4 Equation = Sij Procedure = "ElmerIceSolvers" "ComputeDevStress" Variable = -nooutput "Sij" Variable DOFs = 1 Exported Variable 1 = -dofs 4 Stress Stress Variable Name = String "Stress" Flow Solver Name = String "Flow Solution" Linear System Solver = Direct Linear System Direct Method = umfpack End Solver 5 Exec Solver = "After TimeStep" Equation = String "Free Surface Evolution" Variable = "FS" Variable DOFs = 1 Procedure = "FreeSurfaceSolver" "FreeSurfaceSolver" Apply Dirichlet = Logical True Linear System Solver = Iterative Linear System Iterative Method = BiCGStab Linear System Max Iterations = 100 Linear System Preconditioning = ILU1 Linear System Convergence Tolerance = 1.0e-8 Nonlinear System Max Iterations = 100 Nonlinear System Min Iterations = 2 Nonlinear System Convergence Tolerance = 1.0e-6 Steady State Convergence Tolerance = 1.0e-4 Stabilization Method = Bubbles Flow Solution Name = String "Flow Solution" Exported Variable 1 = FS Residual Exported Variable 1 DOFS = 1 Exported Variable 2 = ReferenceFS Exported Variable 2 DOFS = 1 End Solver 6 Exec Solver = "after timestep" Equation = "MapCoordinate" Procedure = "StructuredMeshMapper" "StructuredMeshMapper" Active Coordinate = Integer 2 ! the mesh-update is y-direction ! For time being this is currently externally allocated Mesh Velocity Variable = String "Mesh Velocity 2" ! The 1st value is special as the mesh velocity could be unrelistically high Mesh Velocity First Zero = Logical True ! Top Surface Variable = String "Zs" Dot Product Tolerance = Real 0.01 End ! we use m-yr-MPa system 1 yr = 31556926.0 sec Material 1 Name = "ice-ice-baby" Density = Real $910.0*1.0E-06*(31556926.0)^(-2.0) Viscosity Model = String "Glen" ! Viscosity has to be set to a dummy value ! to avoid warning output from Elmer Viscosity = Real 1.0 Glen Exponent = Real 3.0 Critical Shear Rate = Real 1.0e-10 ! Rate factors (Paterson value in MPa^-3a^-1) Rate Factor 1 = Real 1.258e13 Rate Factor 2 = Real 6.046e28 ! these are in SI units - no problem, as long as ! the gas constant also is Activation Energy 1 = Real 60e3 Activation Energy 2 = Real 139e3 Glen Enhancement Factor = Real 1.0 ! the variable taken to evaluate the Arrhenius law ! in general this should be the temperature relative ! to pressure melting point. The suggestion below plugs ! in the correct value obtained with TemperateIceSolver ! Temperature Field Variable = String "Temp Homologous" ! the temperature to switch between the ! two regimes in the flow law Limit Temperature = Real 0.0 ! In case there is no temperature variable (which here is the case) Constant Temperature = Real MATC "230-273.15" Cauchy = Logical False End Material 2 Density = Real $910.0*1.0E-06*(31556926.0)^(-2.0) Min Zs = Real 0.0e0 Max Zs = Real 1000 End Body Force 1 Name = "BodyForce1" Heat Source = 1 Flow BodyForce 1 = Real 0.0 Flow BodyForce 2 = Real -3.711e15 !MPa - a - m End Body Force 2 Zs Accumulation Flux 1 = Real 0.0e0 Zs Accumulation Flux 2 = Real 0.0e0 End Boundary Condition 1 Name = "bedrock" Target Boundaries = 1 Compute Normals = Logical True ! include the bedrock DEM, which has two colums Bottom Surface = Variable Coordinate 1 Real cubic include "groundwide.dat" End Velocity 1 = Real 0.0e0 Velocity 2 = Real 0.0e0 End Boundary Condition 2 Name = "sides" Target Boundaries(2) = 3 4 ! combine left and right boundary Velocity 1 = 0.0 End Boundary Condition 3 Name = "surface" Target Boundaries = 2 Body ID = 2 ! include the surface DEM, which has two colums Top Surface = Variable Coordinate 1 Real cubic include "icescarp70degground.dat" End ! Top Surface = Equals Zs Depth = Real 0.0 End