(*)+Restore USE_WRIGHT_2ND_DERIV_BUG functionality

[Hallberg-NOAA]

This commit restores the effectiveness of the runtime parameter USE_WRIGHT_2ND_DERIV_BUG in determining whether a bug is corrected in the calculation of two of the second derivative terms returned by calculate_density_second_derivs_elem() with the "WRIGHT" equation of state, recreating the behavior (and answers) that are currently on the main branch of MOM6. To do this, it adds and calls the new routine set_params_buggy_Wright() when appropriate, and adds the new element "three" to the buggy_Wright_EOS type. When the bug is fixed, buggy_Wright_EOS%three = 3, but ...%three = 2 to recreate the bug. This commit does change answers for cases using the "WRIGHT" equation of state and one of the "USE_STANLEY_..." parameterizations from those on the dev/gfdl branch of MOM6, but in so doing it restores the answers on the main branch of MOM6. There is also a new publicly visible subroutine.

[alperaltuntas]

This PR fixes the baseline comparison failures for all but one CESM test: All of the tests using our workhorse configuration are now passing, but a quarter-degree test is failing the bitwise reproducibility test. The failing test configuration was initially based on GFDL's quarter-degree setup, with some changes over time. Note, however, that changing DEFAULT_ANSWER_DATE = 99991231 to DEFAULT_ANSWER_DATE = 20231231 fixes the baseline comparison failure for this test as well.

[Hallberg-NOAA]

This has passed pipeline testing at https://gitlab.gfdl.noaa.gov/ogrp/mom6ci/MOM6/-/pipelines/24117. Because that test is comparing MOM_parameter_doc files with the results from dev/gfdl, it is falsely reporting parameter changes, but the answers are bitwise identical.