
Number of active cells: 40 (on 1 levels)
Number of degrees of freedom: 1,985 (486+202+486+82+243+243+243)

*** Timestep 0:  t=0 seconds
   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.3822e-16, 3.39362e-17, 2.35518e-17, 1
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1

   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.3822e-16, 3.39362e-17, 2.35518e-17, 1.07096e-08
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.07096e-08


   Postprocessing:
     RMS, max velocity:        2e-11 m/s, 2e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00000

*** Timestep 1:  t=6.25e+14 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 18 iterations.
   Solving peridotite system ... 18 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.29988e-16, 0.0977234, 0.132953, 0.163961
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.163961

   Solving temperature system... 0 iterations.
   Solving porosity system ... 18 iterations.
   Solving peridotite system ... 18 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.59084e-16, 0.00268964, 0.00235709, 0.413162
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.413162

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.20809e-16, 2.78069e-11, 2.96114e-12, 2.36272e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.36272e-08


   Postprocessing:
     RMS, max velocity:        2.31e-11 m/s, 2.66e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00001

*** Timestep 2:  t=1.09488e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 15 iterations.
   Solving peridotite system ... 15 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.4678e-16, 0.0474001, 0.063079, 9.56028
      Relative nonlinear residual (total system) after nonlinear iteration 1: 9.56028

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 13 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.3006e-16, 0.000785074, 0.000306137, 0.0477706
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0477706

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.63044e-16, 1.03433e-11, 1.72909e-12, 2.52309e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.52309e-08


   Postprocessing:
     RMS, max velocity:        2.28e-11 m/s, 2.59e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00002

*** Timestep 3:  t=1.57707e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 15 iterations.
   Solving peridotite system ... 15 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.91402e-16, 0.0401071, 0.0463243, 2.67332
      Relative nonlinear residual (total system) after nonlinear iteration 1: 2.67332

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 13 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.65102e-16, 0.000460872, 0.000163195, 0.0136498
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0136498

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.35545e-16, 4.65946e-12, 1.14912e-12, 2.7471e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.7471e-08


   Postprocessing:
     RMS, max velocity:        2.27e-11 m/s, 2.55e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00003

*** Timestep 4:  t=2.06638e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 15 iterations.
   Solving peridotite system ... 15 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 4+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.56017e-16, 0.0346474, 0.0413239, 1.31839
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.31839

   Solving temperature system... 0 iterations.
   Solving porosity system ... 14 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.70792e-16, 0.000457554, 0.000142169, 0.0208133
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0208133

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.8359e-16, 4.50542e-12, 1.07411e-12, 6.14779e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 6.14779e-09


   Postprocessing:
     RMS, max velocity:        2.25e-11 m/s, 2.52e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00004

*** Timestep 5:  t=2.56171e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 15 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.35467e-16, 0.0319796, 0.0391285, 1.38118
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.38118

   Solving temperature system... 0 iterations.
   Solving porosity system ... 14 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.49273e-16, 0.000417912, 0.000120768, 0.0269993
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0269993

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.35481e-16, 3.35233e-12, 7.3181e-13, 8.26523e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 8.26523e-09


   Postprocessing:
     RMS, max velocity:        2.24e-11 m/s, 2.5e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00005

*** Timestep 6:  t=3.06202e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 15 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.86455e-16, 0.0314643, 0.0365334, 1.23056
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1.23056

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 13 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.19957e-16, 0.000362142, 0.000106821, 0.0274715
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0274715

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.15401e-16, 3.15127e-12, 3.85172e-13, 3.72203e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.72203e-08


   Postprocessing:
     RMS, max velocity:        2.23e-11 m/s, 2.48e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00006

*** Timestep 7:  t=3.56642e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 15 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.27258e-16, 0.0316403, 0.0352049, 0.924719
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.924719

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.15981e-16, 0.000305358, 9.43437e-05, 0.0225622
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0225622

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.12341e-16, 2.47389e-12, 1.06044e-12, 2.46442e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.46442e-09


   Postprocessing:
     RMS, max velocity:        2.22e-11 m/s, 2.46e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00007

