Check Keywords "Warn" Header Mesh DB "." "mesh2D" End !************************************************************** ! Constant !************************************************************** $gravity = 9.7696e15 $yearinsec = 365.25*24*60*60 $rhow = 1000.0/(1.0e6*yearinsec^2) $rhoi = 910.0/(1.0e6*yearinsec^2) $name = "steady" !************************************************************** ! Simulation !************************************************************** Simulation Simulation Type = steady Steady State Min Iterations = 1 Steady State Max Iterations = 1 Output Intervals(1) = 1 !Post File = case.ep Max Output Level = 42 Extruded Mesh Levels = Integer 10 Restart Before Initial Conditions = Logical True Restart File = "gl.result" Restart Position = 0 Output File = $name".result" Post File = $name".ep" Initialize Dirichlet Conditions = Logical False End !************************************************************** ! Body !************************************************************** Body 1!Main ice body Name = "glacier" Equation = 1 Material = 1 Initial Condition = 1 Body Force = 1 End !*************************************************** !body force !*************************************************** Body Force 1 Flow BodyForce 1 = Real 0.0 Flow BodyForce 2 = Real 0.0 Flow BodyForce 3 = Real -$gravity End !*************************************************** !Initial condition !*************************************************** Initial Condition 1 Pressure = Real 0.0 End !*************************************************** !Material !*************************************************** Material 1 Density = Real $rhoi ! Glen's flow law (using Glen) !----------------- ! viscosity stuff !---------------- Viscosity Model = String "Glen" Viscosity = Real 2.5 !---- Parametres pour Geln Viscosityinverse.sif !the temperature to switch between the two regimes in the flow law !Temperature Field Variable = String "Temp Homologous" Glen Exponent = Real 3.0 Glen Enhancement Factor = Real 2.5e00 Limit Temperature = Real -10.0e00 Critical Shear Rate = Real 1.0e-12 ! Rate factors (Paterson value in MPa^-3a^-1) Rate Factor 1 = Real 1.257543501100000e+13 !$3.985e-13*spy*1.0e18 Activation Energy 1 = Real 60.0e03 Rate Factor 2 = Real 6.046307021600000e+28 !$1.916e03*spy*1.0e18 Activation Energy 2 = Real 139.0e03 Sea level = Real 0.0 End !**************************************************************** !Equation !**************************************************************** Equation 1 !main ice Active Solvers(3) = 1 2 3 Convection = Computed Flow Solution Name = String "Flow Solution" NS Convect = Logical False End !*************************************************** !Solver !*************************************************** !!!! Navier-Stokes Solution Solver 1 Exec Solver = Never Equation = "Navier-Stokes" Stabilize = logical True flow model = Stokes !mandatory to save bulk stiffness matrix calculate loads = Logical True Linear System Solver = Direct Linear System Direct Method = mumps mumps percentage increase working space = integer 200 !Linear System Solver = Iterative ! Linear System Iterative Method = GMRES ! Linear System GMRES Restart = 100 ! Linear System Preconditioning= ILU0 ! Linear System Convergence Tolerance= 1.0e-08 ! Linear System Max Iterations = 1000 !system self adjoint if Newton is used for the last iterations Nonlinear System Max Iterations = Integer 100!100 Nonlinear System Convergence Tolerance = Real 1.0e-7 Nonlinear System Newton After Iterations = Integer 10!15 Nonlinear System Newton After Tolerance = Real 1.0e-04 Nonlinear System Relaxation Factor = Real 1.0 Nonlinear System Reset Newton = Logical True Steady State Convergence Tolerance = Real 1.0e-7 End Solver 2 Exec Solver = Before simulation Equation = "ScatteredInter" procedure = "Scattered2DDataInterpolator" "Scattered2DDataInterpolator" CheckNaN = Logical True ! Default True; check is interpolation method gives NaN ! By default replace NaN by nearest available value ! We can replace NaNs by Real value with following flag !Replace NaN by = Real -9999999 !!!!! NNI or linear (nn-c library) ! Default Sibson interpolation Variable 1 = String "beta" Variable 1 data file = File "beta.dat" Variable 1 W = Real 0. ! W restricts extrapolation by assigning minimal allowed ! weight for a vertex (normally "-1" or so; lower ! values correspond to lower reliability; "0" means ! no extrapolation) ! Default W=-HUGE(RealNumber);i.e. extrapolation allowed !Variable 1 method = String "ns" ! method Non-Sibsonian interpolation (nn-c); W can be applied here too Variable 1 method = String "li" ! method linear interpolation (nn-c); W no effect here !Variable 1 method = String "cs" ! method cubic spline (csa) Variable 1 nppc = integer 5 !set the average number of points per cell (default = 5, !works best for uniform data. Decrease to get smaller ! cells or increase to get larger cells) Variable 1 k = integer 140 ! set the spline sensitivity (default = 140, reduce to get ! smoother results) End Solver 3 Exec Solver = "After Timestep" Equation = "ResultOutput" Procedure = File "ResultOutputSolve" "ResultOutputSolver" Output File Name = file "Output_$name"" Vtu Format = logical true Binary Output = True Single Precision = True Save Geometry Ids = Logical True Vector Field 1 = "Velocity" Scalar Field 1 = speed Scalar Field 2 = surface_ele Scalar Field 3 = beta Scalar Field 4 = Height Scalar Field 5 = Depth End !************************************************************* !boundary condition !************************************************************* Boundary Condition 1 Name = "front" Target Boundaries(1) = 1 !sea pressure on the grounding line External Pressure = Variable Coordinate 3 Real Procedure "ElmerIceUSF" "SeaPressure" Compute Sea Pressure = Logical True Free Moving = True End Boundary Condition 2 Name = "side" Target Boundaries(1) = 2 Normal-Tangential Velocity = Logical True Mass Consistent Normals = Logical True Normal-Tangential Velocity Rotate = Logical False Slip Coefficient 1 = Real 10.0 Slip Coefficient 2 = Real 0.00001 Slip Coefficient 3 = Real 0.00001 External Pressure = Variable depth Real MATC "1.0 * rhoi * gravity * tx" End Boundary Condition 3 Name = "bed" Target Boundaries = 3 !Flow Force BC = Logical True Mass Consistent Normals = Logical True Normal-Tangential Velocity = Logical True Velocity 1 = Real 0.0 Slip Coefficient 2 = Variable Beta Real MATC "10^tx(0)" Slip Coefficient 3 = Variable Beta Real MATC "10^tx(0)" height = Real 0.0 End Boundary Condition 4 Name = "surface" Target Boundaries = 4 Depth = Real 0.0 !Top Surface = Equals FSnew End