MOM6: +Add ACCURATE_NONBOUS_THICK_CELLO

- NOAA-GFDL/MOM6@e1932f6e7 +Add ACCURATE_NONBOUS_THICK_CELLO
- NOAA-GFDL/MOM6@4706b3e16 +Extract find_dz_for_eta from find_eta
- NOAA-GFDL/MOM6@5d2eba77d +Enhance ice_shelf_diag_mediator like framework
- NOAA-GFDL/MOM6@6fb5dc746 +Fix bugs in ice shelf top slope calculation
- NOAA-GFDL/MOM6@054621cc1 +Add & use optional maskT arg to initialize_masks
- NOAA-GFDL/MOM6@666ff94e0 Add subroutine allocate_segment_field_data
- NOAA-GFDL/MOM6@7df3b18be Move alloc buffer_dst to initialize_segment_data
- NOAA-GFDL/MOM6@5fab3c00b Refactor update_OBC_segment_data loop indices

Author: mom6bot

ice_ocean_SIS2/Baltic_OM5nb_025/MOM_parameter_doc.all

@@ -1849,6 +1849,10 @@ INTWAVE_REMAPPING_USE_OM4_SUBCELLS = True !   [Boolean] default = True
                                 ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT
                                 ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                 ! this option to false.
+ACCURATE_NONBOUS_THICK_CELLO = False !   [Boolean] default = False
+                                ! If true, use the same careful integrals to find the diagnosed non-Boussinesq
+                                ! layer thicknesses as are used to find the free surface height, instead of
+                                ! using an approximate thickness based on division by the mid-layer density.
 
 ! === module MOM_diabatic_driver ===
 ! The following parameters are used for diabatic processes.

ocean_only/nonBous_global/MOM_parameter_doc.all

@@ -1507,6 +1507,10 @@ INTWAVE_REMAPPING_USE_OM4_SUBCELLS = True !   [Boolean] default = True
                                 ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT
                                 ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                 ! this option to false.
+ACCURATE_NONBOUS_THICK_CELLO = False !   [Boolean] default = False
+                                ! If true, use the same careful integrals to find the diagnosed non-Boussinesq
+                                ! layer thicknesses as are used to find the free surface height, instead of
+                                ! using an approximate thickness based on division by the mid-layer density.
 
 ! === module MOM_diabatic_driver ===
 ! The following parameters are used for diabatic processes.