There are shell/Pythons scripts commonly used in Gaussian software calculations

LIIC is linearly interpolated internal coordinate (liic) script writte by python

MECI is This is script which can optimize Conical Intersections without Derivative Coupling Vectors with GAUSSIAN16 interface. ref: "J. Phys. Chem. B 2008, 112, 405-413"