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:
The sif and mesh files are attached. Also, it may worth mentioning that the solver is restarted from the previous solution (unfortunately, the restart file is too large to be attached). I will also provide additional information, if needed. I do not know what is incorrect and how to obtain convergence. So, I would be very grateful for any help or advice.
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