Volume conservation in numerical methods for open channel of arbitrary cross section