Added function to find the max value of the flux at each timestep.
Authored by: guidrymw 2012-11-14
Parent: [r78]
Child: [r80]