***,O2
print,orbitals,civector    !print orbitals and ci-coefficients
geometry={                 !geometry specification, using z-matrix
o1
o2,o1,r
}

r=2.2 bohr                 !bond distance

basis=vtz                  !triple-zeta basis set

{hf                         !invoke RHF program
wf,16,4,2                  !define wavefunction: 16 electrons,
symmetry 4,
                           !triplet
occ,3,1,1,,2,1,1           !number of occupied orbitals in each symmetry
open,1.6,1.7}             !define open shell orbitals

casscf                     !casscf using full valence active space
mrci                       !mrci using full valence casscf reference function
{mrci                       !mrci with only 2p orbitals active, 2s closed
                           !in reference
closed,2,,,,2}             !inactive orbitals in the reference function.