*** Timestep 8:  t=4.07414e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 15 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.67164e-16, 0.0320587, 0.0349608, 0.635931
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.635931

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.01781e-16, 0.00025389, 8.07423e-05, 0.0168946
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0168946

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.98285e-16, 1.66341e-12, 4.96437e-13, 3.29916e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.29916e-09


   Postprocessing:
     RMS, max velocity:        2.22e-11 m/s, 2.45e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00008

*** Timestep 9:  t=4.58453e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.33962e-16, 0.0325536, 0.0342598, 0.415632
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.415632

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.42763e-16, 0.000203212, 6.9242e-05, 0.0115908
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0115908

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.45933e-16, 3.94072e-12, 1.45881e-12, 3.74458e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.74458e-09


   Postprocessing:
     RMS, max velocity:        2.21e-11 m/s, 2.44e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00009

*** Timestep 10:  t=5.09708e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.28308e-16, 0.0329901, 0.0337229, 0.2627
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.2627

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.0708e-16, 0.000162166, 5.87107e-05, 0.00756804
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00756804

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.20327e-16, 5.25168e-12, 1.26606e-12, 2.25009e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.25009e-08


   Postprocessing:
     RMS, max velocity:        2.21e-11 m/s, 2.43e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00010

*** Timestep 11:  t=5.61134e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.08478e-16, 0.0333404, 0.0334707, 0.160061
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.160061

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.40769e-16, 0.000129884, 4.85666e-05, 0.00492007
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00492007

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.18953e-16, 5.9003e-12, 2.30253e-12, 6.12643e-10
      Relative nonlinear residual (total system) after nonlinear iteration 3: 6.12643e-10


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.42e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00011

*** Timestep 12:  t=6.12696e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.472e-16, 0.0336004, 0.0332104, 0.093983
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.093983

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.13361e-16, 0.000103611, 3.95398e-05, 0.00326266
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00326266

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.04643e-16, 4.54033e-12, 1.81392e-12, 6.00069e-10
      Relative nonlinear residual (total system) after nonlinear iteration 3: 6.00069e-10


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.42e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00012

*** Timestep 13:  t=6.64367e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.75085e-16, 0.0337759, 0.0329617, 0.0535587
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0535587

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.00555e-16, 8.21064e-05, 3.16991e-05, 0.00225746
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00225746

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.02215e-16, 2.79124e-12, 1.15737e-12, 3.76081e-10
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.76081e-10


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.42e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00013

*** Timestep 14:  t=7.16123e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.17208e-16, 0.0338852, 0.0327774, 0.0305636
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0338852

   Solving temperature system... 0 iterations.
   Solving porosity system ... 13 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.95679e-16, 6.46137e-05, 2.50868e-05, 0.00165406
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00165406

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.03986e-16, 4.91209e-12, 1.79444e-12, 7.11963e-10
      Relative nonlinear residual (total system) after nonlinear iteration 3: 7.11963e-10


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.41e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00014

*** Timestep 15:  t=7.67947e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.45029e-16, 0.0339495, 0.032615, 0.0182209
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0339495

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 12 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.11854e-16, 5.04459e-05, 1.96708e-05, 0.00124588
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00124588

   Solving temperature system... 0 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.28739e-16, 3.39095e-11, 1.75566e-11, 1.00561e-07
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.00561e-07


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.41e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00015

*** Timestep 16:  t=8.19825e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.57174e-16, 0.0339875, 0.0324875, 0.0121604
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0339875

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.2178e-16, 3.9129e-05, 1.53154e-05, 0.000946562
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000946562

   Solving temperature system... 0 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.66279e-16, 2.82882e-11, 1.56251e-11, 8.04019e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 8.04019e-08


   Postprocessing:
     RMS, max velocity:        2.2e-11 m/s, 2.41e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00016

