An implementation of a two-component all-electron ( j j/ ) treatment of both scalar and spinorbit relativistic effects in the MOLFDIR program suite is presented. Relativity is accounted for by Douglas-Kroll transformed one-electron operators: scalar ͑spin-free͒ and so called mean-field spinorbit terms. The interelectronic interaction is represented by the nonrelativistic Coulomb operator. High-level correlated calculations of properties of several systems ͑FO, ClO, Cl, O 2 ϩ , O 2 Ϫ , Tl, and TlH͒ where spin-orbit effects play a dominant role are presented and compared with other data. Agreement with Dirac-Coulomb͑-Gaunt͒ reference values is in general very good.