Initialize solution_rate,instantiate interior_rate
Use solution rates now in mass_residual

This is an API-breaking change for many DiffSystem-based codes