def METHOD_NAME(inp_h2s, datadir):
    mol = psi4.geometry("""
            0 1
            S
            H   1   1.350490
            H   1   1.350490   2   96.061977
            """)
    psi4.set_options({
        'scf_type': 'pk',
        'e_convergence': 10
        })
    psi4.set_options(inp_h2s['options'])
    molden_file = f"{inp_h2s['name']}.molden"
    ref = datadir.join(f"{inp_h2s['name']}.ref")
    e, wfn = psi4.energy('scf', return_wfn=True, molecule=mol)
    wfn.write_molden(molden_file, do_virtual=True, use_natural=False)
    assert psi4.compare_moldenfiles(ref, molden_file)