Contents:
Function calculates the hydrodynamic source terms for s.max -> assumes that I can callculate the hydrodynamic source term for the product and divide it by Sig1