*** Timestep 17:  t=8.71744e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.01267e-16, 0.0340118, 0.0323611, 0.00974862
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340118

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.02684e-16, 3.01786e-05, 1.19097e-05, 0.000721586
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000721586

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.31933e-16, 2.3272e-11, 1.33806e-11, 9.82268e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 9.82268e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.41e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00017

*** Timestep 18:  t=9.23696e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.36018e-16, 0.0340293, 0.0322284, 0.00825934
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340293

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.91819e-16, 2.31964e-05, 9.28052e-06, 0.000552073
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000552073

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.23602e-16, 1.95351e-11, 1.17996e-11, 4.74088e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 4.74088e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00018

*** Timestep 19:  t=9.75673e+15 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.0052e-16, 0.0340433, 0.0321254, 0.00693213
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340433

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.25771e-16, 1.77944e-05, 7.2535e-06, 0.000424806
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000424806

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.24894e-16, 1.64858e-11, 1.0265e-11, 3.70615e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.70615e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00019

*** Timestep 20:  t=1.02767e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.47071e-16, 0.0340549, 0.0320306, 0.00572705
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340549

   Solving temperature system... 0 iterations.
   Solving porosity system ... 12 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.37258e-16, 1.36398e-05, 5.68528e-06, 0.000329395
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000329395

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.21597e-16, 1.35759e-11, 8.67777e-12, 2.91516e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.91516e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00020

*** Timestep 21:  t=1.07969e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.6201e-16, 0.0340644, 0.0319427, 0.00456778
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340644

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 11 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.50563e-16, 1.04561e-05, 4.45964e-06, 0.000257171
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000257171

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.88094e-16, 1.13097e-11, 7.39509e-12, 2.28914e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.28914e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00021

*** Timestep 22:  t=1.13171e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.28045e-16, 0.0340723, 0.0318661, 0.00356017
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340723

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 10 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.48666e-16, 8.01853e-06, 3.50006e-06, 0.000201691
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000201691

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.58707e-16, 9.25863e-12, 6.21023e-12, 1.80625e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.80625e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00022

*** Timestep 23:  t=1.18375e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.07311e-16, 0.0340786, 0.0317826, 0.00277297
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340786

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 10 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.13121e-16, 6.15518e-06, 2.748e-06, 0.000158486
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000158486

   Solving temperature system... 0 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.23049e-16, 7.46777e-12, 5.09945e-12, 1.59057e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.59057e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00023

*** Timestep 24:  t=1.23579e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.56292e-16, 0.0340837, 0.0317038, 0.00217877
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340837

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 10 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.16709e-16, 4.73161e-06, 2.15719e-06, 0.000124557
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000124557

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.41259e-16, 5.63263e-12, 3.89417e-12, 2.99819e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.99819e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00024

*** Timestep 25:  t=1.28784e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.25918e-16, 0.0340876, 0.0316337, 0.00171835
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340876

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 10 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.38334e-16, 3.64578e-06, 1.69269e-06, 9.78431e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 9.78431e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.34033e-16, 3.10895e-12, 2.21516e-12, 6.03811e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 6.03811e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00025

*** Timestep 26:  t=1.3399e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.26758e-16, 0.0340907, 0.0315666, 0.00134802
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340907

   Solving temperature system... 0 iterations.
   Solving porosity system ... 11 iterations.
   Solving peridotite system ... 9 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.22444e-16, 2.80911e-06, 1.32771e-06, 7.69286e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 7.69286e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.36201e-16, 3.60402e-12, 2.73187e-12, 7.00585e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 7.00585e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00026

*** Timestep 27:  t=1.39196e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.24354e-16, 0.0340931, 0.0315079, 0.0010431
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340931

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 9 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.36965e-16, 2.1675e-06, 1.04005e-06, 6.05219e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 6.05219e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.4299e-16, 2.99291e-12, 2.19469e-12, 5.94377e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 5.94377e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00027

*** Timestep 28:  t=1.44402e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.18351e-16, 0.0340951, 0.0314513, 0.000794254
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340951

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 9 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.16913e-16, 1.67285e-06, 8.141e-07, 4.76445e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 4.76445e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.25411e-16, 2.53428e-12, 1.83491e-12, 9.81738e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 9.81738e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00028

