# 1: Time step number
# 2: Time (years)
# 3: Time step size (years)
# 4: Number of mesh cells
# 5: Number of Stokes degrees of freedom
# 6: Number of temperature degrees of freedom
# 7: Number of degrees of freedom for all compositions
# 8: Number of nonlinear iterations
# 9: Iterations for temperature solver
# 10: Iterations for Stokes solver
# 11: Velocity iterations in Stokes preconditioner
# 12: Schur complement iterations in Stokes preconditioner
# 13: RMS velocity (m/year)
# 14: Max. velocity (m/year)
# 15: Outward mass flux through boundary with indicator 0 ("left") (kg/yr)
# 16: Outward mass flux through boundary with indicator 1 ("right") (kg/yr)
# 17: Outward mass flux through boundary with indicator 2 ("bottom") (kg/yr)
# 18: Outward mass flux through boundary with indicator 3 ("top") (kg/yr)
# 19: Visualization file name
# 20: Number of advected particles
0 0.000000000000e+00 0.000000000000e+00 100 1003 441 441 1 0 6 40 8 4.08248564e-04 6.91147019e-04 1.65115500e+05 1.65115500e+05 -1.65115500e+05 -1.65115500e+05 output-visco_plastic_yield_strain_weakening_particles/solution/solution-00000 100000 
1 1.000000000000e+00 1.000000000000e+00 100 1003 441 441 1 0 4 33 6 4.08249191e-04 6.91115882e-04 1.65115500e+05 1.65115500e+05 -1.65115500e+05 -1.65115500e+05                                                                            "" 100000 
