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

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

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

   Solving temperature system... 0 iterations.
   Skipping porosity composition solve because RHS is zero.
   Skipping peridotite composition solve because RHS is zero.
   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.35084e-16, 0, 0, 0.000116155
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.000116155

   Solving temperature system... 0 iterations.
   Skipping porosity composition solve because RHS is zero.
   Skipping peridotite composition solve because RHS is zero.
   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.78591e-16, 0, 0, 1.1321e-06
      Relative nonlinear residual (total system) after nonlinear iteration 4: 1.1321e-06


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

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

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

   Solving temperature system... 0 iterations.
   Skipping porosity composition solve because RHS is zero.
   Skipping peridotite composition solve because RHS is zero.
   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.35084e-16, 0, 0, 0.000116155
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.000116155

   Solving temperature system... 0 iterations.
   Skipping porosity composition solve because RHS is zero.
   Skipping peridotite composition solve because RHS is zero.
   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.78591e-16, 0, 0, 1.1321e-06
      Relative nonlinear residual (total system) after nonlinear iteration 4: 1.1321e-06


   Postprocessing:
     Compositions min/max/mass:          0/0/0 // 0/0/0
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.843e+07 kg/yr
     Heating rate (average/total):       -1.022e-08 W/kg, -7.919e+04 W
     Temperature min/avg/max:            273 K, 613.4 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00000

*** Timestep 1:  t=2399.9 years
   Solving composition reactions in 11 substep(s).
   Solving temperature system... 18 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.163314, 1, 1, 0.00125006
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1

   Solving temperature system... 13 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000143659, 0.016095, 6.86261e-05, 3.44302e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.016095

   Solving temperature system... 9 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 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.05812e-07, 8.65975e-05, 3.20295e-07, 2.21279e-07
      Relative nonlinear residual (total system) after nonlinear iteration 3: 8.65975e-05

   Solving temperature system... 5 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): 1.03361e-09, 6.04679e-07, 1.92488e-09, 4.04903e-09
      Relative nonlinear residual (total system) after nonlinear iteration 4: 6.04679e-07


   Postprocessing:
     Compositions min/max/mass:          0/0.003666/2.292e+05 // 0/0.003727/2.329e+05
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.841e+07 kg/yr
     Heating rate (average/total):       -1.058e-08 W/kg, -8.187e+04 W
     Temperature min/avg/max:            281.8 K, 634.9 K, 2200 K

*** Timestep 2:  t=4799.9 years
   Solving composition reactions in 12 substep(s).
   Solving temperature system... 18 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.112247, 0.468267, 0.468524, 0.00256933
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.468524

   Solving temperature system... 13 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000101675, 0.000543115, 3.31757e-06, 2.39523e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000543115

   Solving temperature system... 8 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): 1.82287e-07, 4.48436e-05, 2.84584e-07, 5.09574e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 4.48436e-05

   Solving temperature system... 3 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.60909e-10, 3.45271e-17, 3.42082e-17, 5.0942e-08
      Relative nonlinear residual (total system) after nonlinear iteration 4: 5.0942e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.006254/4.475e+05 // 0/0.006357/4.547e+05
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.84e+07 kg/yr
     Heating rate (average/total):       -1.093e-08 W/kg, -8.454e+04 W
     Temperature min/avg/max:            289 K, 656.9 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00001

*** Timestep 3:  t=7199.9 years
   Solving composition reactions in 12 substep(s).
   Solving temperature system... 18 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0853807, 0.31666, 0.316546, 0.000787362
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.31666

   Solving temperature system... 12 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 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): 4.10534e-05, 0.000553082, 5.27187e-07, 2.16979e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000553082

   Solving temperature system... 7 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): 1.00076e-07, 3.02714e-05, 2.51612e-07, 6.04797e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.02714e-05

   Solving temperature system... 3 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.26077e-10, 1.54082e-16, 2.02784e-16, 6.04755e-08
      Relative nonlinear residual (total system) after nonlinear iteration 4: 6.04755e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.008008/6.731e+05 // 0/0.008137/6.84e+05
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.839e+07 kg/yr
     Heating rate (average/total):       -1.128e-08 W/kg, -8.722e+04 W
     Temperature min/avg/max:            296.4 K, 678.9 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00002

*** Timestep 4:  t=9599.9 years
   Solving composition reactions in 12 substep(s).
   Solving temperature system... 18 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 6+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0779469, 0.24331, 0.243454, 0.000639252
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.243454

   Solving temperature system... 13 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 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): 5.18974e-05, 0.000578924, 5.90815e-07, 2.09839e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000578924

   Solving temperature system... 8 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): 1.14741e-07, 2.43373e-05, 2.41752e-07, 3.71228e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.43373e-05

   Solving temperature system... 3 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.03743e-10, 1.09287e-16, 1.91441e-17, 3.71189e-08
      Relative nonlinear residual (total system) after nonlinear iteration 4: 3.71189e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.008949/9.019e+05 // 0/0.009087/9.168e+05
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.838e+07 kg/yr
     Heating rate (average/total):       -1.164e-08 W/kg, -8.989e+04 W
     Temperature min/avg/max:            304 K, 700.8 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00003

*** Timestep 5:  t=10000 years
   Solving composition reactions in 10 substep(s).
   Solving temperature system... 9 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0133082, 0.0604291, 0.0599697, 0.000131895
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0604291

   Solving temperature system... 5 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 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.48842e-06, 0.000547463, 5.19197e-06, 6.00707e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000547463

   Solving temperature system... 3 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): 5.81443e-09, 1.48388e-06, 1.38033e-08, 6.85355e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.48388e-06


   Postprocessing:
     Compositions min/max/mass:          -1.93e-22/0.009538/9.644e+05 // -3.273e-19/0.00968/9.802e+05
     RMS, max velocity:                  1.5 m/year, 1.5 m/year
     Mass fluxes through boundary parts: 0 kg/yr, 0 kg/yr, -3.739e+07 kg/yr, 3.837e+07 kg/yr
     Heating rate (average/total):       -1.17e-08 W/kg, -9.033e+04 W
     Temperature min/avg/max:            305.3 K, 704.4 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00004

Termination requested by criterion: end time


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

