psi: (
  label = "6-31G** SCF H2O"
  jobtype = freq
  wfn = scf
  reference = rhf
  multp = 1
  dertype = second
  basis = "6-31G**"
  zmat = (
    o
    h 1 roh
    h 1 roh 2 ahoh
  )
  zvars = (
    roh    0.94305663
    ahoh 105.96876874
  )
)

