Upwind Formulation Finite-Difference

pull from PCHinR