This python script calculates the stage volume relationship.
It is bit in crude form right now, needed to conver the calculation into function