
Loading shared library <./libparticle_melt_advection.so>

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.
   Solving peridotite system ... 0 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 22+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.77121e-16, 1.4111e-16, 1.4111e-16, 0.168325
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.168325

   Solving temperature system... 0 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 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.77121e-16, 1.4111e-16, 1.4111e-16, 6.14908e-08
      Relative nonlinear residual (total system) after nonlinear iteration 2: 6.14908e-08


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00000
     Compositions min/max/mass: 1.438e-28/0.2/1.257e+08 // 1.438e-28/0.2/1.257e+08
     RMS, max velocity:         0.00292 m/year, 0.00958 m/year
     Temperature min/avg/max:   293 K, 1519 K, 1625 K
     Writing depth average:     output-particle_melt_advection/depth_average.gnuplot
     Writing particle output:   output-particle_melt_advection/particles/particle-00000

*** Timestep 1:  t=3221.68 years
   Solving temperature system... 7 iterations.
   Solving porosity system ... 4 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 18+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 0.000236169, 0.0291878, 0.00606522, 0.0180344
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0291878

   Solving temperature system... 5 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 14+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.80631e-07, 0.00179259, 0.000110766, 0.000714395
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.00179259

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 10+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.26763e-08, 8.26168e-05, 3.9685e-06, 2.30822e-05
      Relative nonlinear residual (total system) after nonlinear iteration 3: 8.26168e-05

   Solving temperature system... 3 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.35566e-09, 5.2739e-06, 2.4841e-07, 2.03704e-06
      Relative nonlinear residual (total system) after nonlinear iteration 4: 5.2739e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00001
     Compositions min/max/mass: -6.192e-07/0.1974/1.257e+08 // -2.807e-09/0.1976/1.256e+08
     RMS, max velocity:         0.00293 m/year, 0.00964 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00001

*** Timestep 2:  t=6396.64 years
   Solving temperature system... 7 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 15+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.47535e-05, 0.00281221, 0.0030857, 0.00444249
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00444249

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 11+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.81014e-08, 0.000127046, 5.59085e-06, 6.94992e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 0.000127046

   Solving temperature system... 3 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 6+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.29853e-09, 6.30943e-06, 3.6724e-07, 2.65502e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 6.30943e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00002
     Compositions min/max/mass: -2.406e-06/0.1939/1.257e+08 // -1.94e-08/0.1944/1.256e+08
     RMS, max velocity:         0.00293 m/year, 0.00966 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00002

*** Timestep 3:  t=9540.55 years
   Solving temperature system... 6 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 15+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 9.18843e-06, 0.00141345, 0.00155608, 0.00416081
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00416081

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 11+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.40067e-08, 9.63422e-05, 3.03379e-06, 5.76048e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 9.63422e-05

   Solving temperature system... 3 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 1.39156e-09, 3.80954e-06, 2.17462e-07, 1.73423e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 3.80954e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00003
     Compositions min/max/mass: -4.284e-06/0.1902/1.257e+08 // -2.254e-08/0.1906/1.256e+08
     RMS, max velocity:         0.00294 m/year, 0.00965 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00003

*** Timestep 4:  t=12647.8 years
   Solving temperature system... 6 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 15+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.44241e-06, 0.00149899, 0.000419693, 0.00403189
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00403189

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 11+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.94299e-08, 9.05579e-05, 2.39467e-06, 5.54046e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 9.05579e-05

   Solving temperature system... 3 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 9.07644e-10, 2.98627e-06, 1.57809e-07, 1.43802e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.98627e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00004
     Compositions min/max/mass: -7.613e-06/0.1868/1.257e+08 // -3.478e-08/0.1867/1.255e+08
     RMS, max velocity:         0.00294 m/year, 0.00964 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00004

*** Timestep 5:  t=15718.5 years
   Solving temperature system... 6 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 14+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.90867e-06, 0.00159914, 0.000195392, 0.0039633
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0039633

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 10+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.98838e-08, 8.93856e-05, 2.31255e-06, 5.65609e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 8.93856e-05

   Solving temperature system... 3 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 5+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 5.32824e-10, 2.67155e-06, 1.29802e-07, 1.40222e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.67155e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00005
     Compositions min/max/mass: -1.52e-05/0.1838/1.257e+08 // -4.951e-08/0.1831/1.255e+08
     RMS, max velocity:         0.00294 m/year, 0.0097 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00005

*** Timestep 6:  t=18754.4 years
   Solving temperature system... 6 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 16+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 2.68937e-06, 0.00163532, 0.00034793, 0.00393132
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.00393132

   Solving temperature system... 4 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 10+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.91019e-08, 8.6807e-05, 2.29183e-06, 5.71671e-05
      Relative nonlinear residual (total system) after nonlinear iteration 2: 8.6807e-05

   Solving temperature system... 2 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 4+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 3.95955e-10, 2.51001e-06, 1.15786e-07, 1.40684e-06
      Relative nonlinear residual (total system) after nonlinear iteration 3: 2.51001e-06


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00006
     Compositions min/max/mass: -2.693e-05/0.1812/1.257e+08 // -6.669e-08/0.18/1.254e+08
     RMS, max velocity:         0.00295 m/year, 0.00973 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00006

*** Timestep 7:  t=20000 years
   Solving temperature system... 6 iterations.
   Solving porosity system ... 3 iterations.
   Solving peridotite system ... 3 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 14+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 7.29072e-07, 0.000482218, 0.000180815, 0.0011325
      Relative nonlinear residual (total system) after nonlinear iteration 1: 0.0011325

   Solving temperature system... 3 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 2 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 7+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.80229e-09, 1.23207e-05, 3.46821e-07, 8.50963e-06
      Relative nonlinear residual (total system) after nonlinear iteration 2: 1.23207e-05

   Solving temperature system... 2 iterations.
   Solving porosity system ... 2 iterations.
   Solving peridotite system ... 1 iterations.
   Rebuilding Stokes preconditioner...
   Solving Stokes system... 1+0 iterations.
   Solving for u_f in 12 iterations.
      Relative nonlinear residuals (temperature, compositional fields, Stokes system): 4.47937e-11, 1.72201e-07, 7.68501e-09, 1.14394e-07
      Relative nonlinear residual (total system) after nonlinear iteration 3: 1.72201e-07


   Postprocessing:
     Writing graphical output:  output-particle_melt_advection/solution/solution-00007
     Compositions min/max/mass: -4.328e-05/0.1802/1.257e+08 // -7.47e-08/0.1789/1.254e+08
     RMS, max velocity:         0.00295 m/year, 0.00974 m/year
     Temperature min/avg/max:   293 K, 1518 K, 1625 K
     Writing particle output:   output-particle_melt_advection/particles/particle-00007

Termination requested by criterion: end time


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

