Port tardis_opacities script from tardisanalysis (#1580)

* Create opacities.py

Ported over opacities script from old tardis analysis

* Black formatting

* Remove header information