Convergence criterion: norm(du)/norm(u) 1.07113 < 10 %
Convergence criterion: max(du) 6.24711e-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: 3.8
Allocating matrices (Kff stiffness matrix size: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.66645e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84265e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.106356 < 1 %
Convergence criterion: norm(du)/norm(u) 0.998493 < 10 %
Convergence criterion: max(du) 5.75025e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.72593e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.86958e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.104206 < 1 %
Convergence criterion: norm(du)/norm(u) 0.937125 < 10 %
Convergence criterion: max(du) 5.3407e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81237e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.98568e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.102251 < 1 %
Convergence criterion: norm(du)/norm(u) 0.884247 < 10 %
Convergence criterion: max(du) 4.98966e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.66645e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.1083e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.100465 < 1 %
Convergence criterion: norm(du)/norm(u) 0.838226 < 10 %
Convergence criterion: max(du) 4.68095e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.89015e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.78994e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0988228 < 1 %
Convergence criterion: norm(du)/norm(u) 0.797679 < 10 %
Convergence criterion: max(du) 4.40552e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.64083e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.92886e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0973049 < 1 %
Convergence criterion: norm(du)/norm(u) 0.761639 < 10 %
Convergence criterion: max(du) 4.15807e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.6161e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.92949e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.09589 < 1 %
Convergence criterion: norm(du)/norm(u) 0.729262 < 10 %
Convergence criterion: max(du) 3.93508e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.57808e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.83235e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0945622 < 1 %
Convergence criterion: norm(du)/norm(u) 0.699987 < 10 %
Convergence criterion: max(du) 3.73383e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.78475e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.98087e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0933088 < 1 %
Convergence criterion: norm(du)/norm(u) 0.673373 < 10 %
Convergence criterion: max(du) 3.55194e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.76339e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.96118e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0921186 < 1 %
Convergence criterion: norm(du)/norm(u) 0.649074 < 10 %
Convergence criterion: max(du) 3.38733e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.01082e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.47349e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0909765 < 1 %
Convergence criterion: norm(du)/norm(u) 0.626782 < 10 %
Convergence criterion: max(du) 3.23787e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.90637e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.8339e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0898938 < 1 %
Convergence criterion: norm(du)/norm(u) 0.606266 < 10 %
Convergence criterion: max(du) 3.1018e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.91639e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.00393e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0888679 < 1 %
Convergence criterion: norm(du)/norm(u) 0.587377 < 10 %
Convergence criterion: max(du) 2.97749e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.87348e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.82308e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0878845 < 1 %
Convergence criterion: norm(du)/norm(u) 0.570048 < 10 %
Convergence criterion: max(du) 2.86347e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.43458e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.76412e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0869359 < 1 %
Convergence criterion: norm(du)/norm(u) 0.554064 < 10 %
Convergence criterion: max(du) 2.75845e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.78768e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84021e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0860176 < 1 %
Convergence criterion: norm(du)/norm(u) 0.539173 < 10 %
Convergence criterion: max(du) 2.6613e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.44756e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.61044e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0851273 < 1 %
Convergence criterion: norm(du)/norm(u) 0.525203 < 10 %
Convergence criterion: max(du) 2.5711e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.85222e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.8854e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0842647 < 1 %
Convergence criterion: norm(du)/norm(u) 0.51205 < 10 %
Convergence criterion: max(du) 2.48706e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.50559e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.80201e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0834299 < 1 %
Convergence criterion: norm(du)/norm(u) 0.499628 < 10 %
Convergence criterion: max(du) 2.40854e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.74048e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.10949e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0826229 < 1 %
Convergence criterion: norm(du)/norm(u) 0.487841 < 10 %
Convergence criterion: max(du) 2.33498e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.66361e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.89524e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.081843 < 1 %
Convergence criterion: norm(du)/norm(u) 0.47662 < 10 %
Convergence criterion: max(du) 2.2659e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.91709e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.6372e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0810856 < 1 %
Convergence criterion: norm(du)/norm(u) 0.465877 < 10 %
Convergence criterion: max(du) 2.20091e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.68362e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.2144e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0803506 < 1 %
Convergence criterion: norm(du)/norm(u) 0.455531 < 10 %
Convergence criterion: max(du) 2.13964e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.76118e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.12521e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0796359 < 1 %
Convergence criterion: norm(du)/norm(u) 0.445517 < 10 %
Convergence criterion: max(du) 2.08176e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.91479e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.93277e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.07894 < 1 %
Convergence criterion: norm(du)/norm(u) 0.435768 < 10 %
Convergence criterion: max(du) 2.02701e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.75359e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.11492e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0782614 < 1 %
Convergence criterion: norm(du)/norm(u) 0.42629 < 10 %
Convergence criterion: max(du) 1.97513e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96235e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.0283e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0775988 < 1 %
Convergence criterion: norm(du)/norm(u) 0.417072 < 10 %
Convergence criterion: max(du) 1.9259e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.65201e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.78005e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0769509 < 1 %
Convergence criterion: norm(du)/norm(u) 0.408118 < 10 %
Convergence criterion: max(du) 1.8791e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96853e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.87199e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0763168 < 1 %
Convergence criterion: norm(du)/norm(u) 0.399417 < 10 %
Convergence criterion: max(du) 1.83457e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.90218e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.85079e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0756956 < 1 %
Convergence criterion: norm(du)/norm(u) 0.390975 < 10 %
Convergence criterion: max(du) 1.79214e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.79429e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.55963e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0750865 < 1 %
Convergence criterion: norm(du)/norm(u) 0.382779 < 10 %
Convergence criterion: max(du) 1.75173e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.84098e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.91218e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0744891 < 1 %
Convergence criterion: norm(du)/norm(u) 0.374844 < 10 %
Convergence criterion: max(du) 1.71313e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.86192e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.97845e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0739079 < 1 %
Convergence criterion: norm(du)/norm(u) 0.362827 < 10 %
Convergence criterion: max(du) 1.6621e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.84787e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.02521e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0733373 < 1 %
Convergence criterion: norm(du)/norm(u) 0.350628 < 10 %
Convergence criterion: max(du) 1.60939e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.6243e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.26051e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0727769 < 1 %
Convergence criterion: norm(du)/norm(u) 0.33879 < 10 %
Convergence criterion: max(du) 1.55893e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.7575e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.85945e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0722264 < 1 %
Convergence criterion: norm(du)/norm(u) 0.327423 < 10 %
Convergence criterion: max(du) 1.51167e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.52899e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84472e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0716855 < 1 %
Convergence criterion: norm(du)/norm(u) 0.316536 < 10 %
Convergence criterion: max(du) 1.4676e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.89639e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.10068e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0711543 < 1 %
Convergence criterion: norm(du)/norm(u) 0.306126 < 10 %
Convergence criterion: max(du) 1.42643e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81594e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.16957e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0706323 < 1 %
Convergence criterion: norm(du)/norm(u) 0.296179 < 10 %
Convergence criterion: max(du) 1.38775e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.74623e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.85091e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0701194 < 1 %
Convergence criterion: norm(du)/norm(u) 0.286654 < 10 %
Convergence criterion: max(du) 1.35119e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96029e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.16526e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0696153 < 1 %
Convergence criterion: norm(du)/norm(u) 0.277549 < 10 %
Convergence criterion: max(du) 1.3164e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.67499e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.99476e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0691199 < 1 %
Convergence criterion: norm(du)/norm(u) 0.268864 < 10 %
Convergence criterion: max(du) 1.28311e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.36332e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.77345e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.068633 < 1 %
Convergence criterion: norm(du)/norm(u) 0.260586 < 10 %
Convergence criterion: max(du) 1.25112e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.70329e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.89116e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0681543 < 1 %
Convergence criterion: norm(du)/norm(u) 0.252703 < 10 %
Convergence criterion: max(du) 1.22026e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.4485e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.75858e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0676835 < 1 %
Convergence criterion: norm(du)/norm(u) 0.245188 < 10 %
Convergence criterion: max(du) 1.19043e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.20605e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.45013e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0672203 < 1 %
Convergence criterion: norm(du)/norm(u) 0.237927 < 10 %
Convergence criterion: max(du) 1.16154e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.5413e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.0776e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0667652 < 1 %
Convergence criterion: norm(du)/norm(u) 0.231092 < 10 %
Convergence criterion: max(du) 1.13352e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96574e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84215e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0663172 < 1 %
Convergence criterion: norm(du)/norm(u) 0.22462 < 10 %
Convergence criterion: max(du) 1.10633e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.89178e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.07247e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0658764 < 1 %
Convergence criterion: norm(du)/norm(u) 0.218499 < 10 %
Convergence criterion: max(du) 1.0799e-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.6
Allocating matrices (Kff stiffness matrix size: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.99862e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.04564e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0654425 < 1 %
Convergence criterion: norm(du)/norm(u) 0.212706 < 10 %
Convergence criterion: max(du) 1.0542e-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.8
Allocating matrices (Kff stiffness matrix size: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.61593e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.23745e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0650153 < 1 %
Convergence criterion: norm(du)/norm(u) 0.20725 < 10 %
Convergence criterion: max(du) 1.02919e-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: 14
Allocating matrices (Kff stiffness matrix size: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81028e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71923e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0645949 < 1 %
Convergence criterion: norm(du)/norm(u) 0.202082 < 10 %
Convergence criterion: max(du) 1.00482e-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: 14.2
Allocating matrices (Kff stiffness matrix size: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.96433e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.19009e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0641809 < 1 %
Convergence criterion: norm(du)/norm(u) 0.197233 < 10 %
Convergence criterion: max(du) 9.81074e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.8346e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.16581e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0637733 < 1 %
Convergence criterion: norm(du)/norm(u) 0.192578 < 10 %
Convergence criterion: max(du) 9.57908e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.62684e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.67504e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.063372 < 1 %
Convergence criterion: norm(du)/norm(u) 0.188381 < 10 %
Convergence criterion: max(du) 9.3529e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.92677e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.89323e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0629767 < 1 %
Convergence criterion: norm(du)/norm(u) 0.18449 < 10 %
Convergence criterion: max(du) 9.13182e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.5288e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.21461e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0625876 < 1 %
Convergence criterion: norm(du)/norm(u) 0.180899 < 10 %
Convergence criterion: max(du) 8.91577e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.59969e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.13388e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0622044 < 1 %
Convergence criterion: norm(du)/norm(u) 0.177604 < 10 %
Convergence criterion: max(du) 8.70459e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.50977e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.72801e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0618266 < 1 %
Convergence criterion: norm(du)/norm(u) 0.174598 < 10 %
Convergence criterion: max(du) 8.49801e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.82491e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.89715e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0614541 < 1 %
Convergence criterion: norm(du)/norm(u) 0.171868 < 10 %
Convergence criterion: max(du) 8.42481e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.93443e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.81123e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0610872 < 1 %
Convergence criterion: norm(du)/norm(u) 0.169371 < 10 %
Convergence criterion: max(du) 8.5016e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.7539e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.88578e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0607257 < 1 %
Convergence criterion: norm(du)/norm(u) 0.166878 < 10 %
Convergence criterion: max(du) 8.56847e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.78423e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.85961e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0603694 < 1 %
Convergence criterion: norm(du)/norm(u) 0.164926 < 10 %
Convergence criterion: max(du) 8.62499e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.59393e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.64344e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0600192 < 1 %
Convergence criterion: norm(du)/norm(u) 0.163357 < 10 %
Convergence criterion: max(du) 8.72759e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81506e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.80604e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0596741 < 1 %
Convergence criterion: norm(du)/norm(u) 0.162021 < 10 %
Convergence criterion: max(du) 8.81599e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.43831e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.71441e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0593326 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160712 < 10 %
Convergence criterion: max(du) 8.84127e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.01431e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.78188e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0589956 < 1 %
Convergence criterion: norm(du)/norm(u) 0.159568 < 10 %
Convergence criterion: max(du) 8.8477e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.54324e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.6623e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0586632 < 1 %
Convergence criterion: norm(du)/norm(u) 0.158695 < 10 %
Convergence criterion: max(du) 8.83843e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.69777e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.75149e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0583353 < 1 %
Convergence criterion: norm(du)/norm(u) 0.158032 < 10 %
Convergence criterion: max(du) 8.81648e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.50955e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.99126e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0580117 < 1 %
Convergence criterion: norm(du)/norm(u) 0.157615 < 10 %
Convergence criterion: max(du) 8.78242e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.92843e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.67196e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0576926 < 1 %
Convergence criterion: norm(du)/norm(u) 0.157574 < 10 %
Convergence criterion: max(du) 8.73655e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.88558e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84008e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0573778 < 1 %
Convergence criterion: norm(du)/norm(u) 0.157866 < 10 %
Convergence criterion: max(du) 8.67918e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.83062e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.55012e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0570669 < 1 %
Convergence criterion: norm(du)/norm(u) 0.158338 < 10 %
Convergence criterion: max(du) 8.61067e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.7026e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.89936e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0567599 < 1 %
Convergence criterion: norm(du)/norm(u) 0.158938 < 10 %
Convergence criterion: max(du) 8.53142e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.43009e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.77822e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0564556 < 1 %
Convergence criterion: norm(du)/norm(u) 0.159689 < 10 %
Convergence criterion: max(du) 8.44192e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81336e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.01067e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0561557 < 1 %
Convergence criterion: norm(du)/norm(u) 0.160629 < 10 %
Convergence criterion: max(du) 8.34268e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.70326e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.84373e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.05586 < 1 %
Convergence criterion: norm(du)/norm(u) 0.161758 < 10 %
Convergence criterion: max(du) 8.23427e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.82575e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=3.18165e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0555678 < 1 %
Convergence criterion: norm(du)/norm(u) 0.162987 < 10 %
Convergence criterion: max(du) 8.11726e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.81781e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.81048e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0552798 < 1 %
Convergence criterion: norm(du)/norm(u) 0.164309 < 10 %
Convergence criterion: max(du) 7.9923e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.67573e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.93209e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0549959 < 1 %
Convergence criterion: norm(du)/norm(u) 0.166002 < 10 %
Convergence criterion: max(du) 7.86001e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.15226e-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: 9300 x 9300)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=2.90962e-16
Merging solution vector from fset to gset
checking convergence
mechanical equilibrium convergence criterion 0.0547161 < 1 %
Convergence criterion: norm(du)/norm(u) 0.167956 < 10 %
Convergence criterion: max(du) 7.74989e-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: 4650 x 4650)
Assembling matrices
Create nodal constraints
Dirichlet lifting applied to load vector
Solving matrix system
solver residue: norm(KU-F)/norm(F)=1.87871e-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.116679
Total Core solution elapsed time: 29.9845
Linear solver elapsed time: 15.3461 (51%)
Total elapsed time: 0 hrs 0 min 30 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 (maci64)
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 (maci64)
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 (maci64)
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 = 3007
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) = 83691.82 | 83608.87 82.9538
x = 1 | f(x) = 71282.22 | 71202.56 79.65915
====================== step 2/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 71434.03 | 71354.35 79.67283
x = 1 | f(x) = 61944.77 | 61868.32 76.4542
====================== step 3/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 61661.08 | 61584.65 76.42878
x = 1 | f(x) = 55044.64 | 54971.17 73.46996
====================== step 4/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 54783.23 | 54709.79 73.44058
x = 1 | f(x) = 49978.08 | 49907.35 70.7277
====================== step 5/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 49772.38 | 49701.68 70.70017
x = 1 | f(x) = 46134.88 | 46066.67 68.21092
====================== step 6/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 45971.44 | 45903.26 68.18565
x = 1 | f(x) = 43246.98 | 43181.03 65.95402
====================== step 7/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 43116.42 | 43050.49 65.93162
x = 1 | f(x) = 41186.24 | 41122.2 64.04128
====================== step 8/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 41086.83 | 41022.81 64.02207
x = 1 | f(x) = 39691.01 | 39628.64 62.36093
====================== step 9/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 39617.72 | 39555.37 62.34435
x = 1 | f(x) = 38647.74 | 38586.78 60.96179
====================== step 10/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 38542.7 | 38481.76 60.93448
x = 1 | f(x) = 37742.73 | 37683.01 59.72204
====================== step 11/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 37651.01 | 37591.32 59.69631
x = 1 | f(x) = 36980.6 | 36922 58.60387
====================== step 12/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 36901.71 | 36843.13 58.58021
x = 1 | f(x) = 36328.21 | 36270.64 57.57771
====================== step 13/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 36259.35 | 36201.79 57.55583
x = 1 | f(x) = 35760.77 | 35704.15 56.6245
====================== step 14/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 35699.86 | 35643.25 56.60405
x = 1 | f(x) = 35257.21 | 35201.48 55.7326
====================== step 15/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 35202.45 | 35146.74 55.71261
x = 1 | f(x) = 34832.95 | 34778.07 54.88295
====================== step 16/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 34786.61 | 34731.75 54.86544
x = 1 | f(x) = 34479.86 | 34425.78 54.08727
x = 0.381966 | f(x) = 34579.35 | 34524.82 54.53269
x = 0.618034 | f(x) = 34536.62 | 34482.25 54.36682
x = 0.763932 | f(x) = 34495.21 | 34440.96 54.2533
x = 0.854102 | f(x) = 34462.49 | 34408.31 54.17991
x = 0.90983 | f(x) = 34438.41 | 34384.27 54.1331
====================== step 17/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 34428.06 | 34373.94 54.12786
x = 1 | f(x) = 34186.62 | 34133.22 53.39879
x = 0.381966 | f(x) = 34281.41 | 34227.59 53.82491
x = 0.618034 | f(x) = 34245.75 | 34192.08 53.66858
x = 0.763932 | f(x) = 34211.4 | 34157.84 53.56115
x = 0.854102 | f(x) = 34184.23 | 34130.74 53.49158
x = 0.90983 | f(x) = 34164.21 | 34110.77 53.4472
====================== step 18/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 34155.78 | 34102.33 53.44227
x = 1 | f(x) = 33946.97 | 33894.23 52.74752
x = 0.381966 | f(x) = 34030.65 | 33977.49 53.1537
x = 0.618034 | f(x) = 34000.42 | 33947.42 53.00515
x = 0.763932 | f(x) = 33970.44 | 33917.53 52.90262
x = 0.854102 | f(x) = 33946.39 | 33893.55 52.83617
x = 0.90983 | f(x) = 33928.5 | 33875.7 52.79378
====================== step 19/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 33920.82 | 33868.03 52.789
x = 1 | f(x) = 33740.96 | 33688.83 52.12422
x = 0.381966 | f(x) = 33811.12 | 33758.61 52.51346
x = 0.618034 | f(x) = 33784.82 | 33732.45 52.37163
x = 0.763932 | f(x) = 33759.33 | 33707.05 52.27353
x = 0.854102 | f(x) = 33739.07 | 33686.86 52.20992
x = 0.90983 | f(x) = 33724.1 | 33671.93 52.16935
====================== step 20/20 ===============================
x = 0 | computing velocities
computing adjoint
saving results
f(x) = 33717.79 | 33665.62 52.16487
x = 1 | f(x) = 33583.19 | 33531.66 51.52454
x = 0.381966 | f(x) = 33635.61 | 33583.71 51.90026
x = 0.618034 | f(x) = 33615.74 | 33563.98 51.76344
x = 0.763932 | f(x) = 33597 | 33545.33 51.66875
x = 0.854102 | f(x) = 33582.39 | 33530.78 51.60748
x = 0.90983 | f(x) = 33571.73 | 33520.16 51.56847
preparing final solution
computing new velocity
write lock file:
FemModel initialization elapsed time: 0.075516
Total Core solution elapsed time: 11.8754
Linear solver elapsed time: 8.1477 (69%)
Total elapsed time: 0 hrs 0 min 11 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 (maci64)
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 (maci64)
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 = 2567
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 1572
Construction of a mesh from a given geometry
Anisotropic mesh adaptation
new number of triangles = 4616
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 (maci64)
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 = 2085
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)= 10690 │ 51.3 │ 1475 9214 3.271e-11│
│ 2 │ f(x)= 3165.2 │ 6.8 │ 68.4 3097 5.082e-10│
│ 3 │ f(x)= 2709.5 │ 5.43 │ 51.99 2658 6.208e-10│
│ 4 │ f(x)= 1527 │ 2.59 │ 19.93 1507 1.367e-09│
│ 5 │ f(x)= 948.76 │ 1.43 │ 8.873 939.9 2.507e-09│
│ 6 │ f(x)= 538.96 │ 0.704 │ 3.362 535.6 4.587e-09│
│ 7 │ f(x)= 320.85 │ 0.352 │ 1.654 319.2 7.502e-09│
│ 8 │ f(x)= 198.61 │ 0.154 │ 1.284 197.3 1.178e-08│
│ 9 │ f(x)= 157.58 │ 0.319 │ 1.336 156.2 1.612e-08│
│ 10 │ f(x)= 134.5 │ 0.116 │ 1.478 133 1.929e-08│
│ 11 │ f(x)= 129.99 │ 0.0779 │ 1.476 128.5 1.983e-08│
│ 12 │ f(x)= 122.97 │ 0.0862 │ 1.467 121.5 2.173e-08│
│ 13 │ f(x)= 119.93 │ 0.0654 │ 1.461 118.5 2.225e-08│
│ 14 │ f(x)= 114.49 │ 0.0559 │ 1.408 113.1 2.332e-08│
│ 15 │ f(x)= 111.94 │ 0.0419 │ 1.394 110.5 2.418e-08│
│ 16 │ f(x)= 110.08 │ 0.0324 │ 1.399 108.7 2.511e-08│
│ 17 │ f(x)= 106.57 │ 0.0347 │ 1.392 105.2 2.754e-08│
│ 18 │ f(x)= 104.78 │ 0.0645 │ 1.351 103.4 3.038e-08│
│ 19 │ f(x)= 103.26 │ 0.0712 │ 1.319 101.9 3.056e-08│
│ 20 │ f(x)= 103.49 │ 0.0697 │ 1.293 102.2 3.191e-08│
│ 21 │ f(x)= 102.85 │ 0.0361 │ 1.296 101.6 3.106e-08│
│ 22 │ f(x)= 102.51 │ 0.0854 │ 1.293 101.2 3.086e-08│
└────┴─────────────────┴────────────┴──────────────────────────────┘
Maximum number of iterations exceeded
preparing final solution
write lock file:
FemModel initialization elapsed time: 0.081722
Total Core solution elapsed time: 5.39931
Linear solver elapsed time: 3.36486 (62%)
Total elapsed time: 0 hrs 0 min 5 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.09567
Total Core solution elapsed time: 1.77715
Linear solver elapsed time: 1.07022 (60%)
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 (maci64)
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.056993
Total Core solution elapsed time: 8.27012
Linear solver elapsed time: 4.80078 (58%)
Total elapsed time: 0 hrs 0 min 8 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.05695
Total Core solution elapsed time: 8.7395
Linear solver elapsed time: 5.01 (57%)
Total elapsed time: 0 hrs 0 min 8 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.054609
Total Core solution elapsed time: 8.39182
Linear solver elapsed time: 4.86908 (58%)
Total elapsed time: 0 hrs 0 min 8 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.055077
Total Core solution elapsed time: 8.27935
Linear solver elapsed time: 4.81606 (58%)
Total elapsed time: 0 hrs 0 min 8 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 (maci64)
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}
