Allow output functions to take a time step and use it to output initial values

All initial values not just the t=0 one, useful for debugging.
Add x^3 and x^2 lines to scatter-err-dts