I am trying to introduce basal sliding into a thermo-mechanically coupled model for a polythermal mountain glacier. I use snow/firn rheology (Porous Solver). If I apply no-slip conditions at the bedrock, the case works. But when I set the basal boundary condition as described in Elmer/Ice Wiki (http://elmerfem.org/elmerice/wiki/doku. ... ers:porous), no convergence is seen even after 10 000 linear system solver iterations. The error message is:
Code: Select all
9990 0.9374E+00
10000 0.9391E+00
10001 0.9391E+00
NUMERICAL ERROR:: IterSolve: Too many iterations was needed.
Note: The following floating-point exceptions are signalling: IEEE_DENORMAL
Regards,
Gleb