total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 3.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.94752e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 20/100 time [yr]: 4.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.51729e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.113057 < 1 %
Convergence criterion: norm(du)/norm(u) 1.06219 < 10 %
Convergence criterion: max(du) 5.36896e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.02278e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 21/100 time [yr]: 4.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 4.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.83565e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.11096 < 1 %
Convergence criterion: norm(du)/norm(u) 1.01971 < 10 %
Convergence criterion: max(du) 5.02071e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 4.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.67186e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 22/100 time [yr]: 4.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 4.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.56222e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.109028 < 1 %
Convergence criterion: norm(du)/norm(u) 0.980716 < 10 %
Convergence criterion: max(du) 4.71242e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 4.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.55355e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 23/100 time [yr]: 4.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 4.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.80966e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.107249 < 1 %
Convergence criterion: norm(du)/norm(u) 0.944488 < 10 %
Convergence criterion: max(du) 4.43423e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 4.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.86739e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 24/100 time [yr]: 4.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 4.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.85174e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.105604 < 1 %
Convergence criterion: norm(du)/norm(u) 0.910484 < 10 %
Convergence criterion: max(du) 4.18077e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 4.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.86171e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 25/100 time [yr]: 5.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 5
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.58625e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.104068 < 1 %
Convergence criterion: norm(du)/norm(u) 0.878442 < 10 %
Convergence criterion: max(du) 3.94913e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 5
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81196e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 26/100 time [yr]: 5.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 5.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.59457e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.102627 < 1 %
Convergence criterion: norm(du)/norm(u) 0.84817 < 10 %
Convergence criterion: max(du) 3.73747e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 5.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.31159e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 27/100 time [yr]: 5.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 5.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.9639e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.101269 < 1 %
Convergence criterion: norm(du)/norm(u) 0.819596 < 10 %
Convergence criterion: max(du) 3.54432e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 5.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.97032e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 28/100 time [yr]: 5.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 5.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.68091e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.099986 < 1 %
Convergence criterion: norm(du)/norm(u) 0.792599 < 10 %
Convergence criterion: max(du) 3.36823e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 5.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.54016e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 29/100 time [yr]: 5.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 5.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.63083e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0987683 < 1 %
Convergence criterion: norm(du)/norm(u) 0.766973 < 10 %
Convergence criterion: max(du) 3.2077e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 5.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.66299e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 30/100 time [yr]: 6.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.87659e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.097609 < 1 %
Convergence criterion: norm(du)/norm(u) 0.742606 < 10 %
Convergence criterion: max(du) 3.06117e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.82249e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 31/100 time [yr]: 6.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 6.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71514e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0965021 < 1 %
Convergence criterion: norm(du)/norm(u) 0.719506 < 10 %
Convergence criterion: max(du) 2.92896e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 6.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.18904e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 32/100 time [yr]: 6.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 6.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.74799e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0954442 < 1 %
Convergence criterion: norm(du)/norm(u) 0.697603 < 10 %
Convergence criterion: max(du) 2.81174e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 6.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.83118e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 33/100 time [yr]: 6.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 6.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.5928e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0944289 < 1 %
Convergence criterion: norm(du)/norm(u) 0.676853 < 10 %
Convergence criterion: max(du) 2.70888e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 6.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.84957e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 34/100 time [yr]: 6.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 6.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.97756e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0934533 < 1 %
Convergence criterion: norm(du)/norm(u) 0.657234 < 10 %
Convergence criterion: max(du) 2.61777e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 6.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.04121e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 35/100 time [yr]: 7.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 7
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.75162e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0925141 < 1 %
Convergence criterion: norm(du)/norm(u) 0.63862 < 10 %
Convergence criterion: max(du) 2.53451e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 7
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.72367e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 36/100 time [yr]: 7.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 7.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.63812e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0916078 < 1 %
Convergence criterion: norm(du)/norm(u) 0.620847 < 10 %
Convergence criterion: max(du) 2.45591e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 7.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.86689e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 37/100 time [yr]: 7.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 7.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72928e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0907312 < 1 %
Convergence criterion: norm(du)/norm(u) 0.603774 < 10 %
Convergence criterion: max(du) 2.38007e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 7.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.99077e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 38/100 time [yr]: 7.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 7.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72189e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0898818 < 1 %
Convergence criterion: norm(du)/norm(u) 0.587275 < 10 %
Convergence criterion: max(du) 2.30628e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 7.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.64158e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 39/100 time [yr]: 7.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 7.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.7118e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0890573 < 1 %
Convergence criterion: norm(du)/norm(u) 0.571272 < 10 %
Convergence criterion: max(du) 2.23457e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 7.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.58669e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 40/100 time [yr]: 8.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72467e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.088256 < 1 %
Convergence criterion: norm(du)/norm(u) 0.555722 < 10 %
Convergence criterion: max(du) 2.16533e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69748e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 41/100 time [yr]: 8.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 8.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.74636e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0874763 < 1 %
Convergence criterion: norm(du)/norm(u) 0.540621 < 10 %
Convergence criterion: max(du) 2.09897e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 8.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.83434e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 42/100 time [yr]: 8.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 8.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.67431e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0867169 < 1 %
Convergence criterion: norm(du)/norm(u) 0.525963 < 10 %
Convergence criterion: max(du) 2.03579e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 8.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.80343e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 43/100 time [yr]: 8.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 8.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.75385e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0859763 < 1 %
Convergence criterion: norm(du)/norm(u) 0.511744 < 10 %
Convergence criterion: max(du) 1.97592e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 8.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.07689e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 44/100 time [yr]: 8.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 8.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72686e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0852535 < 1 %
Convergence criterion: norm(du)/norm(u) 0.49796 < 10 %
Convergence criterion: max(du) 1.91932e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 8.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.76805e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 45/100 time [yr]: 9.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 9
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.56113e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0845475 < 1 %
Convergence criterion: norm(du)/norm(u) 0.484607 < 10 %
Convergence criterion: max(du) 1.86589e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 9
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96218e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 46/100 time [yr]: 9.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 9.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71756e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0838577 < 1 %
Convergence criterion: norm(du)/norm(u) 0.471685 < 10 %
Convergence criterion: max(du) 1.81543e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 9.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.56477e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 47/100 time [yr]: 9.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 9.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.69827e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0831831 < 1 %
Convergence criterion: norm(du)/norm(u) 0.459161 < 10 %
Convergence criterion: max(du) 1.76773e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 9.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69552e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 48/100 time [yr]: 9.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 9.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.64417e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.082523 < 1 %
Convergence criterion: norm(du)/norm(u) 0.44707 < 10 %
Convergence criterion: max(du) 1.72258e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 9.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.99405e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 49/100 time [yr]: 9.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 9.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.73017e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0818764 < 1 %
Convergence criterion: norm(du)/norm(u) 0.435391 < 10 %
Convergence criterion: max(du) 1.67977e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 9.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.07658e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 50/100 time [yr]: 10.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 10
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.73979e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0812437 < 1 %
Convergence criterion: norm(du)/norm(u) 0.424117 < 10 %
Convergence criterion: max(du) 1.63912e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 10
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.15705e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 51/100 time [yr]: 10.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 10.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.59819e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0806238 < 1 %
Convergence criterion: norm(du)/norm(u) 0.413234 < 10 %
Convergence criterion: max(du) 1.60048e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 10.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.08663e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 52/100 time [yr]: 10.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 10.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.55582e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0800205 < 1 %
Convergence criterion: norm(du)/norm(u) 0.397316 < 10 %
Convergence criterion: max(du) 1.552e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 10.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.80841e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 53/100 time [yr]: 10.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 10.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.5006e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0794292 < 1 %
Convergence criterion: norm(du)/norm(u) 0.381177 < 10 %
Convergence criterion: max(du) 1.50189e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 10.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.97641e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 54/100 time [yr]: 10.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 10.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.79874e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0788494 < 1 %
Convergence criterion: norm(du)/norm(u) 0.365536 < 10 %
Convergence criterion: max(du) 1.45386e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 10.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.74228e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 55/100 time [yr]: 11.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 11
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.60333e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0782806 < 1 %
Convergence criterion: norm(du)/norm(u) 0.350543 < 10 %
Convergence criterion: max(du) 1.40887e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 11
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.44546e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 56/100 time [yr]: 11.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 11.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.70797e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0777227 < 1 %
Convergence criterion: norm(du)/norm(u) 0.336218 < 10 %
Convergence criterion: max(du) 1.36692e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 11.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.64262e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 57/100 time [yr]: 11.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 11.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.68942e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0771752 < 1 %
Convergence criterion: norm(du)/norm(u) 0.322546 < 10 %
Convergence criterion: max(du) 1.32775e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 11.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.71963e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 58/100 time [yr]: 11.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 11.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.58694e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0766378 < 1 %
Convergence criterion: norm(du)/norm(u) 0.309493 < 10 %
Convergence criterion: max(du) 1.291e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 11.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.68995e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 59/100 time [yr]: 11.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 11.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.70303e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0761106 < 1 %
Convergence criterion: norm(du)/norm(u) 0.29704 < 10 %
Convergence criterion: max(du) 1.25632e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 11.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.80298e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 60/100 time [yr]: 12.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 12
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71223e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0755928 < 1 %
Convergence criterion: norm(du)/norm(u) 0.285065 < 10 %
Convergence criterion: max(du) 1.22337e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 12
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.53757e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 61/100 time [yr]: 12.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 12.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.67325e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0750838 < 1 %
Convergence criterion: norm(du)/norm(u) 0.273806 < 10 %
Convergence criterion: max(du) 1.19189e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 12.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.83086e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 62/100 time [yr]: 12.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 12.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.60244e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0745836 < 1 %
Convergence criterion: norm(du)/norm(u) 0.263097 < 10 %
Convergence criterion: max(du) 1.16165e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 12.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.4489e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 63/100 time [yr]: 12.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 12.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.70679e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0740922 < 1 %
Convergence criterion: norm(du)/norm(u) 0.253041 < 10 %
Convergence criterion: max(du) 1.13249e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 12.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.68089e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 64/100 time [yr]: 12.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 12.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.76898e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0736093 < 1 %
Convergence criterion: norm(du)/norm(u) 0.243569 < 10 %
Convergence criterion: max(du) 1.10429e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 12.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.61613e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 65/100 time [yr]: 13.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 13
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71003e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0731343 < 1 %
Convergence criterion: norm(du)/norm(u) 0.234663 < 10 %
Convergence criterion: max(du) 1.07697e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 13
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.84941e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 66/100 time [yr]: 13.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 13.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.5544e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.072667 < 1 %
Convergence criterion: norm(du)/norm(u) 0.22627 < 10 %
Convergence criterion: max(du) 1.05046e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 13.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.6062e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 67/100 time [yr]: 13.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 13.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.78849e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0722078 < 1 %
Convergence criterion: norm(du)/norm(u) 0.218445 < 10 %
Convergence criterion: max(du) 1.02472e-06
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 13.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.61362e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 68/100 time [yr]: 13.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 13.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.68389e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0717558 < 1 %
Convergence criterion: norm(du)/norm(u) 0.21117 < 10 %
Convergence criterion: max(du) 9.99696e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 13.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69417e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 69/100 time [yr]: 13.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 13.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.81617e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0713112 < 1 %
Convergence criterion: norm(du)/norm(u) 0.204457 < 10 %
Convergence criterion: max(du) 9.75361e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 13.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.67264e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 70/100 time [yr]: 14.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 14
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72209e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0708733 < 1 %
Convergence criterion: norm(du)/norm(u) 0.198238 < 10 %
Convergence criterion: max(du) 9.5168e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 14
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.55908e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 71/100 time [yr]: 14.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 14.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.01056e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0704422 < 1 %
Convergence criterion: norm(du)/norm(u) 0.192454 < 10 %
Convergence criterion: max(du) 9.28272e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 14.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69649e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 72/100 time [yr]: 14.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 14.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.77195e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0700178 < 1 %
Convergence criterion: norm(du)/norm(u) 0.18729 < 10 %
Convergence criterion: max(du) 9.05091e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 14.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.56669e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 73/100 time [yr]: 14.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 14.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.62856e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0696 < 1 %
Convergence criterion: norm(du)/norm(u) 0.182675 < 10 %
Convergence criterion: max(du) 8.82239e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 14.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.68257e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 74/100 time [yr]: 14.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 14.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.73803e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0691883 < 1 %
Convergence criterion: norm(du)/norm(u) 0.178566 < 10 %
Convergence criterion: max(du) 8.59861e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 14.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.58075e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 75/100 time [yr]: 15.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 15
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.73997e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0687826 < 1 %
Convergence criterion: norm(du)/norm(u) 0.17494 < 10 %
Convergence criterion: max(du) 8.38073e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 15
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.27078e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 76/100 time [yr]: 15.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 15.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.54708e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0683825 < 1 %
Convergence criterion: norm(du)/norm(u) 0.171801 < 10 %
Convergence criterion: max(du) 8.16926e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 15.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.68689e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 77/100 time [yr]: 15.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 15.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.74154e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0679884 < 1 %
Convergence criterion: norm(du)/norm(u) 0.168896 < 10 %
Convergence criterion: max(du) 7.96414e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 15.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.8653e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 78/100 time [yr]: 15.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 15.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.61721e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0675996 < 1 %
Convergence criterion: norm(du)/norm(u) 0.166233 < 10 %
Convergence criterion: max(du) 7.76498e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 15.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.55188e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 79/100 time [yr]: 15.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 15.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.68542e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0672163 < 1 %
Convergence criterion: norm(du)/norm(u) 0.164213 < 10 %
Convergence criterion: max(du) 7.57129e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 15.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.46967e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 80/100 time [yr]: 16.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 16
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.74069e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0668382 < 1 %
Convergence criterion: norm(du)/norm(u) 0.162755 < 10 %
Convergence criterion: max(du) 7.38264e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 16
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.14431e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 81/100 time [yr]: 16.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 16.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.66356e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0664657 < 1 %
Convergence criterion: norm(du)/norm(u) 0.161561 < 10 %
Convergence criterion: max(du) 7.19865e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 16.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.64352e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 82/100 time [yr]: 16.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 16.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72958e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0660982 < 1 %
Convergence criterion: norm(du)/norm(u) 0.16078 < 10 %
Convergence criterion: max(du) 7.01901e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 16.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.07041e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 83/100 time [yr]: 16.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 16.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.64482e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0657357 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160396 < 10 %
Convergence criterion: max(du) 6.84347e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 16.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.58551e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 84/100 time [yr]: 16.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 16.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.75727e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0653783 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160366 < 10 %
Convergence criterion: max(du) 6.6718e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 16.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69339e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 85/100 time [yr]: 17.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 17
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.7468e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0650256 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160586 < 10 %
Convergence criterion: max(du) 6.5038e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 17
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.6756e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 86/100 time [yr]: 17.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 17.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.81505e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0646777 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160763 < 10 %
Convergence criterion: max(du) 6.34082e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 17.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.45107e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 87/100 time [yr]: 17.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 17.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71945e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0643345 < 1 %
Convergence criterion: norm(du)/norm(u) 0.161498 < 10 %
Convergence criterion: max(du) 6.18039e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 17.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.31165e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 88/100 time [yr]: 17.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 17.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.64909e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0639956 < 1 %
Convergence criterion: norm(du)/norm(u) 0.162772 < 10 %
Convergence criterion: max(du) 6.02325e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 17.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.71812e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 89/100 time [yr]: 17.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 17.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.64007e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0636613 < 1 %
Convergence criterion: norm(du)/norm(u) 0.164501 < 10 %
Convergence criterion: max(du) 5.86897e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 17.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.00274e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 90/100 time [yr]: 18.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 18
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.62871e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0633313 < 1 %
Convergence criterion: norm(du)/norm(u) 0.166412 < 10 %
Convergence criterion: max(du) 5.93924e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 18
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.63782e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 91/100 time [yr]: 18.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 18.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.65159e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.063006 < 1 %
Convergence criterion: norm(du)/norm(u) 0.168526 < 10 %
Convergence criterion: max(du) 6.18265e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 18.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.79145e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 92/100 time [yr]: 18.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 18.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.77171e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0626849 < 1 %
Convergence criterion: norm(du)/norm(u) 0.170745 < 10 %
Convergence criterion: max(du) 6.42703e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 18.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.74956e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 93/100 time [yr]: 18.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 18.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.62776e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0623676 < 1 %
Convergence criterion: norm(du)/norm(u) 0.173169 < 10 %
Convergence criterion: max(du) 6.67301e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 18.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.71158e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 94/100 time [yr]: 18.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 18.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.86864e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0620542 < 1 %
Convergence criterion: norm(du)/norm(u) 0.175781 < 10 %
Convergence criterion: max(du) 6.92124e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 18.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.64655e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 95/100 time [yr]: 19.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 19
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.54908e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0617442 < 1 %
Convergence criterion: norm(du)/norm(u) 0.178592 < 10 %
Convergence criterion: max(du) 7.17233e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 19
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.62561e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 96/100 time [yr]: 19.20 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 19.2
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.70731e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0614387 < 1 %
Convergence criterion: norm(du)/norm(u) 0.181552 < 10 %
Convergence criterion: max(du) 7.42693e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 19.2
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.9228e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 97/100 time [yr]: 19.40 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 19.4
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.6987e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.061137 < 1 %
Convergence criterion: norm(du)/norm(u) 0.184572 < 10 %
Convergence criterion: max(du) 7.68569e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 19.4
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.83829e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 98/100 time [yr]: 19.60 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 19.6
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.80897e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0608386 < 1 %
Convergence criterion: norm(du)/norm(u) 0.187808 < 10 %
Convergence criterion: max(du) 7.94928e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 19.6
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.62546e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 99/100 time [yr]: 19.80 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 19.8
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.73574e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0605439 < 1 %
Convergence criterion: norm(du)/norm(u) 0.19133 < 10 %
Convergence criterion: max(du) 8.21837e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 19.8
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.19089e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
iteration 100/100 time [yr]: 20.00 (time step: 0.20)
computing smb
computing new velocity
Updating constraints and active domain of analysis StressbalanceAnalysis for time: 20
Get solution from inputs
Reduce vector from g to f set
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
Allocating matrices (Kff stiffness matrix size: 9332 x 9332)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.66458e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0602527 < 1 %
Convergence criterion: norm(du)/norm(u) 0.194863 < 10 %
Convergence criterion: max(du) 8.49223e-07
Input updates from constant
Updating inputs from solution for StressbalanceAnalysis
number of unstable constraints: 0
total number of iterations: 1
computing basal mass balance
computing mass transport
call computational core
Updating constraints and active domain of analysis MasstransportAnalysis for time: 20
Allocating matrices (Kff stiffness matrix size: 4666 x 4666)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.51008e-16
Merging solution vector from fset to gset
Updating inputs from solution for MasstransportAnalysis
updating vertices positions
computing transient requested outputs
saving temporary results
write lock file:
FemModel initialization elapsed time: 0.026101
Total Core solution elapsed time: 13.9052
Linear solver elapsed time: 8.99336 (65%)
Total elapsed time: 0 hrs 0 min 13 sec
loading results from cluster
Step 5: Plotting
SUCCESS
finished: IceBridge
starting: IceflowModels
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
boundary conditions for stressbalance model: spc set as zero
no smb.mass_balance specified: values set as zero
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
Dot indexing is not supported for variables of this type.
Error in model/extrude (line 1276)
md.mmemasstransport=md.mmemasstransport.extrude(md);
Error in runme (line 6)
md=extrude(md,3,1);FAILURE
finished: IceflowModels
starting: Inversion
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
creating thickness
creating drag
creating flow law paramter
creating boundary conditions
boundary conditions for stressbalance model: spc set as zero
no smb.mass_balance specified: values set as zero
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
no thermal boundary conditions created: no observed temperature found
checking model consistency
model not consistent: field 'mmemasstransport' is not an object
Error using ismodelselfconsistent
Model not consistent, see messages above
Error in solve (line 100)
ismodelselfconsistent(md),
Error in runme (line 14)
md = solve(md,'Stressbalance');FAILURE
finished: Inversion
starting: Jakobshavn
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
Step 1: Mesh creation
Anisotropic mesh adaptation
WARNING: mesh present but no geometry found. Reconstructing...
new number of triangles = 3046
Step 2: Parameterization
Loading SeaRISE data from NetCDF
Interpolating thicknesses
Interpolating bedrock topography
Constructing surface elevation
Interpolating velocities
Interpolating temperatures
Interpolating surface mass balance
Construct basal friction parameters
Construct ice rheological properties
Set other boundary conditions
boundary conditions for stressbalance model: spc set as observed velocities
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
Step 3: Control method friction
checking model consistency
INFO: the outlog will look better if only md.verbose.control is turned on
marshalling file Jakobshavn.bin
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
preparing initial solution
x | Cost function f(x) | List of contributions
====================== step 1/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 89363.89 | 89280.58 83.31433
x = 1 | f(x) = 76234.74 | 76154.74 80.00512
====================== step 2/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 76416.09 | 76336.07 80.02122
x = 1 | f(x) = 66667.06 | 66590.22 76.84629
====================== step 3/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 66341.94 | 66265.12 76.82057
x = 1 | f(x) = 59236.12 | 59162.26 73.86055
====================== step 4/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 58920.94 | 58847.11 73.82904
x = 1 | f(x) = 53598.23 | 53527.14 71.08546
====================== step 5/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 53344.96 | 53273.91 71.05464
x = 1 | f(x) = 49643.94 | 49575.14 68.7982
====================== step 6/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 49456.08 | 49387.31 68.77171
x = 1 | f(x) = 46751.11 | 46684.22 66.88241
====================== step 7/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 46608.41 | 46541.55 66.85948
x = 1 | f(x) = 44738.76 | 44673.54 65.21319
====================== step 8/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 44634.79 | 44569.6 65.19433
x = 1 | f(x) = 43374.35 | 43310.55 63.80219
====================== step 9/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 43229.81 | 43166.04 63.77154
x = 1 | f(x) = 42186.1 | 42123.54 62.5627
====================== step 10/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 42059.62 | 41997.08 62.53554
x = 1 | f(x) = 41202.93 | 41141.47 61.45245
====================== step 11/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 41096.1 | 41034.68 61.42691
x = 1 | f(x) = 40360.09 | 40299.65 60.44179
====================== step 12/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 40267.04 | 40206.62 60.41795
x = 1 | f(x) = 39627.41 | 39567.9 59.50847
====================== step 13/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 39545.72 | 39486.23 59.48616
x = 1 | f(x) = 39003.04 | 38944.4 58.63879
====================== step 14/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 38933.02 | 38874.41 58.61774
x = 1 | f(x) = 38463.76 | 38405.94 57.81945
====================== step 15/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 38402.42 | 38344.62 57.79982
x = 1 | f(x) = 38037.59 | 37980.54 57.04745
x = 0.381966 | f(x) = 38149.66 | 38092.18 57.47378
x = 0.618034 | f(x) = 38099.72 | 38042.41 57.31471
x = 0.763932 | f(x) = 38051.17 | 37993.96 57.20521
x = 0.854102 | f(x) = 38012.55 | 37955.41 57.13404
x = 0.90983 | f(x) = 37983.95 | 37926.86 57.08845
====================== step 16/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 37971.42 | 37914.34 57.08297
x = 1 | f(x) = 37677.27 | 37620.88 56.3812
x = 0.381966 | f(x) = 37794.33 | 37737.55 56.78912
x = 0.618034 | f(x) = 37751.99 | 37695.36 56.63927
x = 0.763932 | f(x) = 37709.47 | 37652.93 56.53592
x = 0.854102 | f(x) = 37675.2 | 37618.73 56.46872
x = 0.90983 | f(x) = 37649.62 | 37593.2 56.42567
====================== step 17/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 37638.54 | 37582.12 56.42059
x = 1 | f(x) = 37369.28 | 37313.53 55.75399
x = 0.381966 | f(x) = 37475.4 | 37419.26 56.14166
x = 0.618034 | f(x) = 37437.25 | 37381.25 55.99942
x = 0.763932 | f(x) = 37398.42 | 37342.52 55.90123
x = 0.854102 | f(x) = 37366.93 | 37311.09 55.83739
x = 0.90983 | f(x) = 37343.35 | 37287.56 55.7965
====================== step 18/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 37333.02 | 37277.23 55.79169
x = 1 | f(x) = 37082.09 | 37026.94 55.15374
x = 0.381966 | f(x) = 37181.19 | 37125.66 55.5253
x = 0.618034 | f(x) = 37145.62 | 37090.23 55.38915
x = 0.763932 | f(x) = 37109.3 | 37054.01 55.29505
x = 0.854102 | f(x) = 37079.8 | 37024.57 55.23385
x = 0.90983 | f(x) = 37057.68 | 37002.49 55.19466
====================== step 19/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 37047.95 | 36992.76 55.19004
x = 1 | f(x) = 36822.28 | 36767.7 54.57745
x = 0.381966 | f(x) = 36906.9 | 36851.96 54.93407
x = 0.618034 | f(x) = 36874.3 | 36819.5 54.80343
x = 0.763932 | f(x) = 36842.53 | 36787.81 54.71314
x = 0.854102 | f(x) = 36816.93 | 36762.28 54.65441
x = 0.90983 | f(x) = 36797.79 | 36743.18 54.6168
====================== step 20/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 36789.31 | 36734.7 54.61236
x = 1 | f(x) = 36588.69 | 36534.67 54.02496
x = 0.381966 | f(x) = 36667.29 | 36612.93 54.36737
x = 0.618034 | f(x) = 36638.86 | 36584.62 54.24222
x = 0.763932 | f(x) = 36609.83 | 36555.67 54.15549
x = 0.854102 | f(x) = 36586.19 | 36532.09 54.09908
x = 0.90983 | f(x) = 36568.42 | 36514.36 54.06297
preparing final solution
computing new velocity
write lock file:
FemModel initialization elapsed time: 0.011652
Total Core solution elapsed time: 4.96431
Linear solver elapsed time: 4.08766 (82%)
Total elapsed time: 0 hrs 0 min 4 sec
loading results from cluster
Plotting
SUCCESS
finished: Jakobshavn
starting: LcurveAnalysis
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
creating thickness
creating drag
creating flow law paramter
creating boundary conditions
boundary conditions for stressbalance model: spc set as zero
no smb.mass_balance specified: values set as zero
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
no thermal boundary conditions created: no observed temperature found
checking model consistency
model not consistent: field 'mmemasstransport' is not an object
Error using ismodelselfconsistent
Model not consistent, see messages above
Error in solve (line 100)
ismodelselfconsistent(md),
Error in runme (line 14)
md = solve(md,'Stressbalance');FAILURE
finished: LcurveAnalysis
starting: Mesh
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
Construction of a mesh from a given geometry
Construction of a mesh from a given geometry
Construction of a mesh from a given geometry
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 2587
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 1622
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 4550
SUCCESS
finished: Mesh
starting: Pig
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 2075
Loading SeaRISE data from NetCDF
Interpolating surface and ice base
Constructing thickness
Interpolating temperatures
Loading velocities data from NetCDF
Set observed velocities
Set Pressure
Construct ice rheological properties
Interpolating surface mass balance
Set geothermal heat flux
Construct basal friction parameters
Set boundary conditions
boundary conditions for stressbalance model: spc set as observed velocities
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
Initialize M1QN3 parameters
Computing initial solution
┌────┬─────────────────┬────────────┬──────────────────────────────┐
│Iter│ Cost function │ Grad. norm │ List of contributions │
├────┼─────────────────┼────────────┼──────────────────────────────┤
│ 1 │ f(x)= 10553 │ 49.9 │ 1481 9071 3.153e-11│
│ 2 │ f(x)= 3074.6 │ 6.54 │ 59.58 3015 5.056e-10│
│ 3 │ f(x)= 2636.1 │ 5.27 │ 45.38 2591 6.23e-10│
│ 4 │ f(x)= 1451.8 │ 2.43 │ 17.34 1434 1.459e-09│
│ 5 │ f(x)= 913.05 │ 1.36 │ 7.811 905.2 2.589e-09│
│ 6 │ f(x)= 519.54 │ 0.681 │ 3.009 516.5 4.766e-09│
│ 7 │ f(x)= 300.11 │ 0.321 │ 1.582 298.5 8.128e-09│
│ 8 │ f(x)= 187.72 │ 0.181 │ 1.306 186.4 1.252e-08│
│ 9 │ f(x)= 185.22 │ 0.843 │ 1.484 183.7 1.736e-08│
│ 10 │ f(x)= 142.19 │ 0.164 │ 1.455 140.7 1.785e-08│
│ 11 │ f(x)= 137.81 │ 0.137 │ 1.484 136.3 1.83e-08│
│ 12 │ f(x)= 125.8 │ 0.0796 │ 1.535 124.3 2.039e-08│
│ 13 │ f(x)= 121.31 │ 0.0608 │ 1.518 119.8 2.132e-08│
│ 14 │ f(x)= 114.45 │ 0.057 │ 1.44 113 2.267e-08│
│ 15 │ f(x)= 110.88 │ 0.0557 │ 1.419 109.5 2.36e-08│
│ 16 │ f(x)= 108.01 │ 0.0509 │ 1.389 106.6 2.419e-08│
│ 17 │ f(x)= 105.92 │ 0.0371 │ 1.368 104.6 2.483e-08│
│ 18 │ f(x)= 104.32 │ 0.0319 │ 1.355 103 2.562e-08│
│ 19 │ f(x)= 102.04 │ 0.0404 │ 1.342 100.7 2.802e-08│
│ 20 │ f(x)= 101.22 │ 0.0607 │ 1.321 99.9 2.957e-08│
│ 21 │ f(x)= 100.39 │ 0.0737 │ 1.304 99.09 2.924e-08│
└────┴─────────────────┴────────────┴──────────────────────────────┘
Maximum number of iterations exceeded
preparing final solution
write lock file:
FemModel initialization elapsed time: 0.012924
Total Core solution elapsed time: 9.18539
Linear solver elapsed time: 7.13816 (78%)
Total elapsed time: 0 hrs 0 min 9 sec
Extruding mesh
Using HO Ice Flow Model
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
write lock file:
FemModel initialization elapsed time: 0.04758
Total Core solution elapsed time: 1.57551
Linear solver elapsed time: 1.12583 (71%)
Total elapsed time: 0 hrs 0 min 1 sec
SUCCESS
finished: Pig
starting: PigSensitivity
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
write lock file:
FemModel initialization elapsed time: 0.019999
Total Core solution elapsed time: 16.4712
Linear solver elapsed time: 11.7367 (71%)
Total elapsed time: 0 hrs 0 min 16 sec
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
write lock file:
FemModel initialization elapsed time: 0.019232
Total Core solution elapsed time: 6.72305
Linear solver elapsed time: 4.39759 (65%)
Total elapsed time: 0 hrs 0 min 6 sec
NOTE: using observed velocities to create constraints along new boundary
boundary conditions for stressbalance model: spc set as observed velocities
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
write lock file:
FemModel initialization elapsed time: 0.011799
Total Core solution elapsed time: 6.47448
Linear solver elapsed time: 4.29222 (66%)
Total elapsed time: 0 hrs 0 min 6 sec
uploading input file and queuing script
launching solution sequence on remote cluster
Ice-sheet and Sea-level System Model (ISSM) version 4.24
(website: http://issm.jpl.nasa.gov forum: https://issm.ess.uci.edu/forum/)
call computational core:
write lock file:
FemModel initialization elapsed time: 0.026176
Total Core solution elapsed time: 5.13369
Linear solver elapsed time: 3.45816 (67%)
Total elapsed time: 0 hrs 0 min 5 sec
SUCCESS
finished: PigSensitivity
starting: SquareIceShelf
< M A T L A B (R) >
Copyright 1984-2023 The MathWorks, Inc.
R2023b Update 6 (23.2.0.2485118) 64-bit (maca64)
December 28, 2023
To get started, type doc.
For product information, visit www.mathworks.com.
ISSM development path correctly loaded
creating thickness
creating drag
initial velocity
creating flow law parameter
creating boundary conditions
boundary conditions for stressbalance model: spc set as zero
no smb.mass_balance specified: values set as zero
no basalforcings.groundedice_melting_rate specified: values set as zero
no basalforcings.floatingice_melting_rate specified: values set as zero
no balancethickness.thickening_rate specified: values set as zero
no thermal boundary conditions created: no observed temperature found
checking model consistency
model not consistent: field 'mmemasstransport' is not an object
Error using ismodelselfconsistent
Model not consistent, see messages above
Error in solve (line 100)
ismodelselfconsistent(md),
Error in runme (line 9)
md=solve(md,'Stressbalance');FAILURE
finished: SquareIceShelf
-----------End of matlab_log.log-----------
Build step 'Execute shell' marked build as failure
Recording test results
Finished: FAILURE
![[Jenkins]](/jenkins/static/e550d957/images/svgs/logo.svg){kind=logo}
