Add nonlinear loss to diags

[mnlevy1981]

@cpetrik has asked us to provide the nonlinear zoo_loss term integrated over the top 150m of the ocean to use as a forcing field in FEISTY, rather than trying to estimate it based on the integral of total zoo_loss. In this PR, I break up the zoo_loss computation into zoo_loss_linear and zoo_loss_nonlinear (which will introduce round-off level changes because each of those terms is multiplied by Tfunc and then zoo_loss is the sum of the two, whereas before zoo_loss was the sum of the two without Tfunc and then the sum itself was multipled), and also introduce four new per-zoo diagnostics:

1. The full 3D linear contribution to zoo_loss
2. The vertical integral of the linear contribution to zoo_loss over the top 150 meters
3. The full 3D nonlinear contribution to zoo_loss
4. The vertical integral of the nonlinear contribution to zoo_loss over the top 150 meters

By default, the 150m integral of the nonlinear term is in the high frequency output, but neither the integral of the linear term nor either of the 3D fields are written out.

This branch is from a CESM 2.2.2 run, so it's pretty out of date. I'll merge in the latest development and resolve conflicts, and I'll probably also add the 100m vertical integrals for completeness since we have zoo_loss_zint_100m defined; it's omitted in the initial commits because we don't need them for the 0.1° run we are doing.