! Simulation title $title = "1997-2007" ! total years of simulation $totalYears = 10 Header Mesh DB "." “glacier-mesh” End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Constants Water Density = Real $ 1000.0 ! DEM files SurfaceDEM = String "../../mesh/extrusion/dem.xyz" BedrockDEM = String "../../mesh/extrusion/bedrock_xyz_utm32n.xyz" ! Change Arrhenius factor, balance and basal velocity after ! this time after each the time timesteps SwitchPeriod = Real 6 !!! Arrhenius factor as function of T ! (Patterson interpolated value) WinterArrheniusFactor = Real $ 2.4*10^-24 !WinterArrheniusFactor = Real $ 6.4*10^-25 ! (Patterson T=0°C - 24*10^-25 s-1*Pa-3) 75.74 MPa-3*a-1 SummerArrheniusFactor = Real $ 2.4*10^-24 !!! Basal kinematic condition ! No sliding before SwitchPeriod BasalWinterVelocity = Real $ 0.0 ! Sliding on the bedrock BasalSummerVelocity = Real $ 0.0 !BasalSummerVelocity = Real 0.0 ! winter accumulation file WinterAccumulationFile = String balance_winter.dat ! summer accumulation file SummerAccumulationFile = String balance_summer.dat ! min ice thickness MinIceThickness = Real 0.5 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Simulation Coordinate System = "Cartesian 3D" Simulation Type = Transient Timestepping Method = "bdf" BDF Order = 1 TimeStep Intervals = Integer $totalYears*12 Timestep Sizes = Real $ 30*86400 Output Intervals = 1 Steady State Min Iterations = 1 Steady State Max Iterations = 10 Output File = "$title".result" Post File = "$title".ep" Max Output Level = 3 Extruded Mesh Levels = Integer 10 !Restart File = "$title".result" !Restart Position = 4 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Body 1 Equation = 1 Material = 1 Body Force = 1 Initial condition = 1 End Body 2 Equation = 2 Body Force = 2 Material = 1 Initial Condition = 2 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Body Force 1 Flow BodyForce 1 = Real 0.0 Flow BodyForce 2 = Real 0.0 Flow Bodyforce 3 = Real $ -9.81 End Body Force 2 !Zs Accumulation = Real 0.0 !Zs Accumulation = Real $ -0.5 / (86400*60) 3 metri anno Zs Accumulation = Variable Coordinate 1 Real Procedure "./utils_eros_fede" "getBalance" End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Initial Condition 1 Pressure = Real 0.00 Velocity 1 = Real 0.00 Velocity 2 = Real 0.00 Velocity 3 = Real 0.00 End Initial Condition 2 Zs = Variable Coordinate 1 Real Procedure "./utils_eros_fede" "getSurfaceElevation" End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Boundary Condition 1 ! lateral side Target Boundaries = 1 Save Scalars = Logical True Velocity 1 = real 0.0 Velocity 2 = real 0.0 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Boundary Condition 2 ! bedrock !Normal-Tangential Velocity = Logical True !Velocity 1 = Real 0.0 !Velocity 2 = real Procedure "./utils_eros_fede" "getBottomVelocity" !gets value previously set in BasalSummerVelocity or BasalWinterVelocity !Velocity 3 = real Procedure "./utils_eros_fede" "getBottomVelocity" Flow Force BC = Logical True Normal-Tangential Velocity = Logical True Velocity 1 = Real 0.0 Slip Coefficient 2 = Real 0.1 Slip Coefficient 3 = Real 0.1 Bottom Surface = Variable Coordinate 1 Real Procedure "./utils_eros_fede" "getBedrockElevation" Save Scalars = Logical True End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Boundary Condition 3 ! upper surface Body Id = 2 Save Scalars = Logical True End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Material 1 Density = Real $ 917.0 Viscosity Model = String "glen" Viscosity = 1.0 ! Avoid warning output. This value will not be used. ! (Patterson T=0°C - 24*10^-25 s-1*Pa-3) 75.74 MPa-3*a-1 Set Arrhenius Factor = Logical True Arrhenius Factor = Real Procedure "./utils_eros_fede" "getArrheniusFactor" Critical Shear Rate = Real $ 1.0e-7 ! Bed condition Min Zs = Variable Coordinate 1 Real Procedure "./utils_eros_fede" "getMinSurfaceElevation" Max Zs = Real +1.0e10 End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Solver 1 Equation = "MapCoordinate" Procedure = "StructuredMeshMapper" "StructuredMeshMapper" Active Coordinate = Integer 3 Mesh Velocity Variable = String "dSdt" Mesh Update Variable = String "dS" Mesh Velocity First Zero = Logical True Top Surface Variable Name = String "Zs" Displacement Mode = Logical False Correct Surface = Logical True Minimum Height = Real 0.5 End Solver 2 Equation = "HeightDepth" Procedure = "StructuredProjectToPlane" "StructuredProjectToPlane" Active Coordinate = Integer 3 Operator 1 = depth Operator 2 = height Operator 3 = thickness End Solver 3 Equation = Navier-Stokes Stabilization Method = String Stabilized Flow Model = Stokes Exported Variable 1 = -dofs 1 "dSdt" Exported Variable 2 = -dofs 1 "dS" Linear System Solver = Direct Linear System Direct Method = umfpack Nonlinear System Max Iterations = 50 Nonlinear System Convergence Tolerance = 1.0e-5 ! first 5 iterations or convergence < 1.e-2 before moving from Picard to Newton Nonlinear System Newton After Iterations = 5 Nonlinear System Newton After Tolerance = 1.0e-02 !Nonlinear System Relaxation Factor = 1.00 ! reset newton to false each new time step Nonlinear System Reset Newton = Logical True Steady State Convergence Tolerance = Real 1.0e-3 End Solver 4 Equation = Sij Procedure = "ElmerIceSolvers" "ComputeDevStress" Variable = -nooutput "Sij" Variable DOFs = 1 Exported Variable 1 = -dofs 6 Stress Stress Variable Name = String "Stress" Flow Solver Name = String "Flow Solution" Linear System Solver = Direct Linear System Direct Method = umfpack End Solver 5 Equation = "EigenStresses" Procedure = "ElmerIceSolvers" "ComputeEigenValues" Variable = -nooutput dumy Variable DOFs = 1 ! The 3 eigenvalues Exported Variable 1 = EigenStress Exported Variable 1 DOFS = 3 ! The 3 eigenvectors (Option) Exported Variable 2 = EigenVector1 Exported Variable 2 DOFS = 3 Exported Variable 3 = EigenVector2 Exported Variable 3 DOFS = 3 Exported Variable 4 = EigenVector3 Exported Variable 4 DOFS = 3 End Solver 6 Equation = "Free Surface Top" Variable = String "Zs" Variable DOFs = 1 Exported Variable 1 = String "Zs Residual" Exported Variable 1 DOFs = 1 Procedure = "FreeSurfaceSolver" "FreeSurfaceSolver" Before Linsolve = "EliminateDirichlet" "EliminateDirichlet" Linear System Solver = Iterative Linear System Max Iterations = 1500 Linear System Iterative Method = BiCGStab Linear System Preconditioning = ILU1 Linear System Convergence Tolerance = Real 1.0e-5 Linear System Abort Not Converged = True Linear System Residual Output = 3 Nonlinear System Max Iterations = 20 Nonlinear System Convergence Tolerance = 1.0e-6 Nonlinear System Relaxation Factor = 0.5 Steady State Convergence Tolerance = 1.0e-03 Stabilization Method = Bubbles Apply Dirichlet = Logical False !No limiter needed for the upper surface ! How much the free surface is relaxed Relaxation Factor = Real 1.00 ! Is there a maximum step-size for the displacement use/or not accumulation Use Accumulation = Logical True ! take accumulation to be given normal to surface/as vector Normal Flux = Logical False End Solver 7 Equation = SaveScalars Procedure = "SaveData" "SaveScalars" Exec Solver = After Timestep File Append = Logical True Filename = "$title".SaveScalars.dat" ! Compute convective flux. Add coefficient C (density) to ! compute the mass discharge Variable 1 = String "Flow Solution" Coefficient 1 = String "Density" Operator 1 = String "Convective Flux" Variable 2 = String "Flow Solution" Operator 2 = String "Volume" End Solver 8 Exec Solver = After TimeStep Exec Interval = 1 Equation = "result output" Procedure = "ResultOutputSolve" "ResultOutputSolver" Output File Name = File "$title".vtu" Output Format = String vtu End !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Equation 1 Active Solvers(7) = 1 2 3 4 5 7 8 End Equation 2 Active Solvers(1) = 6 Flow Solution Name = String "Flow Solution" Convection = String Computed End