Application of a multi-step computation scheme to an intratidal estuarine water quality model