Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In the present work, the accuracy of a Fractional step method, which is widely used in LES simulation, has been increased to the fourth-order accurate compact Pade discretization. To validate the present code, the flow-field past a cylinder was simulated and compared with experiment. A good agreement with experiment was achieved.