A finite difference technique for solving the Oldroyd-B model for 3D-unsteady free surface flows