*** Timestep 29:  t=1.49608e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.83383e-16, 0.0340966, 0.0313983, 0.00059951
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340966

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 9 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.8132e-16, 1.30158e-06, 6.36704e-07, 3.749e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 3.749e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.9878e-16, 2.0432e-12, 1.47094e-12, 3.38592e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.38592e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00029

*** Timestep 30:  t=1.54815e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.07182e-16, 0.0340978, 0.031352, 0.000450821
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340978

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 9 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.88992e-16, 9.98422e-07, 4.98684e-07, 2.95249e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 2.95249e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.02451e-16, 1.20247e-12, 8.15253e-13, 1.76365e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.76365e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00030

*** Timestep 31:  t=1.60021e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.72109e-16, 0.0340987, 0.0313082, 0.000341327
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340987

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 8 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.67437e-16, 7.71351e-07, 3.899e-07, 2.31687e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 2.31687e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.60469e-16, 1.30197e-12, 1.18288e-12, 2.00277e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.00277e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00031

*** Timestep 32:  t=1.65228e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.40512e-16, 0.0340994, 0.031269, 0.000264222
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340994

   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 8 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.3664e-16, 5.95608e-07, 3.04446e-07, 1.81144e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.81144e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.41436e-16, 1.02738e-12, 9.23712e-13, 1.57552e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.57552e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00032

*** Timestep 33:  t=1.70435e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.5132e-16, 0.0340999, 0.0312309, 0.000208774
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0340999

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 8 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.59297e-16, 4.59515e-07, 2.37537e-07, 1.41144e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.41144e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.84773e-16, 1.15239e-12, 7.1583e-13, 1.37519e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.37519e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00033

*** Timestep 34:  t=1.75642e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.95683e-16, 0.0341002, 0.0311931, 0.000166517
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341002

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 8 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.52946e-16, 3.54296e-07, 1.85221e-07, 1.09721e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.09721e-05

   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.12122e-16, 9.00761e-13, 5.51509e-13, 2.55688e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.55688e-09


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00034

*** Timestep 35:  t=1.80849e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.22008e-16, 0.0341006, 0.0311579, 0.000132828
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341006

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 7 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.16094e-16, 2.73229e-07, 1.44338e-07, 8.52245e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 8.52245e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00035

*** Timestep 36:  t=1.86056e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.55068e-16, 0.0341008, 0.0311228, 0.000105275
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341008

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 7 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.07204e-16, 2.10999e-07, 1.12401e-07, 6.62114e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 6.62114e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00036

*** Timestep 37:  t=1.91263e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.06877e-16, 0.034101, 0.0310894, 8.28176e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.034101

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 7 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.75545e-16, 1.63237e-07, 8.74826e-08, 5.1482e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.1482e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00037

*** Timestep 38:  t=1.9647e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.04605e-16, 0.0341012, 0.0310573, 6.48154e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341012

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 7 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.25916e-16, 1.26246e-07, 6.81012e-08, 4.00596e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 4.00596e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00038

*** Timestep 39:  t=2.01677e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.5064e-16, 0.0341013, 0.031025, 5.0587e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341013

   Solving temperature system... 0 iterations.
   Solving porosity system ... 9 iterations.
   Solving peridotite system ... 7 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.4571e-16, 9.70444e-08, 5.30919e-08, 3.12431e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 3.12431e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00039

*** Timestep 40:  t=2.06884e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.8071e-16, 0.0341014, 0.0309929, 3.9442e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341014

   Solving temperature system... 0 iterations.
   Solving porosity system ... 8 iterations.
   Solving peridotite system ... 6 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.01951e-16, 7.36068e-08, 4.15295e-08, 2.44629e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 2.44629e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00040

*** Timestep 41:  t=2.12092e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.29583e-16, 0.0341014, 0.0309591, 3.07655e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341014

   Solving temperature system... 0 iterations.
   Solving porosity system ... 8 iterations.
   Solving peridotite system ... 6 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.8847e-16, 5.49411e-08, 3.26659e-08, 1.92741e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.92741e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00041

