add a function for mass flow interpolation at a station, but not finished