# t            : simulation time
# mass_yt      : total simulation mass integrated with yt roughly as : 
#                   [....]
#                   dd = pf.h.all_data()
#                   res = dd.quantities['TotalQuantity']('CylindricalCellMass')[0]
#                   [...]
# mass_flash   : total simulation mass as given by FLASH
# mass_err_rel : relative error

             t    mass_yt  mass_flash  mass_err_rel
0    0.000E+00  1.329E-03   1.416E-03   937.996E-03
1   25.002E-09  1.367E-03   1.412E-03   968.114E-03
2   50.003E-09  1.400E-03   1.400E-03  1000.000E-03
3   75.010E-09  1.412E-03   1.367E-03     1.033E+00
4  100.013E-09  1.439E-03   1.329E-03     1.083E+00