A particle simulation of an interchange instability was performed by taking into account the ion finite Larmor radius (FLR) effects. It is found that the interchange instability with large FLR grows in two phases, that is, linearly growing phase and the nonlinear phase subsequent to the linear phase, where the instability grows exponentially in both phases. The linear growth rates observed in the simulation agree well with the theoretical calculation. The effects of FLR are usually taken in the fluid simulation through the gyroviscosity, the effects of which are verified in the particle simulation with large FLR regime. The gyroviscous cancellation phenomenon observed in the particle simulation causes the drifts in the direction of ion diamagnetic drifts.