We need to save opottemppmdiff as part of ocean_model_z rather than ocean_model. We presently save osaltpmdiff correctly on ocean_model_z. Do you agree? https://github.com/NOAA-GFDL/MOM6-examples/blob/7e5d936adf3879189129607c6f47c005e5a7b04b/ice_ocean_SIS2/OM4_025/diag_table.MOM6#L482