# 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 composition solver 1
# 11: Iterations for composition solver 2
# 12: Iterations for Stokes solver
# 13: Velocity iterations in Stokes preconditioner
# 14: Schur complement iterations in Stokes preconditioner
# 15: RMS velocity (m/year)
# 16: Max. velocity (m/year)
# 17: Outward mass flux through boundary with indicator 0 ("left") (kg/yr)
# 18: Outward mass flux through boundary with indicator 1 ("right") (kg/yr)
# 19: Outward mass flux through boundary with indicator 2 ("bottom") (kg/yr)
# 20: Outward mass flux through boundary with indicator 3 ("top") (kg/yr)
# 21: Visualization file name
0 0.000000000000e+00 0.000000000000e+00 100 1003 441 882 1 0 0 0 17 18 36 4.08254250e-04 6.91117759e-04 1.50678150e+05 1.50678150e+05 -1.65115500e+05 -1.35094500e+05 output-visco_plastic_yield_strain_weakening_compositions/solution/solution-00000 
1 1.000000000000e+00 1.000000000000e+00 100 1003 441 882 1 0 5 4 14 15 30 4.08262817e-04 6.91043793e-04 1.50678150e+05 1.50678150e+05 -1.65115500e+05 -1.35094500e+05                                                                               "" 
