Corrections to thermosphere heating diagnostics

[fvitt]

Issue Type

Other (please describe below)

Issue Description

From Hanli Liu:

Suggested changes to QEI/QEN/QIN history fields In src/physics/waccmx/ion_electron_temp.F90:

1. qei/qen/qin are all with unit ev/g/s. We may want to add a comment about their unit at the beginning of the subroutine when they are declared.

2. Add cpairv to use physconst, only : rairv, mbarv ! Constituent dependent rair and mbar

3. Before adding qei/qen/qin to history, convert them to K/s to be consistent with other heating/cooling terms. qei_out = qei/sToQConv/cpairv qen_out = qen/sToQConv/cpairv qin_out = qin/sToQConv/cpairv

4. Change/add unit for QEN/QEI/QIN to K/s.

Will this change answers?

No

Will you be implementing this yourself?

Yes