A 1D SOL–divertor plasma simulation code by introducing the anisotropic ion temperature with virtual divertor model has been developed. By introducing the anisotropic ion temperature directly, the second-order derivative parallel ion viscosity term in the momentum transport equation can be excluded and the boundary condition at the divertor plate will not be required in the simulation. In order to express the effects of the divertor plate and accompanying sheath implicitly, a virtual divertor model which has artificial sinks for the particle, momentum and energy has been introduced. Periodic boundary condition becomes available by the use of the virtual divertor model. By using this model, SOL–divertor plasmas which satisfy the Bohm condition has been successfully obtained. The dependence of the ion temperature anisotropy on the normalized mean free path of ion and the validity of the parallel ion viscous flux for the Braginskii expression and the limited one are also investigated.