Fully discrete Jacobi-spherical harmonic spectral method for Navier-Stokes equations