run_airglow calculation ignores O(1D) quenching by O

[bharding512]

The run_airglow function calculates 6300 emission. It uses the formulation by Link and Cogger (1988, 1989) which includes quenching by N2 and O2. Recent work (see discussion in Thirupathaiah and Singh (2014)) seems to indicate that quenching by ground state O is actually the dominant loss process at the relevant altitudes. The reaction rate coefficient is 2.5e-11.
O(1D) + O --> O + O

This should be added to the run_airglow function, if it is indeed trustworthy. It will bring VER down significantly, and shift the peak altitude higher.

Link, R., & Cogger, L. L. (1988). A reexamination of the O I 6300-Å nightglow. Journal of Geophysical Research, 93(A9), 9883. https://doi.org/10.1029/JA093iA09p09883

Link, R., & Cogger, L. L. (1989). Correction to ``A reexamination of the O I 6300-Å nightglow’’ by R. Link and L. L. Cogger. Journal of Geophysical Research, 94(A2), 1556. https://doi.org/10.1029/JA094iA02p01556

Thirupathaiah, P., & Singh, V. (2014). An updated model of atomic oxygen redline dayglow emission. Advances in Space Research. https://doi.org/10.1016/j.asr.2014.05.022