  Here, we present simple and efficient numerical scheme to study static and
dynamic properties of spin-1 Bose-Einstein condensates (BECs) with spin-orbit
(SO) coupling by solving three coupled Gross-Pitaevskii equations (CGPEs) in
three-, quasi-two and quasi-one dimensional systems. We provide a set of three
codes developed in FORTRAN 90/95 programming language with user defined '{\em
option}' of imaginary and real-time propagation. We present the numerical
results for energy, chemical potentials, and component densities for the ground
state and compare with the available results from the literature. The results
are presented for both the ferromagnetic and antiferromagnetic spin-1 BECs with
and without SO coupling. To improve the computational speed, all the codes have
the option of OpenMP parallelization. We have also presented the results for
speedup and efficiency of OpenMP parallelization for the three codes with both
imaginary and real-time propagation.
