In this work, we study Crank-Nicolson leap-frog (CNLF) methods with fast-slow wave splittings for Navier-Stokes equations (NSE) with a rotation/Coriolis force term, which is a simplification of geophysical flows. We propose a new stabilized CNLF method where the added stabilization completely removes the method's CFL time step condition. A comprehensive stability and error analysis is given. We also prove that for Oseen equations with the rotation term, the unstable mode (for which u(n+1) + u(n-1) equivalent to 0) of CNLF is asymptotically stable. Numerical results are provided to verify the stability and the convergence of the methods.