I should note that the script computes averages for a given C-H vector over all lipids and time points, but does not average the two values for the two vectors on the same C atom; you would need to that. They are not averaged because the equivalence of the two values is an indication of convergence; for a well converged analysis interval, the values should agree fairly closely.