----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.4.0 -- . using deal.II 9.3.0 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in OPTIMIZED mode -- . running with 32 MPI processes ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.4.0&NewtonSolver=1&sha=&src=code ----------------------------------------------------------------------------- Number of active cells: 26,400 (on 3 levels) Number of degrees of freedom: 876,853 (212,530+26,733+106,265+106,265+106,265+106,265+106,265+106,265) Number of mesh deformation degrees of freedom: 53466 Solving mesh displacement system... 0 iterations. *** Timestep 0: t=0 years, dt=0 years Solving mesh displacement system... 0 iterations. Solving temperature system... 0 iterations. Skipping noninitial_plastic_strain composition solve because RHS is zero. Solving plastic_strain 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 = 3.13029e+15, v = 3.13029e+15, p = 1.17873e+10 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 1, norm of the rhs: 3.13029e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.0104314, norm of the rhs: 3.26533e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.00479001, norm of the rhs: 1.49941e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.00457991, norm of the rhs: 1.43364e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.00335017, norm of the rhs: 1.0487e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.0036867, norm of the rhs: 1.15404e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.00416366, norm of the rhs: 1.30334e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 0.00460419, norm of the rhs: 1.44124e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 0.00275789, norm of the rhs: 8.63298e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 0.00180066, norm of the rhs: 5.63658e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 0.00153968, norm of the rhs: 4.81963e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.00135493, norm of the rhs: 4.24133e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.00173425, norm of the rhs: 5.42868e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 0.00122446, norm of the rhs: 3.83291e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 15: 0.00107937, norm of the rhs: 3.37875e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 16: 0.00102736, norm of the rhs: 3.21592e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 17: 0.0010277, norm of the rhs: 3.21699e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 18: 0.000811313, norm of the rhs: 2.53964e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 19: 0.000814228, norm of the rhs: 2.54877e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 20: 0.000605072, norm of the rhs: 1.89405e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 21: 0.000525556, norm of the rhs: 1.64514e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 22: 0.000467298, norm of the rhs: 1.46278e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 23: 0.000423577, norm of the rhs: 1.32592e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 24: 0.000369614, norm of the rhs: 1.157e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 25: 0.000331565, norm of the rhs: 1.03789e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 26: 0.000295383, norm of the rhs: 9.24634e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 27: 0.000272253, norm of the rhs: 8.52229e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 28: 0.000244877, norm of the rhs: 7.66535e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 29: 0.000217487, norm of the rhs: 6.80798e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 30: 0.000195658, norm of the rhs: 6.12465e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 31: 0.000183774, norm of the rhs: 5.75266e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 32: 0.00017067, norm of the rhs: 5.34246e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 33: 0.000155254, norm of the rhs: 4.8599e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 34: 0.00013887, norm of the rhs: 4.34702e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 35: 0.000131476, norm of the rhs: 4.11559e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 36: 0.000125464, norm of the rhs: 3.9274e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 37: 0.000121031, norm of the rhs: 3.7886e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 38: 0.000113917, norm of the rhs: 3.56592e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 39: 0.00010755, norm of the rhs: 3.36662e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 40: 0.00010286, norm of the rhs: 3.21982e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 41: 9.88349e-05, norm of the rhs: 3.09382e+11, newton_derivative_scaling_factor: 0 Number of active cells: 31,920 (on 4 levels) Number of degrees of freedom: 1,063,943 (257,874+32,447+128,937+128,937+128,937+128,937+128,937+128,937) Number of mesh deformation degrees of freedom: 64894 Solving mesh displacement system... 0 iterations. *** Timestep 0: t=0 years, dt=0 years Solving mesh displacement system... 0 iterations. Solving temperature system... 0 iterations. Skipping noninitial_plastic_strain composition solve because RHS is zero. Solving plastic_strain 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 = 3.13029e+15, v = 3.13029e+15, p = 1.1787e+10 Rebuilding Stokes preconditioner... Solving Stokes system... 0+24 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 1, norm of the rhs: 3.13029e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 0.00808083, norm of the rhs: 2.52953e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 3: 0.00381677, norm of the rhs: 1.19476e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 4: 0.00441869, norm of the rhs: 1.38318e+13, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 5: 0.00238546, norm of the rhs: 7.46717e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 6: 0.00159647, norm of the rhs: 4.9974e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 7: 0.00203518, norm of the rhs: 6.37071e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 8: 0.00174116, norm of the rhs: 5.45033e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 9: 0.00121163, norm of the rhs: 3.79276e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 10: 0.000858767, norm of the rhs: 2.68819e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 11: 0.000860246, norm of the rhs: 2.69282e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 12: 0.000878593, norm of the rhs: 2.75025e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 13: 0.00060987, norm of the rhs: 1.90907e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 14: 0.000678542, norm of the rhs: 2.12403e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+2 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 15: 0.000573307, norm of the rhs: 1.79461e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 16: 0.000495685, norm of the rhs: 1.55163e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 17: 0.000422653, norm of the rhs: 1.32303e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 18: 0.00037736, norm of the rhs: 1.18125e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 19: 0.000350763, norm of the rhs: 1.09799e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 20: 0.000333588, norm of the rhs: 1.04423e+12, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 21: 0.000316254, norm of the rhs: 9.89965e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 22: 0.000300962, norm of the rhs: 9.42096e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 23: 0.000285631, norm of the rhs: 8.94108e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 24: 0.000270882, norm of the rhs: 8.47939e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 25: 0.00025543, norm of the rhs: 7.9957e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 26: 0.000242451, norm of the rhs: 7.5894e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 27: 0.000228167, norm of the rhs: 7.14227e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 28: 0.00021553, norm of the rhs: 6.7467e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 29: 0.000206912, norm of the rhs: 6.47692e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 30: 0.000200885, norm of the rhs: 6.28829e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 31: 0.000188223, norm of the rhs: 5.89191e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 32: 0.000181456, norm of the rhs: 5.68009e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 33: 0.000171437, norm of the rhs: 5.36646e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 34: 0.000161992, norm of the rhs: 5.07083e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 35: 0.000148449, norm of the rhs: 4.64688e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 36: 0.00013698, norm of the rhs: 4.28786e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 37: 0.000132026, norm of the rhs: 4.13279e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 38: 0.000126977, norm of the rhs: 3.97475e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 39: 0.000122895, norm of the rhs: 3.84697e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 40: 0.000118255, norm of the rhs: 3.70173e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 41: 0.000114512, norm of the rhs: 3.58455e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 42: 0.00011037, norm of the rhs: 3.4549e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 43: 0.000107543, norm of the rhs: 3.3664e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 44: 0.000102251, norm of the rhs: 3.20075e+11, newton_derivative_scaling_factor: 0 The linear solver tolerance is set to 0.01. Rebuilding Stokes preconditioner... Solving Stokes system... 0+1 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 45: 9.84162e-05, norm of the rhs: 3.08071e+11, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: 0/0/0 // 0/1.5/4.55e+10 // -0.125/1/2.045e+10 // -0.125/1.125/2e+10 // 0/1.125/6.195e+11 Heat flux densities for boundaries: -2.958 W/m, -2.939 W/m, -0.00112 W/m, 0.2464 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -1.657e+06 W, -1.646e+06 W, -1120 W, 2.464e+05 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -3.259e+07 kg/yr, -3.238e+07 kg/yr, 3.34e+07 kg/yr, 2.431e+07 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -2.963e+07 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 273 K, 1630 K, 1893 K Topography min/max: 0 m, 0 m RMS, max velocity: 0.018 m/year, 0.0308 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00000 *** Timestep 1: t=20000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 10 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.64829e+17, v = 2.64829e+17, p = 9.16372e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+12 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.423491, norm of the rhs: 1.12153e+17 Rebuilding Stokes preconditioner... Solving Stokes system... 0+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 5.01431e-05, norm of the rhs: 1.32793e+13, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0005049/0.004241/7.395e+06 // -0.01359/1.493/4.551e+10 // -0.1065/1.021/2.045e+10 // -0.103/1.107/2e+10 // -0.0206/1.103/6.197e+11 Heat flux densities for boundaries: 2.208 W/m, 2.177 W/m, -0.001346 W/m, -0.128 W/m, 0.4083 W/m, 0.4083 W/m Heat fluxes through boundary parts: 1.236e+06 W, 1.219e+06 W, -1346 W, -1.28e+05 W, 4.083e+04 W, 4.083e+04 W Mass fluxes through boundary parts: 2.443e+07 kg/yr, 2.409e+07 kg/yr, -2.431e+07 kg/yr, -2.329e+07 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.052e+08 Pa, 9.981e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 273 K, 1630 K, 1893 K Topography min/max: -214.9 m, 579.9 m RMS, max velocity: 0.0137 m/year, 0.0245 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00001 *** Timestep 2: t=40000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 12 iterations. Solving noninitial_plastic_strain system ... 28 iterations. Solving plastic_strain system ... 27 iterations. Solving crust_upper system ... 24 iterations. Solving crust_lower system ... 26 iterations. Solving mantle_lithosphere system ... 23 iterations. Initial Newton Stokes residual = 2.69053e+17, v = 2.69053e+17, p = 8.70733e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+14 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0516191, norm of the rhs: 1.38883e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+48 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 2.69172e-05, norm of the rhs: 7.24216e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001672/0.0183/2.088e+07 // -0.01818/1.599/4.578e+10 // -0.1113/1.019/2.068e+10 // -0.109/1.111/1.999e+10 // -0.02044/1.109/6.193e+11 Heat flux densities for boundaries: -1.189 W/m, -1.18 W/m, -0.001135 W/m, 0.1265 W/m, 0.4097 W/m, 0.4097 W/m Heat fluxes through boundary parts: -6.658e+05 W, -6.608e+05 W, -1135 W, 1.265e+05 W, 4.097e+04 W, 4.097e+04 W Mass fluxes through boundary parts: -1.312e+07 kg/yr, -1.302e+07 kg/yr, 1.306e+07 kg/yr, 8.51e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.432e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 260.2 K, 1630 K, 1893 K Topography min/max: -96.87 m, 93.03 m RMS, max velocity: 0.00687 m/year, 0.0131 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00002 *** Timestep 3: t=60000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 10 iterations. Solving noninitial_plastic_strain system ... 20 iterations. Solving plastic_strain system ... 19 iterations. Solving crust_upper system ... 20 iterations. Solving crust_lower system ... 19 iterations. Solving mantle_lithosphere system ... 16 iterations. Initial Newton Stokes residual = 2.63194e+17, v = 2.63194e+17, p = 8.80134e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+13 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0674617, norm of the rhs: 1.77555e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+38 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 3.02054e-05, norm of the rhs: 7.94988e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001465/0.03099/3.245e+07 // -0.02156/1.509/4.571e+10 // -0.1002/1.019/2.06e+10 // -0.09699/1.1/1.999e+10 // -0.01684/1.097/6.195e+11 Heat flux densities for boundaries: 0.05595 W/m, 0.06326 W/m, -0.001298 W/m, 0.01434 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: 3.133e+04 W, 3.542e+04 W, -1298 W, 1.434e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: 6.335e+05 kg/yr, 7.152e+05 kg/yr, -7.65e+05 kg/yr, -3.151e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -7.373e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 271 K, 1630 K, 1893 K Topography min/max: -175.6 m, 260.9 m RMS, max velocity: 0.00182 m/year, 0.00534 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00003 *** Timestep 4: t=80000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 10 iterations. Solving noninitial_plastic_strain system ... 17 iterations. Solving plastic_strain system ... 16 iterations. Solving crust_upper system ... 16 iterations. Solving crust_lower system ... 16 iterations. Solving mantle_lithosphere system ... 14 iterations. Initial Newton Stokes residual = 2.64936e+17, v = 2.64936e+17, p = 8.80645e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+12 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0204192, norm of the rhs: 5.40978e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+35 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.7608e-05, norm of the rhs: 4.66499e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001179/0.03815/4.185e+07 // -0.02144/1.546/4.576e+10 // -0.09682/1.023/2.063e+10 // -0.09402/1.097/1.998e+10 // -0.02123/1.094/6.194e+11 Heat flux densities for boundaries: 0.01612 W/m, 0.01114 W/m, -0.001228 W/m, 0.03716 W/m, 0.4094 W/m, 0.4094 W/m Heat fluxes through boundary parts: 9028 W, 6239 W, -1228 W, 3.716e+04 W, 4.094e+04 W, 4.094e+04 W Mass fluxes through boundary parts: 1.882e+05 kg/yr, 1.334e+05 kg/yr, -3.039e+05 kg/yr, -2.668e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.813e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 267.2 K, 1630 K, 1893 K Topography min/max: -181.2 m, 163.1 m RMS, max velocity: 0.00172 m/year, 0.00534 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00004 *** Timestep 5: t=100000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.64589e+17, v = 2.64589e+17, p = 8.84585e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000498256, norm of the rhs: 1.31833e+14 Rebuilding Stokes preconditioner... Solving Stokes system... 0+24 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.71285e-05, norm of the rhs: 4.53202e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008713/0.03959/4.946e+07 // -0.01834/1.549/4.577e+10 // -0.09337/1.026/2.063e+10 // -0.09089/1.093/1.998e+10 // -0.02434/1.091/6.194e+11 Heat flux densities for boundaries: -0.6239 W/m, -0.631 W/m, -0.001219 W/m, 0.07505 W/m, 0.4095 W/m, 0.4095 W/m Heat fluxes through boundary parts: -3.494e+05 W, -3.533e+05 W, -1219 W, 7.505e+04 W, 4.095e+04 W, 4.095e+04 W Mass fluxes through boundary parts: -6.871e+06 kg/yr, -6.949e+06 kg/yr, 6.937e+06 kg/yr, 3.216e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.081e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.4 K, 1630 K, 1893 K Topography min/max: -184.5 m, 109.1 m RMS, max velocity: 0.00366 m/year, 0.00732 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00005 *** Timestep 6: t=120000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63785e+17, v = 2.63785e+17, p = 8.83543e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0553528, norm of the rhs: 1.46012e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.7141e-05, norm of the rhs: 4.52153e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.00117/0.03951/5.668e+07 // -0.01332/1.514/4.575e+10 // -0.0881/1.027/2.06e+10 // -0.08544/1.088/1.997e+10 // -0.02496/1.086/6.195e+11 Heat flux densities for boundaries: -0.3924 W/m, -0.3923 W/m, -0.001231 W/m, 0.05946 W/m, 0.4093 W/m, 0.4093 W/m Heat fluxes through boundary parts: -2.198e+05 W, -2.197e+05 W, -1231 W, 5.946e+04 W, 4.093e+04 W, 4.093e+04 W Mass fluxes through boundary parts: -4.305e+06 kg/yr, -4.302e+06 kg/yr, 4.398e+06 kg/yr, 9.358e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.575e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 270.6 K, 1630 K, 1893 K Topography min/max: -227.3 m, 163.7 m RMS, max velocity: 0.00248 m/year, 0.00548 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00006 *** Timestep 7: t=140000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63747e+17, v = 2.63747e+17, p = 8.82179e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0048175, norm of the rhs: 1.2706e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.34447e-05, norm of the rhs: 3.54601e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001029/0.03966/6.344e+07 // -0.01346/1.504/4.576e+10 // -0.08443/1.029/2.06e+10 // -0.08175/1.084/1.997e+10 // -0.02678/1.082/6.195e+11 Heat flux densities for boundaries: 0.1261 W/m, 0.1347 W/m, -0.001236 W/m, 0.0258 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: 7.063e+04 W, 7.546e+04 W, -1236 W, 2.58e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: 1.416e+06 kg/yr, 1.512e+06 kg/yr, -1.495e+06 kg/yr, -3.874e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.041e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.4 K, 1630 K, 1893 K Topography min/max: -257.1 m, 185.2 m RMS, max velocity: 0.00185 m/year, 0.00537 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00007 *** Timestep 8: t=160000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63683e+17, v = 2.63683e+17, p = 8.8203e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0393373, norm of the rhs: 1.03726e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.50195e-05, norm of the rhs: 3.96038e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001236/0.04052/7.007e+07 // -0.01437/1.521/4.579e+10 // -0.08257/1.031/2.062e+10 // -0.08032/1.082/1.997e+10 // -0.02958/1.08/6.194e+11 Heat flux densities for boundaries: -0.207 W/m, -0.1927 W/m, -0.001219 W/m, 0.04785 W/m, 0.4094 W/m, 0.4094 W/m Heat fluxes through boundary parts: -1.159e+05 W, -1.079e+05 W, -1219 W, 4.785e+04 W, 4.094e+04 W, 4.094e+04 W Mass fluxes through boundary parts: -2.264e+06 kg/yr, -2.105e+06 kg/yr, 2.187e+06 kg/yr, -7.935e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -7.024e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.9 K, 1630 K, 1893 K Topography min/max: -252.3 m, 151.7 m RMS, max velocity: 0.00169 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00008 *** Timestep 9: t=180000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63229e+17, v = 2.63229e+17, p = 8.82678e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0370522, norm of the rhs: 9.75321e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.28031e-05, norm of the rhs: 3.37015e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001241/0.04056/7.66e+07 // -0.01482/1.506/4.579e+10 // -0.07953/1.032/2.061e+10 // -0.07786/1.08/1.996e+10 // -0.03079/1.078/6.195e+11 Heat flux densities for boundaries: -0.6173 W/m, -0.6071 W/m, -0.001213 W/m, 0.07444 W/m, 0.4094 W/m, 0.4094 W/m Heat fluxes through boundary parts: -3.457e+05 W, -3.4e+05 W, -1213 W, 7.444e+04 W, 4.094e+04 W, 4.094e+04 W Mass fluxes through boundary parts: -6.785e+06 kg/yr, -6.67e+06 kg/yr, 6.889e+06 kg/yr, 2.946e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.655e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 1893 K Topography min/max: -266.7 m, 140.8 m RMS, max velocity: 0.00386 m/year, 0.00707 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00009 *** Snapshot created! *** Timestep 10: t=200000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.62874e+17, v = 2.62874e+17, p = 8.81262e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0883203, norm of the rhs: 2.32172e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.30012e-05, norm of the rhs: 3.41768e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001185/0.04006/8.307e+07 // -0.01535/1.479/4.578e+10 // -0.07572/1.033/2.058e+10 // -0.07409/1.076/1.996e+10 // -0.03081/1.074/6.195e+11 Heat flux densities for boundaries: -0.07297 W/m, -0.07302 W/m, -0.001232 W/m, 0.03836 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -4.086e+04 W, -4.089e+04 W, -1232 W, 3.836e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -7.761e+05 kg/yr, -7.762e+05 kg/yr, 8.353e+05 kg/yr, -2.042e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.938e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.9 K, 1630 K, 1893 K Topography min/max: -309.4 m, 196.5 m RMS, max velocity: 0.00143 m/year, 0.00534 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00010 *** Timestep 11: t=220000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63024e+17, v = 2.63024e+17, p = 8.80993e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.023009, norm of the rhs: 6.05191e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.20675e-05, norm of the rhs: 3.17404e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001293/0.03959/8.959e+07 // -0.01608/1.49/4.58e+10 // -0.07374/1.035/2.06e+10 // -0.07231/1.074/1.995e+10 // -0.03305/1.072/6.195e+11 Heat flux densities for boundaries: 0.06848 W/m, 0.0588 W/m, -0.001223 W/m, 0.03101 W/m, 0.4093 W/m, 0.4093 W/m Heat fluxes through boundary parts: 3.835e+04 W, 3.293e+04 W, -1223 W, 3.101e+04 W, 4.093e+04 W, 4.093e+04 W Mass fluxes through boundary parts: 7.768e+05 kg/yr, 6.694e+05 kg/yr, -7.772e+05 kg/yr, -3.223e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -4.746e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.7 K, 1630 K, 1893 K Topography min/max: -316.1 m, 176.5 m RMS, max velocity: 0.00161 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00011 *** Timestep 12: t=240000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 13 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.62868e+17, v = 2.62868e+17, p = 8.83674e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0179123, norm of the rhs: 4.70857e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+20 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.12351e-05, norm of the rhs: 2.95335e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001329/0.03857/9.603e+07 // -0.01637/1.499/4.582e+10 // -0.0721/1.036/2.06e+10 // -0.07088/1.073/1.995e+10 // -0.03466/1.07/6.195e+11 Heat flux densities for boundaries: -0.5423 W/m, -0.5516 W/m, -0.001204 W/m, 0.07128 W/m, 0.4094 W/m, 0.4094 W/m Heat fluxes through boundary parts: -3.037e+05 W, -3.089e+05 W, -1204 W, 7.129e+04 W, 4.094e+04 W, 4.094e+04 W Mass fluxes through boundary parts: -5.963e+06 kg/yr, -6.066e+06 kg/yr, 6.095e+06 kg/yr, 2.407e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.114e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.1 K, 1630 K, 1893 K Topography min/max: -311.9 m, 118.3 m RMS, max velocity: 0.00337 m/year, 0.00631 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00012 *** Timestep 13: t=260000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.62406e+17, v = 2.62406e+17, p = 8.82298e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.104788, norm of the rhs: 2.7497e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.28998e-05, norm of the rhs: 3.38498e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001368/0.03801/1.024e+08 // -0.0164/1.475/4.58e+10 // -0.06912/1.036/2.058e+10 // -0.06789/1.07/1.995e+10 // -0.03429/1.067/6.195e+11 Heat flux densities for boundaries: -0.3323 W/m, -0.3317 W/m, -0.00122 W/m, 0.05577 W/m, 0.4093 W/m, 0.4093 W/m Heat fluxes through boundary parts: -1.861e+05 W, -1.857e+05 W, -1220 W, 5.577e+04 W, 4.093e+04 W, 4.093e+04 W Mass fluxes through boundary parts: -3.639e+06 kg/yr, -3.632e+06 kg/yr, 3.732e+06 kg/yr, 3.65e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.822e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 270.4 K, 1630 K, 1893 K Topography min/max: -349.4 m, 190.8 m RMS, max velocity: 0.00219 m/year, 0.00543 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00013 *** Timestep 14: t=280000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.6238e+17, v = 2.6238e+17, p = 8.80979e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00908333, norm of the rhs: 2.38329e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+21 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.75493e-06, norm of the rhs: 2.5595e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001326/0.03892/1.088e+08 // -0.02012/1.47/4.581e+10 // -0.06689/1.037/2.057e+10 // -0.06569/1.068/1.994e+10 // -0.03517/1.065/6.195e+11 Heat flux densities for boundaries: 0.1243 W/m, 0.13 W/m, -0.001225 W/m, 0.02666 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: 6.963e+04 W, 7.281e+04 W, -1225 W, 2.666e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: 1.393e+06 kg/yr, 1.457e+06 kg/yr, -1.486e+06 kg/yr, -3.816e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.584e+07 Pa, 9.982e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.8 K, 1630 K, 1893 K Topography min/max: -373 m, 216.3 m RMS, max velocity: 0.00187 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00014 *** Timestep 15: t=300000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.62426e+17, v = 2.62426e+17, p = 8.83195e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0594834, norm of the rhs: 1.561e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.1941e-05, norm of the rhs: 3.13364e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.00122/0.03988/1.152e+08 // -0.02514/1.485/4.584e+10 // -0.06583/1.038/2.059e+10 // -0.06491/1.067/1.994e+10 // -0.03649/1.064/6.195e+11 Heat flux densities for boundaries: -0.3453 W/m, -0.3404 W/m, -0.001202 W/m, 0.05863 W/m, 0.4094 W/m, 0.4094 W/m Heat fluxes through boundary parts: -1.934e+05 W, -1.906e+05 W, -1202 W, 5.863e+04 W, 4.094e+04 W, 4.094e+04 W Mass fluxes through boundary parts: -3.792e+06 kg/yr, -3.738e+06 kg/yr, 3.77e+06 kg/yr, 5.535e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.526e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.2 K, 1630 K, 1893 K Topography min/max: -364.6 m, 137 m RMS, max velocity: 0.00218 m/year, 0.00543 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00015 *** Timestep 16: t=320000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.62065e+17, v = 2.62065e+17, p = 8.80595e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.111239, norm of the rhs: 2.91518e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.95685e-06, norm of the rhs: 2.60935e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001196/0.04081/1.216e+08 // -0.03008/1.469/4.583e+10 // -0.06419/1.037/2.057e+10 // -0.06328/1.065/1.993e+10 // -0.03614/1.062/6.195e+11 Heat flux densities for boundaries: -0.5595 W/m, -0.5665 W/m, -0.001204 W/m, 0.07195 W/m, 0.4093 W/m, 0.4093 W/m Heat fluxes through boundary parts: -3.133e+05 W, -3.172e+05 W, -1204 W, 7.196e+04 W, 4.093e+04 W, 4.093e+04 W Mass fluxes through boundary parts: -6.147e+06 kg/yr, -6.225e+06 kg/yr, 6.319e+06 kg/yr, 2.509e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.193e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 270 K, 1630 K, 1893 K Topography min/max: -387.3 m, 149.1 m RMS, max velocity: 0.00336 m/year, 0.00639 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00016 *** Timestep 17: t=340000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.61926e+17, v = 2.61926e+17, p = 8.78482e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.121061, norm of the rhs: 3.17091e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.06441e-05, norm of the rhs: 2.78797e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001251/0.04126/1.279e+08 // -0.03448/1.468/4.582e+10 // -0.06243/1.037/2.055e+10 // -0.06103/1.064/1.993e+10 // -0.0356/1.06/6.195e+11 Heat flux densities for boundaries: -0.04527 W/m, -0.06492 W/m, -0.001218 W/m, 0.03865 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -2.535e+04 W, -3.636e+04 W, -1218 W, 3.865e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -4.718e+05 kg/yr, -6.902e+05 kg/yr, 6.127e+05 kg/yr, -2.185e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.7e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 1893 K Topography min/max: -425.9 m, 206.9 m RMS, max velocity: 0.00146 m/year, 0.00534 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00017 *** Timestep 18: t=360000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.6208e+17, v = 2.6208e+17, p = 8.79265e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0623198, norm of the rhs: 1.63328e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.88115e-06, norm of the rhs: 2.58965e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001273/0.04142/1.344e+08 // -0.0377/1.466/4.585e+10 // -0.06127/1.038/2.057e+10 // -0.05987/1.062/1.993e+10 // -0.03675/1.059/6.195e+11 Heat flux densities for boundaries: -0.03139 W/m, -0.04528 W/m, -0.001207 W/m, 0.03857 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -1.758e+04 W, -2.536e+04 W, -1207 W, 3.857e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -3.249e+05 kg/yr, -4.792e+05 kg/yr, 3.714e+05 kg/yr, -2.284e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -4.735e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.4 K, 1630 K, 1893 K Topography min/max: -429.5 m, 154.5 m RMS, max velocity: 0.00149 m/year, 0.00534 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00018 *** Timestep 19: t=380000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.6192e+17, v = 2.6192e+17, p = 8.79434e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+4 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00180371, norm of the rhs: 4.72427e+14 Rebuilding Stokes preconditioner... Solving Stokes system... 0+22 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.16081e-06, norm of the rhs: 2.39939e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001273/0.04125/1.407e+08 // -0.0402/1.469/4.587e+10 // -0.06035/1.038/2.057e+10 // -0.05886/1.061/1.992e+10 // -0.03704/1.058/6.195e+11 Heat flux densities for boundaries: -0.549 W/m, -0.5405 W/m, -0.001193 W/m, 0.07146 W/m, 0.4093 W/m, 0.4093 W/m Heat fluxes through boundary parts: -3.074e+05 W, -3.027e+05 W, -1193 W, 7.147e+04 W, 4.093e+04 W, 4.093e+04 W Mass fluxes through boundary parts: -6.033e+06 kg/yr, -5.939e+06 kg/yr, 6.109e+06 kg/yr, 2.346e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -5.877e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.6 K, 1630 K, 1893 K Topography min/max: -430.9 m, 151.7 m RMS, max velocity: 0.0033 m/year, 0.00624 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00019 *** Snapshot created! *** Timestep 20: t=400000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61556e+17, v = 2.61556e+17, p = 8.75977e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.145691, norm of the rhs: 3.81064e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.09981e-05, norm of the rhs: 2.87663e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001259/0.04097/1.47e+08 // -0.04214/1.457/4.585e+10 // -0.05912/1.037/2.055e+10 // -0.05755/1.06/1.992e+10 // -0.03581/1.057/6.195e+11 Heat flux densities for boundaries: -0.386 W/m, -0.3662 W/m, -0.001204 W/m, 0.05962 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -2.161e+05 W, -2.051e+05 W, -1204 W, 5.962e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -4.227e+06 kg/yr, -4.007e+06 kg/yr, 4.28e+06 kg/yr, 7.187e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.998e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 1893 K Topography min/max: -467.8 m, 176 m RMS, max velocity: 0.00243 m/year, 0.00548 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00020 *** Timestep 21: t=420000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61718e+17, v = 2.61718e+17, p = 8.73459e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0214676, norm of the rhs: 5.61847e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.97799e-06, norm of the rhs: 2.08798e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001195/0.04046/1.533e+08 // -0.04278/1.451/4.586e+10 // -0.05758/1.036/2.054e+10 // -0.05703/1.058/1.991e+10 // -0.03562/1.057/6.196e+11 Heat flux densities for boundaries: 0.1336 W/m, 0.1399 W/m, -0.001214 W/m, 0.02602 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: 7.483e+04 W, 7.833e+04 W, -1214 W, 2.602e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: 1.502e+06 kg/yr, 1.572e+06 kg/yr, -1.535e+06 kg/yr, -3.965e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.151e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.1 K, 1630 K, 1893 K Topography min/max: -494.1 m, 187 m RMS, max velocity: 0.00186 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00021 *** Timestep 22: t=440000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61772e+17, v = 2.61772e+17, p = 8.74365e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0714163, norm of the rhs: 1.86948e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+27 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.7213e-06, norm of the rhs: 2.54476e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001138/0.04005/1.598e+08 // -0.0446/1.458/4.589e+10 // -0.05666/1.037/2.056e+10 // -0.05604/1.057/1.991e+10 // -0.03666/1.056/6.195e+11 Heat flux densities for boundaries: -0.2573 W/m, -0.2757 W/m, -0.001193 W/m, 0.05334 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -1.441e+05 W, -1.544e+05 W, -1193 W, 5.335e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -2.813e+06 kg/yr, -3.018e+06 kg/yr, 2.979e+06 kg/yr, -2.211e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.658e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.9 K, 1630 K, 1893 K Topography min/max: -487.8 m, 151.7 m RMS, max velocity: 0.00193 m/year, 0.0054 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00022 *** Timestep 23: t=460000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61467e+17, v = 2.61467e+17, p = 8.74012e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0829498, norm of the rhs: 2.16886e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.12191e-06, norm of the rhs: 2.12361e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001052/0.03947/1.66e+08 // -0.04504/1.445/4.588e+10 // -0.0557/1.036/2.055e+10 // -0.05526/1.056/1.991e+10 // -0.03604/1.055/6.195e+11 Heat flux densities for boundaries: -0.55 W/m, -0.5735 W/m, -0.001193 W/m, 0.07157 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -3.08e+05 W, -3.211e+05 W, -1193 W, 7.157e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -6.035e+06 kg/yr, -6.296e+06 kg/yr, 6.369e+06 kg/yr, 2.431e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.131e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 1893 K Topography min/max: -510.2 m, 153.3 m RMS, max velocity: 0.00352 m/year, 0.0064 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00023 *** Timestep 24: t=480000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.61359e+17, v = 2.61359e+17, p = 8.73802e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.105546, norm of the rhs: 2.75853e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.59201e-06, norm of the rhs: 2.2456e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.000991/0.0386/1.723e+08 // -0.04394/1.432/4.587e+10 // -0.05446/1.035/2.053e+10 // -0.055/1.055/1.99e+10 // -0.03495/1.054/6.196e+11 Heat flux densities for boundaries: -0.00941 W/m, -0.01534 W/m, -0.001212 W/m, 0.03518 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -5270 W, -8588 W, -1212 W, 3.518e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -7.243e+04 kg/yr, -1.383e+05 kg/yr, 1.714e+05 kg/yr, -2.619e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.569e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.6 K, 1630 K, 1893 K Topography min/max: -554.5 m, 207.1 m RMS, max velocity: 0.00143 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00024 *** Timestep 25: t=500000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61509e+17, v = 2.61509e+17, p = 8.74981e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0675114, norm of the rhs: 1.76548e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.9922e-06, norm of the rhs: 2.35154e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009845/0.03778/1.787e+08 // -0.04124/1.44/4.59e+10 // -0.05367/1.035/2.054e+10 // -0.05413/1.054/1.99e+10 // -0.03551/1.053/6.195e+11 Heat flux densities for boundaries: -0.03304 W/m, -0.01505 W/m, -0.001199 W/m, 0.03757 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -1.85e+04 W, -8430 W, -1199 W, 3.758e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -3.422e+05 kg/yr, -1.421e+05 kg/yr, 2.328e+05 kg/yr, -2.438e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.663e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.4 K, 1630 K, 1893 K Topography min/max: -559.4 m, 170.1 m RMS, max velocity: 0.00143 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00025 Termination requested by criterion: end time +----------------------------------------------+------------+------------+ | Total wallclock time elapsed since start | 180s | | | | | | | Section | no. calls | wall time | % of total | +----------------------------------+-----------+------------+------------+ | Assemble Stokes system | 27 | 2.27s | 1.3% | | Assemble Stokes system Picard | 136 | 10.9s | 6.1% | | Assemble Stokes system rhs | 109 | 7.06s | 3.9% | | Assemble composition system | 135 | 21.3s | 12% | | Assemble temperature system | 27 | 4.62s | 2.6% | | Build Stokes preconditioner | 136 | 13.1s | 7.3% | | Build composition preconditioner | 133 | 0.833s | 0.46% | | Build temperature preconditioner | 27 | 0.167s | 0% | | Create snapshot | 2 | 1.39s | 0.77% | | Initialization | 1 | 0.145s | 0% | | Mesh deformation | 27 | 1.13s | 0.63% | | Mesh deformation initialize | 2 | 0.0433s | 0% | | Postprocessing | 26 | 6.91s | 3.8% | | Refine mesh structure, part 1 | 1 | 0.0303s | 0% | | Refine mesh structure, part 2 | 1 | 0.0383s | 0% | | Setup dof systems | 2 | 0.132s | 0% | | Setup initial conditions | 2 | 0.236s | 0.13% | | Setup matrices | 27 | 6.02s | 3.4% | | Solve Stokes system | 136 | 99.9s | 56% | | Solve composition system | 133 | 1.25s | 0.69% | | Solve temperature system | 27 | 0.221s | 0.12% | +----------------------------------+-----------+------------+------------+ -- Total wallclock time elapsed including restarts: 180s ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.4.0&NewtonSolver=1&sha=&src=code ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- This is ASPECT, the Advanced Solver for Problems in Earth's ConvecTion. -- . version 2.4.0 -- . using deal.II 9.3.0 -- . with 32 bit indices and vectorization level 1 (128 bits) -- . using Trilinos 13.2.0 -- . using p4est 2.3.2 -- . running in OPTIMIZED mode -- . running with 32 MPI processes ----------------------------------------------------------------------------- ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.4.0&NewtonSolver=1&sha=&src=code ----------------------------------------------------------------------------- *** Resuming from snapshot! Number of active cells: 31,920 (on 4 levels) Number of degrees of freedom: 1,063,943 (257,874+32,447+128,937+128,937+128,937+128,937+128,937+128,937) Number of mesh deformation degrees of freedom: 64894 Solving mesh displacement system... 0 iterations. *** Timestep 20: t=400000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.64153e+17, v = 2.64153e+17, p = 8.75977e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.145326, norm of the rhs: 3.83884e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+39 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.09605e-05, norm of the rhs: 2.89525e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001259/0.04097/1.47e+08 // -0.04214/1.457/4.585e+10 // -0.06607/1.037/2.055e+10 // -0.06645/1.066/1.992e+10 // -0.03581/1.066/6.195e+11 Heat flux densities for boundaries: -0.3861 W/m, -0.3663 W/m, -0.138 W/m, 0.05963 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -2.162e+05 W, -2.051e+05 W, -1.38e+05 W, 5.963e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -4.228e+06 kg/yr, -4.008e+06 kg/yr, 4.278e+06 kg/yr, 7.203e+05 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -9.038e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 2393 K Topography min/max: -467.8 m, 176 m RMS, max velocity: 0.00243 m/year, 0.00548 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00020 *** Timestep 21: t=420000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.64901e+17, v = 2.64901e+17, p = 8.73638e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0221307, norm of the rhs: 5.86244e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.81484e-06, norm of the rhs: 2.07016e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001195/0.04046/1.533e+08 // -0.04278/1.451/4.586e+10 // -0.06603/1.036/2.054e+10 // -0.06639/1.066/1.991e+10 // -0.0356/1.066/6.196e+11 Heat flux densities for boundaries: 0.1337 W/m, 0.14 W/m, -0.0001638 W/m, 0.02601 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: 7.489e+04 W, 7.842e+04 W, -163.8 W, 2.601e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: 1.503e+06 kg/yr, 1.574e+06 kg/yr, -1.535e+06 kg/yr, -3.966e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -6.07e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.1 K, 1630 K, 2393 K Topography min/max: -494.1 m, 187 m RMS, max velocity: 0.00186 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00021 *** Timestep 22: t=440000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.64753e+17, v = 2.64753e+17, p = 8.74245e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0744296, norm of the rhs: 1.97055e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+33 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.73376e-06, norm of the rhs: 2.57705e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001138/0.04005/1.597e+08 // -0.0446/1.458/4.589e+10 // -0.06375/1.037/2.056e+10 // -0.06411/1.064/1.991e+10 // -0.03668/1.064/6.195e+11 Heat flux densities for boundaries: -0.257 W/m, -0.2755 W/m, -0.02966 W/m, 0.05333 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -1.439e+05 W, -1.543e+05 W, -2.966e+04 W, 5.334e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -2.81e+06 kg/yr, -3.015e+06 kg/yr, 2.972e+06 kg/yr, -2.223e+05 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.61e+07 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.9 K, 1630 K, 2393 K Topography min/max: -487.8 m, 151.6 m RMS, max velocity: 0.00193 m/year, 0.0054 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00022 *** Timestep 23: t=460000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63835e+17, v = 2.63835e+17, p = 8.73643e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0843359, norm of the rhs: 2.22508e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.1398e-06, norm of the rhs: 2.14756e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001052/0.03948/1.66e+08 // -0.04504/1.445/4.588e+10 // -0.06011/1.036/2.055e+10 // -0.06045/1.06/1.991e+10 // -0.03608/1.06/6.195e+11 Heat flux densities for boundaries: -0.5509 W/m, -0.574 W/m, -0.02614 W/m, 0.07162 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -3.085e+05 W, -3.215e+05 W, -2.614e+04 W, 7.163e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -6.045e+06 kg/yr, -6.302e+06 kg/yr, 6.372e+06 kg/yr, 2.438e+06 kg/yr, 1.483e+06 kg/yr, 1.483e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.14e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1630 K, 2393 K Topography min/max: -510.2 m, 153.3 m RMS, max velocity: 0.00352 m/year, 0.00641 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00023 *** Timestep 24: t=480000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.63144e+17, v = 2.63144e+17, p = 8.73285e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.103349, norm of the rhs: 2.71958e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.48055e-06, norm of the rhs: 2.23161e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009908/0.0386/1.723e+08 // -0.04394/1.432/4.587e+10 // -0.05553/1.035/2.053e+10 // -0.05578/1.056/1.99e+10 // -0.03502/1.055/6.196e+11 Heat flux densities for boundaries: -0.008704 W/m, -0.01473 W/m, -0.02244 W/m, 0.03513 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -4874 W, -8247 W, -2.244e+04 W, 3.513e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -6.458e+04 kg/yr, -1.315e+05 kg/yr, 1.643e+05 kg/yr, -2.626e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.587e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.6 K, 1630 K, 2393 K Topography min/max: -554.5 m, 207.3 m RMS, max velocity: 0.00143 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00024 *** Timestep 25: t=500000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.62955e+17, v = 2.62955e+17, p = 8.7436e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0674526, norm of the rhs: 1.7737e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.85676e-06, norm of the rhs: 2.32893e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009839/0.03779/1.787e+08 // -0.04124/1.44/4.59e+10 // -0.05375/1.035/2.054e+10 // -0.05415/1.054/1.99e+10 // -0.03561/1.053/6.195e+11 Heat flux densities for boundaries: -0.03209 W/m, -0.01402 W/m, -0.02143 W/m, 0.03751 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -1.797e+04 W, -7851 W, -2.143e+04 W, 3.752e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -3.319e+05 kg/yr, -1.308e+05 kg/yr, 2.206e+05 kg/yr, -2.447e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.643e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.4 K, 1630 K, 2393 K Topography min/max: -559.4 m, 170.1 m RMS, max velocity: 0.00144 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00025 *** Timestep 26: t=520000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.62523e+17, v = 2.62523e+17, p = 8.75393e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+5 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00360182, norm of the rhs: 9.45561e+14 Rebuilding Stokes preconditioner... Solving Stokes system... 0+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.86072e-06, norm of the rhs: 2.06362e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001007/0.03881/1.851e+08 // -0.03687/1.441/4.591e+10 // -0.05292/1.035/2.054e+10 // -0.0533/1.054/1.989e+10 // -0.03544/1.052/6.195e+11 Heat flux densities for boundaries: -0.5795 W/m, -0.5531 W/m, -0.02068 W/m, 0.07279 W/m, 0.4092 W/m, 0.4092 W/m Heat fluxes through boundary parts: -3.245e+05 W, -3.097e+05 W, -2.068e+04 W, 7.28e+04 W, 4.092e+04 W, 4.092e+04 W Mass fluxes through boundary parts: -6.368e+06 kg/yr, -6.075e+06 kg/yr, 6.364e+06 kg/yr, 2.527e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.48e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1630 K, 2393 K Topography min/max: -564.5 m, 134.8 m RMS, max velocity: 0.00343 m/year, 0.00644 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00026 *** Timestep 27: t=540000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61944e+17, v = 2.61944e+17, p = 8.73109e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.157163, norm of the rhs: 4.11679e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.12423e-05, norm of the rhs: 2.94485e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001002/0.04062/1.913e+08 // -0.03078/1.424/4.589e+10 // -0.05271/1.034/2.052e+10 // -0.05299/1.054/1.989e+10 // -0.0339/1.052/6.196e+11 Heat flux densities for boundaries: -0.2973 W/m, -0.2939 W/m, -0.01865 W/m, 0.05415 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -1.665e+05 W, -1.646e+05 W, -1.865e+04 W, 5.416e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -3.249e+06 kg/yr, -3.212e+06 kg/yr, 3.358e+06 kg/yr, -8861 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.75e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 270 K, 1631 K, 2393 K Topography min/max: -608.9 m, 200.1 m RMS, max velocity: 0.00199 m/year, 0.00542 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00027 *** Timestep 28: t=560000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.6191e+17, v = 2.6191e+17, p = 8.72081e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.00369982, norm of the rhs: 9.69021e+14 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.23219e-06, norm of the rhs: 1.89419e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009541/0.04248/1.976e+08 // -0.0244/1.422/4.59e+10 // -0.05212/1.033/2.052e+10 // -0.05247/1.053/1.989e+10 // -0.03361/1.051/6.196e+11 Heat flux densities for boundaries: 0.1006 W/m, 0.07233 W/m, -0.01759 W/m, 0.03014 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: 5.632e+04 W, 4.051e+04 W, -1.759e+04 W, 3.014e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: 1.136e+06 kg/yr, 8.211e+05 kg/yr, -9.984e+05 kg/yr, -3.469e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.773e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.6 K, 1631 K, 2393 K Topography min/max: -634.4 m, 211.2 m RMS, max velocity: 0.00178 m/year, 0.00537 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00028 *** Timestep 29: t=580000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61855e+17, v = 2.61855e+17, p = 8.74281e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0786672, norm of the rhs: 2.05994e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+34 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.9343e-06, norm of the rhs: 2.33949e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008839/0.04452/2.041e+08 // -0.02254/1.43/4.593e+10 // -0.05129/1.034/2.053e+10 // -0.05164/1.052/1.988e+10 // -0.03414/1.051/6.195e+11 Heat flux densities for boundaries: -0.3489 W/m, -0.3871 W/m, -0.01761 W/m, 0.06053 W/m, 0.4091 W/m, 0.4091 W/m Heat fluxes through boundary parts: -1.954e+05 W, -2.168e+05 W, -1.761e+04 W, 6.053e+04 W, 4.091e+04 W, 4.091e+04 W Mass fluxes through boundary parts: -3.825e+06 kg/yr, -4.25e+06 kg/yr, 4.103e+06 kg/yr, 7.312e+05 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.447e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.3 K, 1631 K, 2393 K Topography min/max: -632.1 m, 144.9 m RMS, max velocity: 0.00231 m/year, 0.00546 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00029 *** Snapshot created! *** Timestep 30: t=600000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.61384e+17, v = 2.61384e+17, p = 8.71998e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.121076, norm of the rhs: 3.16474e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.72819e-06, norm of the rhs: 2.28141e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008579/0.04658/2.103e+08 // -0.01952/1.418/4.592e+10 // -0.05092/1.033/2.051e+10 // -0.0511/1.052/1.988e+10 // -0.03304/1.05/6.196e+11 Heat flux densities for boundaries: -0.5418 W/m, -0.5506 W/m, -0.01669 W/m, 0.07116 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -3.034e+05 W, -3.083e+05 W, -1.669e+04 W, 7.117e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -5.946e+06 kg/yr, -6.044e+06 kg/yr, 6.188e+06 kg/yr, 2.289e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.441e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.8 K, 1631 K, 2393 K Topography min/max: -662.3 m, 167.6 m RMS, max velocity: 0.0033 m/year, 0.00609 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00030 *** Timestep 31: t=620000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 13 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.61144e+17, v = 2.61144e+17, p = 8.70233e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.11214, norm of the rhs: 2.92845e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.93049e-06, norm of the rhs: 2.071e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008693/0.04807/2.165e+08 // -0.01688/1.405/4.591e+10 // -0.05065/1.031/2.05e+10 // -0.05071/1.052/1.987e+10 // -0.03196/1.049/6.196e+11 Heat flux densities for boundaries: -0.03299 W/m, -0.006992 W/m, -0.01544 W/m, 0.03656 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -1.848e+04 W, -3915 W, -1.544e+04 W, 3.656e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -3.339e+05 kg/yr, -4.527e+04 kg/yr, 2.589e+05 kg/yr, -2.55e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.638e+08 Pa, 9.983e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1631 K, 2393 K Topography min/max: -706 m, 211.6 m RMS, max velocity: 0.00147 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00031 *** Timestep 32: t=640000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.61129e+17, v = 2.61129e+17, p = 8.70705e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0766705, norm of the rhs: 2.00209e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+37 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.48365e-06, norm of the rhs: 1.9542e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008486/0.04921/2.229e+08 // -0.01616/1.413/4.594e+10 // -0.04986/1.031/2.051e+10 // -0.04983/1.051/1.987e+10 // -0.03234/1.048/6.196e+11 Heat flux densities for boundaries: -0.07687 W/m, -0.04298 W/m, -0.01534 W/m, 0.04019 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -4.305e+04 W, -2.407e+04 W, -1.534e+04 W, 4.02e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -8.245e+05 kg/yr, -4.478e+05 kg/yr, 6.569e+05 kg/yr, -2.129e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.447e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.3 K, 1631 K, 2393 K Topography min/max: -712.1 m, 169.8 m RMS, max velocity: 0.00144 m/year, 0.00535 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00032 *** Timestep 33: t=660000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.60685e+17, v = 2.60685e+17, p = 8.70979e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0136534, norm of the rhs: 3.55926e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+27 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.46645e-06, norm of the rhs: 1.94639e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009135/0.05025/2.292e+08 // -0.01567/1.413/4.595e+10 // -0.04905/1.031/2.051e+10 // -0.04892/1.05/1.987e+10 // -0.03191/1.048/6.196e+11 Heat flux densities for boundaries: -0.5682 W/m, -0.5643 W/m, -0.01524 W/m, 0.07285 W/m, 0.409 W/m, 0.409 W/m Heat fluxes through boundary parts: -3.182e+05 W, -3.16e+05 W, -1.524e+04 W, 7.286e+04 W, 4.09e+04 W, 4.09e+04 W Mass fluxes through boundary parts: -6.238e+06 kg/yr, -6.195e+06 kg/yr, 6.407e+06 kg/yr, 2.482e+06 kg/yr, 1.482e+06 kg/yr, 1.482e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.398e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1631 K, 2393 K Topography min/max: -720.7 m, 152.5 m RMS, max velocity: 0.00346 m/year, 0.00629 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00033 *** Timestep 34: t=680000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.60173e+17, v = 2.60173e+17, p = 8.6875e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.142976, norm of the rhs: 3.71985e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+32 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.04493e-05, norm of the rhs: 2.71861e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008523/0.05096/2.353e+08 // -0.01504/1.398/4.594e+10 // -0.04856/1.029/2.049e+10 // -0.04821/1.05/1.986e+10 // -0.03035/1.047/6.196e+11 Heat flux densities for boundaries: -0.2589 W/m, -0.2931 W/m, -0.01419 W/m, 0.0529 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -1.45e+05 W, -1.642e+05 W, -1.419e+04 W, 5.291e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -2.82e+06 kg/yr, -3.2e+06 kg/yr, 3.18e+06 kg/yr, -2.336e+05 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.258e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.6 K, 1631 K, 2393 K Topography min/max: -765.9 m, 195.5 m RMS, max velocity: 0.00197 m/year, 0.00542 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00034 *** Timestep 35: t=700000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.60173e+17, v = 2.60173e+17, p = 8.67758e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.000497306, norm of the rhs: 1.29386e+14 Rebuilding Stokes preconditioner... Solving Stokes system... 0+23 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.74217e-06, norm of the rhs: 1.75413e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.000839/0.05155/2.416e+08 // -0.01414/1.397/4.595e+10 // -0.04784/1.029/2.049e+10 // -0.04733/1.049/1.986e+10 // -0.02997/1.046/6.196e+11 Heat flux densities for boundaries: 0.1306 W/m, 0.09639 W/m, -0.01367 W/m, 0.02817 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: 7.314e+04 W, 5.398e+04 W, -1.367e+04 W, 2.817e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: 1.47e+06 kg/yr, 1.09e+06 kg/yr, -1.268e+06 kg/yr, -3.754e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.226e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.7 K, 1631 K, 2393 K Topography min/max: -790.8 m, 201.7 m RMS, max velocity: 0.0018 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00035 *** Timestep 36: t=720000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.60098e+17, v = 2.60098e+17, p = 8.69664e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0757828, norm of the rhs: 1.9711e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.97485e-06, norm of the rhs: 2.07424e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008641/0.05218/2.481e+08 // -0.01425/1.402/4.598e+10 // -0.04704/1.029/2.051e+10 // -0.04651/1.049/1.985e+10 // -0.03048/1.045/6.196e+11 Heat flux densities for boundaries: -0.3773 W/m, -0.373 W/m, -0.01395 W/m, 0.06081 W/m, 0.4089 W/m, 0.4089 W/m Heat fluxes through boundary parts: -2.113e+05 W, -2.089e+05 W, -1.395e+04 W, 6.082e+04 W, 4.089e+04 W, 4.089e+04 W Mass fluxes through boundary parts: -4.137e+06 kg/yr, -4.089e+06 kg/yr, 4.218e+06 kg/yr, 7.609e+05 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.439e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.6 K, 1631 K, 2393 K Topography min/max: -787.4 m, 165.1 m RMS, max velocity: 0.0024 m/year, 0.00547 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00036 *** Timestep 37: t=740000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.59588e+17, v = 2.59588e+17, p = 8.68432e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.116531, norm of the rhs: 3.02502e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.67637e-06, norm of the rhs: 1.99269e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008979/0.05261/2.542e+08 // -0.01438/1.391/4.597e+10 // -0.04626/1.028/2.049e+10 // -0.04546/1.048/1.985e+10 // -0.02936/1.044/6.196e+11 Heat flux densities for boundaries: -0.5154 W/m, -0.4761 W/m, -0.01337 W/m, 0.06747 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -2.886e+05 W, -2.666e+05 W, -1.337e+04 W, 6.749e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -5.655e+06 kg/yr, -5.218e+06 kg/yr, 5.641e+06 kg/yr, 1.805e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.522e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.4 K, 1631 K, 2393 K Topography min/max: -818.9 m, 169.6 m RMS, max velocity: 0.00306 m/year, 0.00562 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00037 *** Timestep 38: t=760000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 13 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.59478e+17, v = 2.59478e+17, p = 8.6774e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0891317, norm of the rhs: 2.31277e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.41449e-06, norm of the rhs: 1.9239e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009247/0.05421/2.604e+08 // -0.01452/1.381/4.596e+10 // -0.04561/1.027/2.047e+10 // -0.04467/1.047/1.985e+10 // -0.02841/1.043/6.196e+11 Heat flux densities for boundaries: 0.01093 W/m, 0.03731 W/m, -0.01262 W/m, 0.03367 W/m, 0.4087 W/m, 0.4087 W/m Heat fluxes through boundary parts: 6119 W, 2.089e+04 W, -1.262e+04 W, 3.368e+04 W, 4.087e+04 W, 4.087e+04 W Mass fluxes through boundary parts: 1.493e+05 kg/yr, 4.42e+05 kg/yr, -2.456e+05 kg/yr, -2.947e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.317e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1631 K, 2393 K Topography min/max: -859.5 m, 222.4 m RMS, max velocity: 0.00155 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00038 *** Timestep 39: t=780000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.59593e+17, v = 2.59593e+17, p = 8.69201e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0851428, norm of the rhs: 2.21025e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+29 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.97354e-06, norm of the rhs: 2.06988e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009678/0.05626/2.668e+08 // -0.01438/1.389/4.599e+10 // -0.04485/1.027/2.049e+10 // -0.04375/1.047/1.984e+10 // -0.0288/1.042/6.196e+11 Heat flux densities for boundaries: -0.1363 W/m, -0.1566 W/m, -0.01281 W/m, 0.04623 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -7.631e+04 W, -8.768e+04 W, -1.281e+04 W, 4.624e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -1.479e+06 kg/yr, -1.705e+06 kg/yr, 1.623e+06 kg/yr, -1.323e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.08e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.1 K, 1631 K, 2393 K Topography min/max: -862.9 m, 175.1 m RMS, max velocity: 0.00148 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00039 *** Snapshot created! *** Timestep 40: t=800000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.59221e+17, v = 2.59221e+17, p = 8.69023e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0478832, norm of the rhs: 1.24123e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.00103e-06, norm of the rhs: 2.07403e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.001018/0.05742/2.73e+08 // -0.01429/1.386/4.599e+10 // -0.04397/1.026/2.048e+10 // -0.04276/1.046/1.984e+10 // -0.02806/1.041/6.196e+11 Heat flux densities for boundaries: -0.5466 W/m, -0.6032 W/m, -0.01274 W/m, 0.07354 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -3.061e+05 W, -3.378e+05 W, -1.274e+04 W, 7.356e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -5.998e+06 kg/yr, -6.628e+06 kg/yr, 6.503e+06 kg/yr, 2.564e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.526e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.1 K, 1631 K, 2393 K Topography min/max: -878 m, 150.5 m RMS, max velocity: 0.00345 m/year, 0.00655 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00040 *** Timestep 41: t=820000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.58959e+17, v = 2.58959e+17, p = 8.67015e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+11 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.152124, norm of the rhs: 3.9394e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+30 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 1.00614e-05, norm of the rhs: 2.6055e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009673/0.05815/2.791e+08 // -0.01399/1.374/4.598e+10 // -0.04311/1.025/2.046e+10 // -0.04188/1.045/1.983e+10 // -0.02658/1.04/6.196e+11 Heat flux densities for boundaries: -0.1863 W/m, -0.2245 W/m, -0.01201 W/m, 0.04874 W/m, 0.4086 W/m, 0.4086 W/m Heat fluxes through boundary parts: -1.043e+05 W, -1.257e+05 W, -1.201e+04 W, 4.875e+04 W, 4.086e+04 W, 4.086e+04 W Mass fluxes through boundary parts: -2.021e+06 kg/yr, -2.446e+06 kg/yr, 2.365e+06 kg/yr, -8.619e+05 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.43e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.7 K, 1631 K, 2393 K Topography min/max: -924.2 m, 208.3 m RMS, max velocity: 0.00165 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00041 *** Timestep 42: t=840000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.59113e+17, v = 2.59113e+17, p = 8.66809e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0211685, norm of the rhs: 5.48505e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.13019e-06, norm of the rhs: 1.84753e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.000969/0.05967/2.854e+08 // -0.01393/1.375/4.6e+10 // -0.04228/1.024/2.046e+10 // -0.04097/1.044/1.983e+10 // -0.02641/1.039/6.196e+11 Heat flux densities for boundaries: 0.02408 W/m, 0.04527 W/m, -0.01186 W/m, 0.03409 W/m, 0.4086 W/m, 0.4087 W/m Heat fluxes through boundary parts: 1.348e+04 W, 2.535e+04 W, -1.186e+04 W, 3.41e+04 W, 4.086e+04 W, 4.087e+04 W Mass fluxes through boundary parts: 2.904e+05 kg/yr, 5.26e+05 kg/yr, -3.952e+05 kg/yr, -3.012e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.475e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.4 K, 1631 K, 2393 K Topography min/max: -944.3 m, 199.5 m RMS, max velocity: 0.0016 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00042 *** Timestep 43: t=860000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.59042e+17, v = 2.59042e+17, p = 8.68083e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+8 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0356554, norm of the rhs: 9.23627e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+31 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.77486e-06, norm of the rhs: 2.01402e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009976/0.06106/2.918e+08 // -0.01395/1.377/4.602e+10 // -0.04155/1.024/2.047e+10 // -0.04021/1.044/1.983e+10 // -0.02666/1.039/6.196e+11 Heat flux densities for boundaries: -0.4761 W/m, -0.4028 W/m, -0.01208 W/m, 0.06523 W/m, 0.4088 W/m, 0.4088 W/m Heat fluxes through boundary parts: -2.666e+05 W, -2.256e+05 W, -1.208e+04 W, 6.524e+04 W, 4.088e+04 W, 4.088e+04 W Mass fluxes through boundary parts: -5.23e+06 kg/yr, -4.414e+06 kg/yr, 4.956e+06 kg/yr, 1.341e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.874e+08 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.6 K, 1631 K, 2393 K Topography min/max: -946.9 m, 165.9 m RMS, max velocity: 0.00272 m/year, 0.00553 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00043 *** Timestep 44: t=880000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.58625e+17, v = 2.58625e+17, p = 8.66368e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.123558, norm of the rhs: 3.19552e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 8.07234e-06, norm of the rhs: 2.08771e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009639/0.06207/2.978e+08 // -0.01389/1.368/4.601e+10 // -0.04054/1.023/2.045e+10 // -0.03906/1.043/1.982e+10 // -0.02549/1.037/6.196e+11 Heat flux densities for boundaries: -0.4587 W/m, -0.4159 W/m, -0.01163 W/m, 0.064 W/m, 0.4086 W/m, 0.4086 W/m Heat fluxes through boundary parts: -2.569e+05 W, -2.329e+05 W, -1.163e+04 W, 6.402e+04 W, 4.086e+04 W, 4.086e+04 W Mass fluxes through boundary parts: -5.029e+06 kg/yr, -4.553e+06 kg/yr, 4.994e+06 kg/yr, 1.255e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.624e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1631 K, 2393 K Topography min/max: -983.3 m, 182.7 m RMS, max velocity: 0.00272 m/year, 0.00553 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00044 *** Timestep 45: t=900000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.58677e+17, v = 2.58677e+17, p = 8.6517e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+6 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0721033, norm of the rhs: 1.86515e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.73533e-06, norm of the rhs: 1.74228e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009156/0.06364/3.04e+08 // -0.01393/1.362/4.601e+10 // -0.03981/1.022/2.044e+10 // -0.03809/1.042/1.982e+10 // -0.02479/1.036/6.197e+11 Heat flux densities for boundaries: 0.05758 W/m, 0.01746 W/m, -0.01115 W/m, 0.03324 W/m, 0.4085 W/m, 0.4085 W/m Heat fluxes through boundary parts: 3.224e+04 W, 9775 W, -1.115e+04 W, 3.325e+04 W, 4.085e+04 W, 4.085e+04 W Mass fluxes through boundary parts: 6.661e+05 kg/yr, 2.2e+05 kg/yr, -3.934e+05 kg/yr, -3.073e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.295e+08 Pa, 9.984e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269 K, 1631 K, 2393 K Topography min/max: -1020 m, 231.1 m RMS, max velocity: 0.0016 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00045 *** Timestep 46: t=920000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.58827e+17, v = 2.58827e+17, p = 8.66644e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0768882, norm of the rhs: 1.99007e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.7246e-06, norm of the rhs: 1.74051e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.000846/0.0658/3.104e+08 // -0.0139/1.368/4.603e+10 // -0.03921/1.023/2.046e+10 // -0.03747/1.041/1.981e+10 // -0.02527/1.036/6.196e+11 Heat flux densities for boundaries: -0.1645 W/m, -0.2566 W/m, -0.01139 W/m, 0.05051 W/m, 0.4086 W/m, 0.4086 W/m Heat fluxes through boundary parts: -9.212e+04 W, -1.437e+05 W, -1.139e+04 W, 5.053e+04 W, 4.086e+04 W, 4.086e+04 W Mass fluxes through boundary parts: -1.785e+06 kg/yr, -2.81e+06 kg/yr, 2.366e+06 kg/yr, -7.537e+05 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.542e+08 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 267.9 K, 1631 K, 2393 K Topography min/max: -1023 m, 169.5 m RMS, max velocity: 0.00171 m/year, 0.00538 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00046 *** Timestep 47: t=940000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.58375e+17, v = 2.58375e+17, p = 8.66008e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+9 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0746339, norm of the rhs: 1.92835e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+25 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.71653e-06, norm of the rhs: 1.73538e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008485/0.06651/3.165e+08 // -0.01378/1.363/4.603e+10 // -0.03913/1.022/2.045e+10 // -0.03649/1.04/1.981e+10 // -0.02445/1.035/6.196e+11 Heat flux densities for boundaries: -0.5288 W/m, -0.5746 W/m, -0.0113 W/m, 0.07195 W/m, 0.4086 W/m, 0.4086 W/m Heat fluxes through boundary parts: -2.961e+05 W, -3.218e+05 W, -1.13e+04 W, 7.197e+04 W, 4.086e+04 W, 4.086e+04 W Mass fluxes through boundary parts: -5.8e+06 kg/yr, -6.31e+06 kg/yr, 6.267e+06 kg/yr, 2.324e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -1.531e+08 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.3 K, 1631 K, 2393 K Topography min/max: -1043 m, 173.9 m RMS, max velocity: 0.00334 m/year, 0.00621 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00047 *** Timestep 48: t=960000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 15 iterations. Solving plastic_strain system ... 14 iterations. Solving crust_upper system ... 14 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 12 iterations. Initial Newton Stokes residual = 2.58133e+17, v = 2.58133e+17, p = 8.64368e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.152598, norm of the rhs: 3.93906e+16 Rebuilding Stokes preconditioner... Solving Stokes system... 0+28 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 9.2355e-06, norm of the rhs: 2.38398e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.000906/0.06671/3.225e+08 // -0.01376/1.353/4.602e+10 // -0.03836/1.021/2.043e+10 // -0.03538/1.04/1.981e+10 // -0.02318/1.034/6.197e+11 Heat flux densities for boundaries: -0.1504 W/m, -0.104 W/m, -0.01073 W/m, 0.04369 W/m, 0.4084 W/m, 0.4085 W/m Heat fluxes through boundary parts: -8.422e+04 W, -5.822e+04 W, -1.073e+04 W, 4.371e+04 W, 4.084e+04 W, 4.085e+04 W Mass fluxes through boundary parts: -1.629e+06 kg/yr, -1.114e+06 kg/yr, 1.489e+06 kg/yr, -1.585e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -9.983e+07 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 269.2 K, 1631 K, 2393 K Topography min/max: -1088 m, 217.5 m RMS, max velocity: 0.00148 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00048 *** Timestep 49: t=980000 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 9 iterations. Solving noninitial_plastic_strain system ... 16 iterations. Solving plastic_strain system ... 15 iterations. Solving crust_upper system ... 15 iterations. Solving crust_lower system ... 15 iterations. Solving mantle_lithosphere system ... 13 iterations. Initial Newton Stokes residual = 2.58322e+17, v = 2.58322e+17, p = 8.65005e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+10 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0294536, norm of the rhs: 7.6085e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+27 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 6.84599e-06, norm of the rhs: 1.76847e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0008893/0.06833/3.289e+08 // -0.01376/1.355/4.604e+10 // -0.03833/1.021/2.044e+10 // -0.03453/1.039/1.98e+10 // -0.02323/1.033/6.197e+11 Heat flux densities for boundaries: -0.07113 W/m, 0.01541 W/m, -0.01076 W/m, 0.03848 W/m, 0.4085 W/m, 0.4085 W/m Heat fluxes through boundary parts: -3.983e+04 W, 8632 W, -1.076e+04 W, 3.85e+04 W, 4.085e+04 W, 4.085e+04 W Mass fluxes through boundary parts: -7.628e+05 kg/yr, 1.988e+05 kg/yr, 3.144e+05 kg/yr, -2.441e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -8.253e+07 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.2 K, 1631 K, 2393 K Topography min/max: -1103 m, 191.7 m RMS, max velocity: 0.00153 m/year, 0.00536 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00049 *** Snapshot created! *** Timestep 50: t=1e+06 years, dt=20000 years Solving mesh surface diffusion Solving mesh displacement system... 8 iterations. Solving temperature system... 8 iterations. Solving noninitial_plastic_strain system ... 14 iterations. Solving plastic_strain system ... 13 iterations. Solving crust_upper system ... 13 iterations. Solving crust_lower system ... 14 iterations. Solving mantle_lithosphere system ... 11 iterations. Initial Newton Stokes residual = 2.58112e+17, v = 2.58112e+17, p = 8.66287e+09 Rebuilding Stokes preconditioner... Solving Stokes system... 0+7 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 1: 0.0114267, norm of the rhs: 2.94936e+15 Rebuilding Stokes preconditioner... Solving Stokes system... 0+26 iterations. Relative nonlinear residual (total Newton system) after nonlinear iteration 2: 7.27285e-06, norm of the rhs: 1.87721e+12, newton_derivative_scaling_factor: 0 Postprocessing: Compositions min/max/mass: -0.0009352/0.06993/3.351e+08 // -0.01379/1.357/4.606e+10 // -0.03864/1.02/2.044e+10 // -0.0339/1.04/1.98e+10 // -0.02299/1.032/6.197e+11 Heat flux densities for boundaries: -0.4952 W/m, -0.4595 W/m, -0.01091 W/m, 0.06767 W/m, 0.4086 W/m, 0.4086 W/m Heat fluxes through boundary parts: -2.773e+05 W, -2.573e+05 W, -1.091e+04 W, 6.769e+04 W, 4.086e+04 W, 4.086e+04 W Mass fluxes through boundary parts: -5.437e+06 kg/yr, -5.041e+06 kg/yr, 5.395e+06 kg/yr, 1.678e+06 kg/yr, 1.481e+06 kg/yr, 1.481e+06 kg/yr Total system matrix memory consumption: 209.11 MB. Total system matrix nnz: 15,404,506 system matrix nnz by block: 8,138,934 1,592,771 32 32 32 32 32 32 1,592,771 236 32 32 32 32 32 32 32 32 2,022,221 32 32 32 32 32 32 32 32 2,055,845 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Total system preconditioner matrix memory consumption: 84.69 MB. Total system preconditioner matrix nnz: 4,360,397 system preconditioner matrix nnz by block: 4,070,562 32 32 32 32 32 32 32 32 288,043 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 32 32 32 32 32 32 32 32 0 Pressure min/avg/max: -9.642e+07 Pa, 9.985e+09 Pa, 2.008e+10 Pa Temperature min/avg/max: 268.8 K, 1631 K, 2393 K Topography min/max: -1111 m, 161.3 m RMS, max velocity: 0.0029 m/year, 0.00556 m/year Writing graphical output: output-TR_check_BC34/solution/solution-00050 Termination requested by criterion: end time +----------------------------------------------+------------+------------+ | Total wallclock time elapsed since start | 179s | | | | | | | Section | no. calls | wall time | % of total | +----------------------------------+-----------+------------+------------+ | Assemble Stokes system | 31 | 2.65s | 1.5% | | Assemble Stokes system Picard | 62 | 5.36s | 3% | | Assemble Stokes system rhs | 31 | 2.17s | 1.2% | | Assemble composition system | 155 | 25.2s | 14% | | Assemble temperature system | 31 | 5.36s | 3% | | Build Stokes preconditioner | 62 | 6.33s | 3.5% | | Build composition preconditioner | 155 | 0.975s | 0.54% | | Build temperature preconditioner | 31 | 0.191s | 0.11% | | Create snapshot | 3 | 1.95s | 1.1% | | Initialization | 1 | 0.124s | 0% | | Mesh deformation | 31 | 1.32s | 0.74% | | Mesh deformation initialize | 1 | 0.0295s | 0% | | Postprocessing | 31 | 8.32s | 4.6% | | Setup dof systems | 1 | 0.0879s | 0% | | Setup matrices | 31 | 6.96s | 3.9% | | Solve Stokes system | 62 | 108s | 60% | | Solve composition system | 155 | 1.46s | 0.81% | | Solve temperature system | 31 | 0.258s | 0.14% | +----------------------------------+-----------+------------+------------+ -- Total wallclock time elapsed including restarts: 324s ----------------------------------------------------------------------------- -- For information on how to cite ASPECT, see: -- https://aspect.geodynamics.org/citing.html?ver=2.4.0&NewtonSolver=1&sha=&src=code -----------------------------------------------------------------------------