“…The Dufort-Frankel finite difference approximations do not have severe stability constraints [52,55,58]. In the present solution, the apparent viscosity is a function of local velocity gradient (3U/3Y)j •• This has some ramifications on the stability of the finite-J difference scheme.…”