
Number of active cells: 512 (on 5 levels)
Number of degrees of freedom: 14,505 (4,290+2,097+4,290+561+2,145+561+561)

*** 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... 38+0 iterations.
   Solving for u_f in 13 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.50805e-16, 1.17203e-16, 0, 0.000956596
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.000956596

   Postprocessing:
     Writing graphical output:           output-average_density_boundary_fluid_pressure/solution/solution-00000.0000
     Compositions min/max/mass:          0/0.07253/8.923e+08 // 0/0/0
     RMS, max velocity:                  8.02e-07 m/year, 4.19e-06 m/year
     Temperature min/avg/max:            293 K, 1210 K, 1624 K
     Mass fluxes through boundary parts: 0 kg/yr, -89.22 kg/yr, 85.73 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 13 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.50805e-16, 1.17203e-16, 0, 7.95357e-09
      Relative nonlinear residual (total system) after nonlinear iteration 2: 7.95357e-09

   Postprocessing:
     Writing graphical output:           output-average_density_boundary_fluid_pressure/solution/solution-00000.0001
     Compositions min/max/mass:          0/0.07253/8.923e+08 // 0/0/0
     RMS, max velocity:                  8.02e-07 m/year, 4.19e-06 m/year
     Temperature min/avg/max:            293 K, 1210 K, 1624 K
     Mass fluxes through boundary parts: 0 kg/yr, -89.22 kg/yr, 85.73 kg/yr, 0 kg/yr


*** Timestep 1:  t=100000 years
   Solving composition reactions in 50 substep(s).
   Solving temperature system... 5 iterations.
   Solving porosity system ... 3 iterations.
   Skipping peridotite composition solve because RHS is zero.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 13+0 iterations.
   Solving for u_f in 13 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00035945, 0.00216443, 0, 7.06593e-06
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00216443

   Postprocessing:
     Writing graphical output:           output-average_density_boundary_fluid_pressure/solution/solution-00001.0000
     Compositions min/max/mass:          -7.948e-05/0.07254/8.922e+08 // 0/0/0
     RMS, max velocity:                  9.25e-07 m/year, 4.58e-06 m/year
     Temperature min/avg/max:            293 K, 1210 K, 1624 K
     Mass fluxes through boundary parts: 0 kg/yr, -366 kg/yr, 369.9 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... 0+0 iterations.
   Solving for u_f in 13 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.69564e-08, 3.61841e-05, 0, 8.82077e-09
      Relative nonlinear residual (total system) after nonlinear iteration 2: 3.61841e-05

   Postprocessing:
     Writing graphical output:           output-average_density_boundary_fluid_pressure/solution/solution-00001.0001
     Compositions min/max/mass:          -7.132e-05/0.07254/8.923e+08 // 0/0/0
     RMS, max velocity:                  9.25e-07 m/year, 4.58e-06 m/year
     Temperature min/avg/max:            293 K, 1210 K, 1624 K
     Mass fluxes through boundary parts: 0 kg/yr, -366 kg/yr, 369.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 13 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.49171e-13, 1.17138e-15, 0, 8.82077e-09
      Relative nonlinear residual (total system) after nonlinear iteration 3: 8.82077e-09

   Postprocessing:
     Writing graphical output:           output-average_density_boundary_fluid_pressure/solution/solution-00001.0002
     Compositions min/max/mass:          -7.132e-05/0.07254/8.923e+08 // 0/0/0
     RMS, max velocity:                  9.25e-07 m/year, 4.58e-06 m/year
     Temperature min/avg/max:            293 K, 1210 K, 1624 K
     Mass fluxes through boundary parts: 0 kg/yr, -366 kg/yr, 369.9 kg/yr, 0 kg/yr


Termination requested by criterion: end time


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

