
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, dt=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 (AMG)... 4+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.75134e-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 (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60814e-16, 0, 0, 0.000359498
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000359498

   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 (AMG)... 3+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.58765e-16, 0, 0, 4.81749e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 4.81749e-06


Skipping mesh refinement, because the mesh did not change.

*** Timestep 0:  t=0 years, dt=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 (AMG)... 4+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.75134e-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 (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.60814e-16, 0, 0, 0.000359498
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000359498

   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 (AMG)... 3+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.58765e-16, 0, 0, 4.81749e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 4.81749e-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, dt=2399.9 years
   Solving composition reactions... in 1 substep(s).
   Solving temperature system... 19 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 4+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.193675, 1, 1, 0.000180061
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1

   Solving temperature system... 14 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000273373, 0.0160998, 5.82834e-05, 2.27544e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0160998

   Solving temperature system... 10 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.84492e-07, 9.41328e-05, 2.8301e-07, 1.63849e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 9.41328e-05

   Solving temperature system... 5 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.12209e-09, 2.00805e-16, 1.09858e-17, 1.63353e-08
      Relative nonlinear residual (total system) after nonlinear iteration 4: 1.63353e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.003104/2.004e+05 // 0/0.003155/2.037e+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.059e-08 W/kg, -8.201e+04 W
     Temperature min/avg/max:            282.7 K, 636 K, 2200 K

*** Timestep 2:  t=4799.9 years, dt=2400 years
   Solving composition reactions... in 1 substep(s).
   Solving temperature system... 16 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0734832, 0.0860908, 0.0872861, 0.00127401
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0872861

   Solving temperature system... 12 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 4+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 8.62655e-05, 0.00920947, 4.00947e-05, 1.1267e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00920947

   Solving temperature system... 7 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.57025e-08, 1.58444e-05, 9.1875e-08, 2.08278e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.58444e-05

   Solving temperature system... 2 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.06744e-11, 1.51499e-16, 1.15832e-18, 2.08281e-08
      Relative nonlinear residual (total system) after nonlinear iteration 4: 2.08281e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.005795/4.14e+05 // 0/0.00589/4.206e+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.095e-08 W/kg, -8.472e+04 W
     Temperature min/avg/max:            289.4 K, 658.3 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00001

*** Timestep 3:  t=7199.9 years, dt=2400 years
   Solving composition reactions... in 1 substep(s).
   Solving temperature system... 16 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0565606, 0.0922834, 0.092121, 0.00234648
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0922834

   Solving temperature system... 10 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 3+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.92024e-05, 0.000765469, 7.04536e-07, 1.82526e-07
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000765469

   Solving temperature system... 6 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.77364e-08, 3.59105e-06, 2.91189e-08, 6.15163e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.59105e-06


   Postprocessing:
     Compositions min/max/mass:          0/0.007724/6.382e+05 // 0/0.007849/6.485e+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.131e-08 W/kg, -8.741e+04 W
     Temperature min/avg/max:            296.4 K, 680.4 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00002

*** Timestep 4:  t=9599.9 years, dt=2400 years
   Solving composition reactions... in 1 substep(s).
   Solving temperature system... 15 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 6+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0403663, 0.0775667, 0.0778858, 0.00045208
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0778858

   Solving temperature system... 10 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.24851e-05, 5.46172e-05, 2.03284e-08, 3.94331e-08
      Relative nonlinear residual (total system) after nonlinear iteration 2: 5.46172e-05

   Solving temperature system... 6 iterations.
   Solving porosity system ... 0 iterations.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.76633e-08, 1.13491e-16, 5.71636e-17, 3.94213e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.94213e-08


   Postprocessing:
     Compositions min/max/mass:          0/0.008836/8.676e+05 // 0/0.008972/8.819e+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.166e-08 W/kg, -9.008e+04 W
     Temperature min/avg/max:            304 K, 702.4 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00003

*** Timestep 5:  t=10000 years, dt=400.102 years
   Solving composition reactions... in 1 substep(s).
   Solving temperature system... 7 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 5+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.0039363, 0.182191, 0.179593, 0.000372748
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.182191

   Solving temperature system... 6 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 4+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 8.27843e-06, 0.00253913, 2.94151e-05, 5.95329e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00253913

   Solving temperature system... 3 iterations.
   Solving porosity system ... 1 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system (AMG)... 0+0 iterations.
   Solving fluid velocity system... 1 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.86209e-09, 2.30132e-06, 1.87265e-08, 2.24129e-08
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.30132e-06


   Postprocessing:
     Compositions min/max/mass:          -1.563e-13/0.01207/1.046e+06 // -2.494e-13/0.01225/1.063e+06
     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.835e+07 kg/yr
     Heating rate (average/total):       -1.172e-08 W/kg, -9.048e+04 W
     Temperature min/avg/max:            305.3 K, 706.1 K, 2200 K
     Writing graphical output:           output-shear_heating_with_melt/solution/solution-00004

Termination requested by criterion: end time



