
Number of active cells: 192 (on 4 levels)
Number of degrees of freedom: 5,223 (1,666+225+833+833+833+833)

*** Timestep 0:  t=0 years, dt=0 years
   Solving temperature system... 0 iterations.
   Solving composition system ... 0 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 0 iterations.
   Initial Newton Stokes residual = 1.22791e+14, v = 1.22791e+14, p = 0

   Solving Stokes system (GMG)... 15+0 iterations.
      Newton system information: Norm of the rhs: 1.22791e+14
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.71633e-16, 1.61486e-16, 0, 1.56133e-16, 1
      Relative nonlinear residual (total system) after nonlinear iteration 1: 1

   Solving temperature system... 0 iterations.
   Solving composition system ... 0 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 0 iterations.
   Solving Stokes system (GMG)... 11+0 iterations.
      Newton system information: Norm of the rhs: 2.23856e+12, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.71633e-16, 1.61486e-16, 0, 1.56133e-16, 0.0182306
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.0182306

   Solving temperature system... 0 iterations.
   Solving composition system ... 0 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 0 iterations.
   Solving Stokes system (GMG)... 11+0 iterations.
      Newton system information: Norm of the rhs: 6.56739e+10, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.71633e-16, 1.61486e-16, 0, 1.56133e-16, 0.000534843
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.000534843

   Solving temperature system... 0 iterations.
   Solving composition system ... 0 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 0 iterations.
   Solving Stokes system (GMG)... 11+0 iterations.
      Newton system information: Norm of the rhs: 1.62309e+09, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.71633e-16, 1.61486e-16, 0, 1.56133e-16, 1.32183e-05
      Relative nonlinear residual (total system) after nonlinear iteration 4: 1.32183e-05

   Solving temperature system... 0 iterations.
   Solving composition system ... 0 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 0 iterations.
   Solving Stokes system (GMG)... 11+0 iterations.
      Newton system information: Norm of the rhs: 4.47762e+07, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.71633e-16, 1.61486e-16, 0, 1.56133e-16, 3.64653e-07
      Relative nonlinear residual (total system) after nonlinear iteration 5: 3.64653e-07


   Postprocessing:
     Writing graphical output:           output-steinberger_projected_density_dcPicard/solution/solution-00000
     RMS, max velocity:                  0.578 m/year, 1.1 m/year
     Temperature min/avg/max:            273 K, 2437 K, 4250 K
     Heat fluxes through boundary parts: 4.993e+05 W, -4.274e+06 W, 0 W, 0 W
     Writing depth average:              output-steinberger_projected_density_dcPicard/depth_average

*** Timestep 1:  t=100000 years, dt=100000 years
   Solving temperature system... 13 iterations.
   Solving composition system ... 11 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 11 iterations.
   Initial Newton Stokes residual = 8.39912e+13, v = 8.31745e+13, p = 1.16846e+13

   Solving Stokes system (GMG)... 13+0 iterations.
      Newton system information: Norm of the rhs: 1.49315e+13
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00377681, 0.0202636, 0, 0.0108039, 0.177774
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.177774

   Solving temperature system... 11 iterations.
   Solving composition system ... 9 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 9 iterations.
   Solving Stokes system (GMG)... 13+0 iterations.
      Newton system information: Norm of the rhs: 2.55278e+11, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.00012418, 0.000414973, 0, 0.000193229, 0.00303934
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00303934

   Solving temperature system... 10 iterations.
   Solving composition system ... 8 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 8 iterations.
   Solving Stokes system (GMG)... 13+0 iterations.
      Newton system information: Norm of the rhs: 2.2233e+10, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 6.8685e-06, 3.11445e-05, 0, 1.43921e-05, 0.000264707
      Relative nonlinear residual (total system) after nonlinear iteration 3: 0.000264707

   Solving temperature system... 8 iterations.
   Solving composition system ... 7 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 7 iterations.
   Solving Stokes system (GMG)... 12+0 iterations.
      Newton system information: Norm of the rhs: 1.00197e+09, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 6.11361e-07, 2.60566e-06, 0, 1.38988e-06, 1.19295e-05
      Relative nonlinear residual (total system) after nonlinear iteration 4: 1.19295e-05

   Solving temperature system... 7 iterations.
   Solving composition system ... 6 iterations.
   Copying properties into prescribed compositional field density_field... done.
   Solving another_field system ... 5 iterations.
   Solving Stokes system (GMG)... 13+0 iterations.
      Newton system information: Norm of the rhs: 9.33525e+07, Derivative scaling factor: 0
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.31558e-08, 1.25138e-07, 0, 8.13074e-08, 1.11146e-06
      Relative nonlinear residual (total system) after nonlinear iteration 5: 1.11146e-06


   Postprocessing:
     Writing graphical output:           output-steinberger_projected_density_dcPicard/solution/solution-00001
     RMS, max velocity:                  0.578 m/year, 1.12 m/year
     Temperature min/avg/max:            273 K, 2436 K, 4250 K
     Heat fluxes through boundary parts: -2.718e+06 W, -1.136e+06 W, 0 W, 0 W
     Writing depth average:              output-steinberger_projected_density_dcPicard/depth_average

Termination requested by criterion: end time



