mathieumorlighem
I double checked my case using as default solver and it seems that thermal model wasn't the one that blew up.
Would you like to have the ready-to-run model file?
call computational core:
iteration 1/10000 time [yr]: 23.00 (time step: 0.00)
updating boundary conditions...
computing enthalpy
solver residue: norm(KU-F)/norm(F)=1.17029e-07
mechanical equilibrium convergence criterion 1.17029e-05 < 5 %
Convergence criterion: norm(du)/norm(u) 0 < 1 %
Convergence criterion: max(du) 0
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)=4.34517e-06
mechanical equilibrium convergence criterion 0.191162 > 0.1 %
Convergence criterion: norm(du)/norm(u) 0.0194689 < 1 %
Convergence criterion: max(du) 1.79719e-07
number of unstable constraints: 0
solver residue: norm(KU-F)/norm(F)=5.7609e-06
mechanical equilibrium convergence criterion 0.0773657 < 0.1 %
Convergence criterion: norm(du)/norm(u) 0.00386815 < 1 %
Convergence criterion: max(du) 3.75391e-08
number of unstable constraints: 0
total number of iterations: 2
computing basal mass balance
computing vertical velocities
solver residue: norm(KU-F)/norm(F)=-nan
[4] ??? Error using ==> ./modules/Solverx/Solverx.cpp:63
[4] checkconvergence error message: Solver residue is NaN
[10] ??? Error using ==> ./modules/Solverx/Solverx.cpp:63
[10] checkconvergence error message: Solver residue is NaN
[8] ??? Error using ==> ./modules/Solverx/Solverx.cpp:63
[8] checkconvergence error message: Solver residue is NaN