mathieumorlighem
Indeed, turning off md.timestepping.cycle_forcing
solves the problem. Although I am also puzzled why it was turned on when loadmodel
converted the mat
model file.
I would like to also ask one question I've had for a while. What does "Linear solver elapsed time" mean, and does that 34% seem all right for this HO case with iterative solver?
Thanks,
Wade
iteration 17/17 time [yr]: 1.00 (time step: 0.04)
updating boundary conditions...
computing enthalpy
** PETSc DEPRECATION WARNING ** : the option -pc_factor_mat_solver_package is deprecated as of version 3.9 and will be removed in a future release. Please use the option -pc_factor_mat_solver_type instead. (Silence this warning with -options_suppress_deprecated_warnings)
solver residue: norm(KU-F)/norm(F)=3.51206e-16
mechanical equilibrium convergence criterion 2.27501e-05 < 5 %
Convergence criterion: norm(du)/norm(u) 0.00373949 < 1 %
Convergence criterion: max(du) 514.134
number of unstable constraints: 0
total number of iterations: 1
depth averaging WaterfractionDrainage
extruding WaterfractionDrainageIntegrated from base...
extruding BasalforcingsGroundediceMeltingRate from base...
computing smb
call Delta18opdParameterization module
call positive degree day module
computing new velocity
solver residue: norm(KU-F)/norm(F)=5.0072e-06
mechanical equilibrium convergence criterion 0.151226 > 0.1 %
Convergence criterion: norm(du)/norm(u) 0.0177775 < 1 %
Convergence criterion: max(du) 5.73974e-08
number of unstable constraints: 0
solver residue: norm(KU-F)/norm(F)=5.196e-06
mechanical equilibrium convergence criterion 0.0640812 < 0.1 %
Convergence criterion: norm(du)/norm(u) 0.00444832 < 1 %
Convergence criterion: max(du) 3.52178e-08
number of unstable constraints: 0
total number of iterations: 2
computing basal mass balance
computing vertical velocities
** PETSc DEPRECATION WARNING ** : the option -pc_factor_mat_solver_package is deprecated as of version 3.9 and will be removed in a future release. Please use the option -pc_factor_mat_solver_type instead. (Silence this warning with -options_suppress_deprecated_warnings)
solver residue: norm(KU-F)/norm(F)=3.69218e-15
computing basal mass balance
computing mass transport
call computational core
depth averaging Vx
depth averaging Vy
** PETSc DEPRECATION WARNING ** : the option -pc_factor_mat_solver_package is deprecated as of version 3.9 and will be removed in a future release. Please use the option -pc_factor_mat_solver_type instead. (Silence this warning with -options_suppress_deprecated_warnings)
solver residue: norm(KU-F)/norm(F)=1.48422e-16
extruding Thickness from base...
extruding Base from base...
extruding Surface from base...
updating vertices positions
computing transient requested outputs
saving temporary results
write lock file:
FemModel initialization elapsed time: 6.47204
Total Core solution elapsed time: 244.825
Linear solver elapsed time: 84.3585 (34%)
Total elapsed time: 0 hrs 4 min 11 sec