Fixes bugs in CMix/KPP interface

[vanroekel]

Two parts of the KPP interface had bugs

1. The computation of surface buoyancy forcing and friction velocity were incorrect. The surfaceBuoyancyForcing subtracted the wrong values from the heat fluxes from thickness fluxes. The surfaceFrictionVelocity had incorrect interpolation to centers, the interpolation itself was unneeded. MPAS was interpolating stress to edges (for momentum forcing, but interpolating back (instead of using coupler values) for mixing friction velocity.
2. The KPP interface for matchBoth is fixed and the indexing is fixed for matching interior to boundary layer diffusivities.

[vanroekel]

zppy results from a 50 year LR spinup B-case are here

https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.vanroekel/E3SMv3/20240918_kppTest_ne30pg2_icos30_anvil/

it includes mpas-analysis and e3sm-diags

the comparison run is

https://web.lcrc.anl.gov/public/e3sm/diagnostic_output/ac.wlin/E3SMv3/20231209.v3.LR.piControl-spinup.chrysalis/

Looking at the diags, the changes to the climate are minor. Some degradation in bias (e.g. N. Atlantic temperature/salinity near surface), some improvement (e.g., drake passage transport).

[vanroekel]

Further discussion and testing is here - E3SM-Ocean-Discussion#107

[bpbond]

Tagging @njeffery for her disciplinary expertise here 🙏

[njeffery]

Approved based on review of mpas-analysis changes to sea ice and ocean metrics. After 50 years, changes are minor.