Added function for mass fluxes
Just realised I meant flow rates and not fluxes