
Number of active cells: 256 (on 5 levels)
Number of degrees of freedom: 7,369 (2,178+1,057+2,178+289+1,089+289+289)

*** Timestep 0:  t=0 years
   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 35+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.49526e-16, 2.30078e-16, 0, 0.000973208
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.000973208

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00000.0000
     Compositions min/max/mass:          0/0.07253/9.981e+08 // 0/0/0
     RMS, max velocity:                  1.65e-06 m/year, 1.19e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -263.6 kg/yr, 273.7 kg/yr, 0 kg/yr

   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.49526e-16, 2.30078e-16, 0, 5.67468e-09
      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.67468e-09

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00000.0001
     Compositions min/max/mass:          0/0.07253/9.981e+08 // 0/0/0
     RMS, max velocity:                  1.65e-06 m/year, 1.19e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -263.6 kg/yr, 273.7 kg/yr, 0 kg/yr


*** Timestep 1:  t=100000 years
   Solving composition reactions in 50 substep(s).
   Solving temperature system... 6 iterations.
   Solving porosity system ... 4 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 13+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000316165, 0.00699099, 0, 1.18977e-05
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00699099

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00001.0000
     Compositions min/max/mass:          -0.0002415/0.07229/9.98e+08 // 0/0/0
     RMS, max velocity:                  1.82e-06 m/year, 1.27e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -977.4 kg/yr, 979 kg/yr, 0 kg/yr

   Solving temperature system... 3 iterations.
   Solving porosity system ... 3 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.85176e-08, 0.000102628, 0, 2.33702e-08
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000102628

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00001.0001
     Compositions min/max/mass:          -0.0002302/0.07229/9.981e+08 // 0/0/0
     RMS, max velocity:                  1.82e-06 m/year, 1.27e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -972.7 kg/yr, 983.9 kg/yr, 0 kg/yr

   Solving temperature system... 2 iterations.
   Solving porosity system ... 2 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60675e-10, 1.20275e-07, 0, 9.8038e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.20275e-07

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00001.0002
     Compositions min/max/mass:          -0.0002302/0.07229/9.981e+08 // 0/0/0
     RMS, max velocity:                  1.82e-06 m/year, 1.27e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -972.7 kg/yr, 983.9 kg/yr, 0 kg/yr

   Solving temperature system... 0 iterations.
   Solving porosity system ... 0 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 0+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.91918e-14, 1.96332e-14, 0, 9.8038e-09
      Relative nonlinear residual (total system) after nonlinear iteration 4: 9.8038e-09

   Postprocessing:
     Writing graphical output:           output-melt_and_traction/solution/solution-00001.0003
     Compositions min/max/mass:          -0.0002302/0.07229/9.981e+08 // 0/0/0
     RMS, max velocity:                  1.82e-06 m/year, 1.27e-05 m/year
     Temperature min/avg/max:            293 K, 1433 K, 1677 K
     Mass fluxes through boundary parts: 0 kg/yr, -972.7 kg/yr, 983.9 kg/yr, 0 kg/yr


Termination requested by criterion: end time


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

