Added mass fill steps after each ssprk3 substep, also altered FCT limiting to use more accurate substep max fluxes
Updated way C is computed, how the masses are computes, and added in sine^4 testing
Changed mfill subroutine to always look to the right when filling mass within elements
Precompute fluxes beforehand to reduce number of overall function calls