Contents:
Function calculates the hydrodynamic source terms.
Parent simulation frame
Surface density to be used if not None
Hydrodynamic source terms