Here I put nontrivial MOLPRO scripts, which I found useful for me and other:

Use the relaxed CCSD one-electron density matrix to calculate the CCSD electrostatic energy (other densities, e.g. MP3 density, can be used as well)

Save correlated densities in the MOLDEN format

Electron affinities through the electron-excitation program on example of EA-EOM-CCSD. "Excitation energies" with "excitation" from a special very diffuse orbital are minus vertical EAs (the second example works with Molpro2010).
If the orbital energy of LUMO is greater than zero (the most common case), use this
If the orbital energy of LUMO is smaller than zero, use this

Electrostatic, induction, 1-order exchange (S²) and exchange-induction components of SAPT in the MCBS (monomer-centered basis set) basis