The nonlinear Schrodinger type equations are of tremendous interest in both theory and applications. Various regimes of pulse propagation in optical fibers are modeled by some form of the nonlinear Schrodinger equation. In this paper we introduce parallel split-step Fourier methods for the numerical simulations of the coupled nonlinear Schrodinger equation that describes the propagation of two orthogonally polarized pulses in a monomode birefringent fibers. These methods are implemented on the Origin 2000 multiprocessor computer. Our numerical experiments have shown that these methods give accurate results and considerable speedup.