----------------------------------------------------------------------------- -- This is ASPECT -- -- The Advanced Solver for Planetary Evolution, Convection, and Tectonics. -- ----------------------------------------------------------------------------- -- . version 2.6.0-pre (main, e3a102e7b) -- . using deal.II 9.5.0 -- . with 32 bit indices and vectorization level 2 (256 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in DEBUG mode -- . running with 4 MPI processes ----------------------------------------------------------------------------- Loading Ascii data boundary file /home/aspect/aspect/build/test-gravity/simple-geometry.txt. ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&NewtonSolver=1&sha=e3a102e7b&src=code ----------------------------------------------------------------------------- Number of active cells: 800 (on 3 levels) Number of degrees of freedom: 20,787 (6,642+861+3,321+3,321+3,321+3,321) Number of mesh deformation degrees of freedom: 1722 Solving mesh displacement system... 0 iterations. Added initial topography to grid *** Timestep 0: t=0 years, dt=0 years Solving mesh displacement system... 0 iterations. Solving temperature system... 0 iterations. Solving crust_upper system ... 0 iterations. Solving crust_lower system ... 0 iterations. Solving mantle_lithosphere system ... 0 iterations. Initial Newton Stokes residual = 1.3689e+12, v = 1.3689e+12, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 95+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 1, norm of the rhs: 1.3689e+12 Rebuilding Stokes preconditioner... Solving Stokes system... 146+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 20.0903, norm of the rhs: 2.75017e+13, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+3 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 17.5963, norm of the rhs: 2.40877e+13, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+15 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 13.5837, norm of the rhs: 1.85947e+13, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 8.96176, norm of the rhs: 1.22678e+13, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 5.76748, norm of the rhs: 7.89512e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 3.37194, norm of the rhs: 4.61585e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 2.37686, norm of the rhs: 3.25369e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 1.72592, norm of the rhs: 2.36261e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 1.31671, norm of the rhs: 1.80245e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 1.04168, norm of the rhs: 1.42595e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.864784, norm of the rhs: 1.18381e+12, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.659122, norm of the rhs: 9.02273e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 0.480804, norm of the rhs: 6.58174e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 15: 0.325347, norm of the rhs: 4.45369e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 16: 0.19051, norm of the rhs: 2.6079e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 17: 0.15074, norm of the rhs: 2.06348e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 18: 0.121945, norm of the rhs: 1.66931e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 19: 0.100443, norm of the rhs: 1.37497e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+36 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 20: 0.0844199, norm of the rhs: 1.15563e+11, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 21: 0.0717441, norm of the rhs: 9.82107e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 22: 0.0617118, norm of the rhs: 8.44774e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 23: 0.0536676, norm of the rhs: 7.34657e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 24: 0.0470753, norm of the rhs: 6.44415e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 25: 0.0415, norm of the rhs: 5.68095e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 26: 0.036938, norm of the rhs: 5.05646e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 27: 0.0329323, norm of the rhs: 4.50811e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 28: 0.0294087, norm of the rhs: 4.02577e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 29: 0.0258603, norm of the rhs: 3.54002e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 30: 0.0230077, norm of the rhs: 3.14953e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 31: 0.0205987, norm of the rhs: 2.81977e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 32: 0.0185379, norm of the rhs: 2.53765e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+36 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 33: 0.0168256, norm of the rhs: 2.30326e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+36 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 34: 0.0152702, norm of the rhs: 2.09034e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 35: 0.0138554, norm of the rhs: 1.89667e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 36: 0.0126018, norm of the rhs: 1.72506e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 37: 0.0114855, norm of the rhs: 1.57225e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 38: 0.0104873, norm of the rhs: 1.43561e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 39: 0.00959156, norm of the rhs: 1.31299e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 40: 0.00878527, norm of the rhs: 1.20262e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 41: 0.00805774, norm of the rhs: 1.10303e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 42: 0.00739987, norm of the rhs: 1.01297e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 43: 0.00680387, norm of the rhs: 9.31383e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 44: 0.00626299, norm of the rhs: 8.57342e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 45: 0.00577135, norm of the rhs: 7.90042e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 46: 0.0053238, norm of the rhs: 7.28776e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 47: 0.00491578, norm of the rhs: 6.72923e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 48: 0.0045433, norm of the rhs: 6.21933e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 49: 0.00420279, norm of the rhs: 5.75321e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 50: 0.0038911, norm of the rhs: 5.32654e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 51: 0.00360543, norm of the rhs: 4.93548e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 52: 0.00334327, norm of the rhs: 4.57661e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 53: 0.0031024, norm of the rhs: 4.24688e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 54: 0.00288083, norm of the rhs: 3.94358e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 55: 0.00267679, norm of the rhs: 3.66426e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 56: 0.00248867, norm of the rhs: 3.40674e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 57: 0.00231505, norm of the rhs: 3.16908e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 58: 0.00215466, norm of the rhs: 2.94952e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 59: 0.00200634, norm of the rhs: 2.74648e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 60: 0.00186905, norm of the rhs: 2.55855e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 61: 0.00174191, norm of the rhs: 2.3845e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 62: 0.00162423, norm of the rhs: 2.22341e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 63: 0.00151511, norm of the rhs: 2.07404e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 64: 0.00141378, norm of the rhs: 1.93532e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 65: 0.00131958, norm of the rhs: 1.80637e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 66: 0.00123195, norm of the rhs: 1.68642e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 67: 0.0011504, norm of the rhs: 1.57478e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 68: 0.00107446, norm of the rhs: 1.47083e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 69: 0.00100371, norm of the rhs: 1.37398e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 70: 0.000937767, norm of the rhs: 1.28371e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 71: 0.000876288, norm of the rhs: 1.19955e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 72: 0.000818947, norm of the rhs: 1.12106e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 73: 0.000765449, norm of the rhs: 1.04782e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 74: 0.000715521, norm of the rhs: 9.79478e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 75: 0.000668912, norm of the rhs: 9.15675e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 76: 0.00062539, norm of the rhs: 8.56098e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 77: 0.000584742, norm of the rhs: 8.00455e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 78: 0.00054677, norm of the rhs: 7.48475e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 79: 0.000511291, norm of the rhs: 6.99908e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 80: 0.000478137, norm of the rhs: 6.54522e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 81: 0.000447149, norm of the rhs: 6.12103e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 82: 0.000418183, norm of the rhs: 5.72451e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 83: 0.000391103, norm of the rhs: 5.35382e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 84: 0.000365784, norm of the rhs: 5.00723e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 85: 0.000342109, norm of the rhs: 4.68314e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 86: 0.00031997, norm of the rhs: 4.38008e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 87: 0.000299266, norm of the rhs: 4.09665e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 88: 0.000279902, norm of the rhs: 3.83158e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 89: 0.00026179, norm of the rhs: 3.58365e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 90: 0.00024485, norm of the rhs: 3.35176e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 91: 0.000229005, norm of the rhs: 3.13485e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 92: 0.000214183, norm of the rhs: 2.93195e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 93: 0.000200318, norm of the rhs: 2.74216e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 94: 0.000187348, norm of the rhs: 2.56462e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 95: 0.000175216, norm of the rhs: 2.39854e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 96: 0.000163867, norm of the rhs: 2.24318e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 97: 0.00015325, norm of the rhs: 2.09785e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 98: 0.000143319, norm of the rhs: 1.9619e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 99: 0.000134029, norm of the rhs: 1.83473e+08, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 100: 0.000125339, norm of the rhs: 1.71577e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.98e-05 m/year, 0.00015 m/year Temperature min/avg/max: 288 K, 890.2 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3416 W, 1.819e+04 W Writing graphical output: test-gravity/output/solution/solution-00000 *** Timestep 1: t=20000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 7 iterations. Solving crust_upper system ... 11 iterations. Solving crust_lower system ... 11 iterations. Solving mantle_lithosphere system ... 10 iterations. Initial Newton Stokes residual = 1.29656e+13, v = 1.29656e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00528621, norm of the rhs: 6.85391e+10 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.00121174, norm of the rhs: 1.57109e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.000875787, norm of the rhs: 1.13551e+10, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.000664427, norm of the rhs: 8.61471e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+42 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.000525444, norm of the rhs: 6.81272e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+40 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.000422542, norm of the rhs: 5.47852e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.000344595, norm of the rhs: 4.46789e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+41 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 0.000283134, norm of the rhs: 3.67101e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+40 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 0.000233958, norm of the rhs: 3.03341e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 0.000194095, norm of the rhs: 2.51657e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 0.00016155, norm of the rhs: 2.09459e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.000134827, norm of the rhs: 1.74812e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.000112792, norm of the rhs: 1.46242e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 9.45595e-05, norm of the rhs: 1.22602e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.95e-05 m/year, 0.000149 m/year Temperature min/avg/max: 288 K, 890.1 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3516 W, 1.814e+04 W Writing graphical output: test-gravity/output/solution/solution-00001 *** Timestep 2: t=40000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 6 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.2886e+13, v = 1.2886e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00182498, norm of the rhs: 2.35168e+10 Rebuilding Stokes preconditioner... Solving Stokes system... 200+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000279741, norm of the rhs: 3.60475e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.000156651, norm of the rhs: 2.01861e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.000103132, norm of the rhs: 1.32896e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 7.47732e-05, norm of the rhs: 9.63529e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.92e-05 m/year, 0.000148 m/year Temperature min/avg/max: 288 K, 890 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3610 W, 1.799e+04 W Writing graphical output: test-gravity/output/solution/solution-00002 *** Timestep 3: t=60000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 6 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.28435e+13, v = 1.28435e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 161+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00134965, norm of the rhs: 1.73342e+10 Rebuilding Stokes preconditioner... Solving Stokes system... 200+27 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000335285, norm of the rhs: 4.30623e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+36 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.000247691, norm of the rhs: 3.18122e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.000197093, norm of the rhs: 2.53137e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.000159103, norm of the rhs: 2.04344e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.000130164, norm of the rhs: 1.67177e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.000106981, norm of the rhs: 1.37401e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 8.83469e-05, norm of the rhs: 1.13468e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.89e-05 m/year, 0.000147 m/year Temperature min/avg/max: 288 K, 889.9 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3721 W, 1.789e+04 W Writing graphical output: test-gravity/output/solution/solution-00003 *** Timestep 4: t=80000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 6 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.28041e+13, v = 1.28041e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000683356, norm of the rhs: 8.74976e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000197636, norm of the rhs: 2.53055e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.000149211, norm of the rhs: 1.91051e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.000121365, norm of the rhs: 1.55397e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 9.8445e-05, norm of the rhs: 1.2605e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.86e-05 m/year, 0.000146 m/year Temperature min/avg/max: 288 K, 889.8 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3835 W, 1.78e+04 W Writing graphical output: test-gravity/output/solution/solution-00004 *** Timestep 5: t=100000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 6 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.27663e+13, v = 1.27663e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0004492, norm of the rhs: 5.73461e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000160109, norm of the rhs: 2.04399e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+40 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.000114925, norm of the rhs: 1.46717e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 9.2322e-05, norm of the rhs: 1.17861e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.84e-05 m/year, 0.000145 m/year Temperature min/avg/max: 288 K, 889.7 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -3944 W, 1.772e+04 W Writing graphical output: test-gravity/output/solution/solution-00005 *** Timestep 6: t=120000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.27305e+13, v = 1.27305e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 199+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000378073, norm of the rhs: 4.81305e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000135329, norm of the rhs: 1.72281e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 8.75055e-05, norm of the rhs: 1.11399e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.81e-05 m/year, 0.000145 m/year Temperature min/avg/max: 288 K, 889.6 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4044 W, 1.765e+04 W Writing graphical output: test-gravity/output/solution/solution-00006 *** Timestep 7: t=140000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.26958e+13, v = 1.26958e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 196+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000335223, norm of the rhs: 4.25593e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.00012178, norm of the rhs: 1.5461e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 7.48999e-05, norm of the rhs: 9.50914e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.79e-05 m/year, 0.000144 m/year Temperature min/avg/max: 288 K, 889.5 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4136 W, 1.759e+04 W Writing graphical output: test-gravity/output/solution/solution-00007 *** Timestep 8: t=160000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.26627e+13, v = 1.26627e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000317032, norm of the rhs: 4.01448e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000106904, norm of the rhs: 1.3537e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+40 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 5.90817e-05, norm of the rhs: 7.48133e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.76e-05 m/year, 0.000143 m/year Temperature min/avg/max: 288 K, 889.4 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4219 W, 1.753e+04 W Writing graphical output: test-gravity/output/solution/solution-00008 *** Timestep 9: t=180000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.26302e+13, v = 1.26302e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000267765, norm of the rhs: 3.38193e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.07304e-05, norm of the rhs: 1.14594e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.74e-05 m/year, 0.000142 m/year Temperature min/avg/max: 288 K, 889.3 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4295 W, 1.747e+04 W Writing graphical output: test-gravity/output/solution/solution-00009 *** Timestep 10: t=200000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 8 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.25981e+13, v = 1.25981e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 196+0 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000284415, norm of the rhs: 3.58307e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000109588, norm of the rhs: 1.3806e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 5.31733e-05, norm of the rhs: 6.69879e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.72e-05 m/year, 0.000141 m/year Temperature min/avg/max: 288 K, 889.2 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4364 W, 1.741e+04 W Writing graphical output: test-gravity/output/solution/solution-00010 *** Timestep 11: t=220000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.25673e+13, v = 1.25673e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000151406, norm of the rhs: 1.90276e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.55778e-05, norm of the rhs: 9.49808e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.69e-05 m/year, 0.000141 m/year Temperature min/avg/max: 288 K, 889.1 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4427 W, 1.736e+04 W Writing graphical output: test-gravity/output/solution/solution-00011 *** Timestep 12: t=240000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.25364e+13, v = 1.25364e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000305771, norm of the rhs: 3.83326e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000100838, norm of the rhs: 1.26414e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 4.84219e-05, norm of the rhs: 6.07035e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.67e-05 m/year, 0.00014 m/year Temperature min/avg/max: 288 K, 889 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4485 W, 1.731e+04 W Writing graphical output: test-gravity/output/solution/solution-00012 *** Timestep 13: t=260000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.25066e+13, v = 1.25066e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000201865, norm of the rhs: 2.52463e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.44019e-05, norm of the rhs: 9.30514e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.65e-05 m/year, 0.000139 m/year Temperature min/avg/max: 288 K, 889 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4539 W, 1.727e+04 W Writing graphical output: test-gravity/output/solution/solution-00013 *** Timestep 14: t=280000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.24765e+13, v = 1.24765e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000170792, norm of the rhs: 2.13089e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.4337e-05, norm of the rhs: 1.177e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.63e-05 m/year, 0.000138 m/year Temperature min/avg/max: 288 K, 888.9 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4589 W, 1.722e+04 W Writing graphical output: test-gravity/output/solution/solution-00014 *** Timestep 15: t=300000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.24473e+13, v = 1.24473e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+3 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000472735, norm of the rhs: 5.88428e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.000100495, norm of the rhs: 1.25089e+09, newton_derivative_scaling_factor: 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 4.83408e-05, norm of the rhs: 6.01713e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.6e-05 m/year, 0.000138 m/year Temperature min/avg/max: 288 K, 888.8 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4635 W, 1.718e+04 W Writing graphical output: test-gravity/output/solution/solution-00015 *** Timestep 16: t=320000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.24188e+13, v = 1.24188e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000132116, norm of the rhs: 1.64073e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.48196e-05, norm of the rhs: 8.04981e+08, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.58e-05 m/year, 0.000137 m/year Temperature min/avg/max: 288 K, 888.7 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4679 W, 1.714e+04 W Writing graphical output: test-gravity/output/solution/solution-00016 *** Timestep 17: t=340000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.239e+13, v = 1.239e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000350317, norm of the rhs: 4.34044e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.13371e-05, norm of the rhs: 1.13167e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.56e-05 m/year, 0.000136 m/year Temperature min/avg/max: 288 K, 888.6 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4719 W, 1.71e+04 W Writing graphical output: test-gravity/output/solution/solution-00017 *** Timestep 18: t=360000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 5 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.23617e+13, v = 1.23617e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000256276, norm of the rhs: 3.16801e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.98192e-05, norm of the rhs: 1.23394e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.54e-05 m/year, 0.000136 m/year Temperature min/avg/max: 288 K, 888.5 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4757 W, 1.706e+04 W Writing graphical output: test-gravity/output/solution/solution-00018 *** Timestep 19: t=380000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 4 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.23338e+13, v = 1.23338e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000222126, norm of the rhs: 2.73967e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.48592e-05, norm of the rhs: 1.16998e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.51e-05 m/year, 0.000135 m/year Temperature min/avg/max: 288 K, 888.4 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4793 W, 1.703e+04 W Writing graphical output: test-gravity/output/solution/solution-00019 *** Timestep 20: t=400000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 1 iterations. Solving temperature system... 4 iterations. Solving crust_upper system ... 7 iterations. Solving crust_lower system ... 8 iterations. Solving mantle_lithosphere system ... 7 iterations. Initial Newton Stokes residual = 1.23063e+13, v = 1.23063e+13, p = 0 Rebuilding Stokes preconditioner... Solving Stokes system... 200+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000607853, norm of the rhs: 7.48042e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 200+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.36782e-05, norm of the rhs: 1.15283e+09, newton_derivative_scaling_factor: 0 Postprocessing: RMS, max velocity: 4.49e-05 m/year, 0.000134 m/year Temperature min/avg/max: 288 K, 888.3 K, 1271 K Heat fluxes through boundary parts: 0 W, 0 W, -4827 W, 1.699e+04 W Writing graphical output: test-gravity/output/solution/solution-00020 Termination requested by criterion: end time +----------------------------------------------+------------+------------+ | Total wallclock time elapsed since start | 583s | | | | | | | Section | no. calls | wall time | % of total | +----------------------------------+-----------+------------+------------+ | Assemble Stokes system | 21 | 5.75s | 0.99% | | Assemble Stokes system Picard | 151 | 140s | 24% | | Assemble Stokes system rhs | 172 | 35.4s | 6.1% | | Assemble composition system | 63 | 31s | 5.3% | | Assemble temperature system | 21 | 11.4s | 2% | | Build Stokes preconditioner | 172 | 131s | 23% | | Build composition preconditioner | 63 | 0.129s | 0% | | Build temperature preconditioner | 21 | 0.0404s | 0% | | Initialization | 1 | 0.444s | 0% | | Mesh deformation | 21 | 2.26s | 0.39% | | Mesh deformation initialize | 1 | 0.061s | 0% | | Postprocessing | 21 | 11.8s | 2% | | Setup dof systems | 1 | 0.186s | 0% | | Setup initial conditions | 1 | 0.473s | 0% | | Setup matrices | 1 | 0.297s | 0% | | Solve Stokes system | 172 | 196s | 34% | | Solve composition system | 63 | 0.152s | 0% | | Solve temperature system | 21 | 0.0467s | 0% | +----------------------------------+-----------+------------+------------+ -- Total wallclock time elapsed including restarts: 583s ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.6.0-pre&NewtonSolver=1&sha=e3a102e7b&src=code -----------------------------------------------------------------------------