*** Timestep 42:  t=2.17299e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.2539e-16, 0.0341015, 0.0309204, 2.40381e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341015

   Solving temperature system... 0 iterations.
   Solving porosity system ... 8 iterations.
   Solving peridotite system ... 6 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.27129e-16, 4.07636e-08, 2.58908e-08, 1.53156e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.53156e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00042

*** Timestep 43:  t=2.22506e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.77842e-16, 0.0341016, 0.030872, 1.88229e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341016

   Solving temperature system... 0 iterations.
   Solving porosity system ... 8 iterations.
   Solving peridotite system ... 6 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.94083e-16, 3.08179e-08, 2.07058e-08, 1.22908e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.22908e-06


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00043

*** Timestep 44:  t=2.27713e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.05625e-16, 0.0341016, 0.0308063, 1.47734e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341016

   Solving temperature system... 0 iterations.
   Solving porosity system ... 7 iterations.
   Solving peridotite system ... 6 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.57276e-16, 2.42189e-08, 1.6705e-08, 9.95788e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 9.95788e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00044

*** Timestep 45:  t=2.3292e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.0024e-16, 0.0341016, 0.0307121, 1.1628e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341016

   Solving temperature system... 0 iterations.
   Solving porosity system ... 7 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.42722e-16, 1.95422e-08, 1.35662e-08, 8.12505e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 8.12505e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00045

*** Timestep 46:  t=2.38127e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.22274e-16, 0.0341016, 0.0305765, 9.18793e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341016

   Solving temperature system... 0 iterations.
   Solving porosity system ... 7 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.93911e-16, 1.56616e-08, 1.10439e-08, 6.65598e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 6.65598e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00046

*** Timestep 47:  t=2.43335e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 3+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.40336e-16, 0.0341017, 0.0303854, 7.3006e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341017

   Solving temperature system... 0 iterations.
   Solving porosity system ... 7 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.53219e-16, 1.22225e-08, 8.96153e-09, 5.40763e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.40763e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00047

*** Timestep 48:  t=2.48542e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.32207e-16, 0.0341017, 0.0301301, 5.82093e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341017

   Solving temperature system... 0 iterations.
   Solving porosity system ... 7 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.01515e-16, 9.35923e-09, 7.20197e-09, 4.35073e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 4.35073e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00048

*** Timestep 49:  t=2.53749e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.27827e-16, 0.0341017, 0.0298196, 4.65751e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341017

   Solving temperature system... 0 iterations.
   Solving porosity system ... 6 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.23705e-16, 7.19512e-09, 5.69561e-09, 3.43491e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 3.43491e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00049

*** Timestep 50:  t=2.58956e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 16 iterations.
   Solving peridotite system ... 16 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.08606e-16, 0.0341017, 0.0294976, 3.72392e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0341017

   Solving temperature system... 0 iterations.
   Solving porosity system ... 6 iterations.
   Solving peridotite system ... 5 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.59731e-16, 5.5915e-09, 4.40742e-09, 2.64507e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 2.64507e-07


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00050

*** Timestep 51:  t=2.6e+16 seconds
   Solving composition reactions in 1 substep(s).
   Solving temperature system... 0 iterations.
   Solving porosity system ... 10 iterations.
   Solving peridotite system ... 10 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 2+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.88475e-16, 0.0068427, 0.00583797, 0.00571631
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0068427

   Solving temperature system... 0 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.76828e-16, 1.10686e-09, 8.60374e-10, 5.38875e-08
      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.38875e-08


   Postprocessing:
     RMS, max velocity:        2.19e-11 m/s, 2.4e-11 m/s
     Writing graphical output: output-mass_reaction_term/solution/solution-00051

Termination requested by criterion: end time


+---------------------------------------------+------------+------------+
+---------------------------------+-----------+------------+------------+
+---------------------------------+-----------+------------+------------+

