Several little things, one is making sponge less verbose.

- Pointing to OBC wiki file from the lateral parameterizations doc.

- Using the MOM6 verbosity to control the time_interp verbosity.

- Making the check for negative water depths more informative.

src/framework/MOM_horizontal_regridding.F90

@@ -740,7 +740,7 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id,  Time, conversion, G, t
 
   if (.not.spongeDataOngrid) then
     if (is_root_pe()) &
-      call time_interp_external(fms_id, Time, data_in, verbose=(verbosity>5), turns=turns)
+      call time_interp_external(fms_id, Time, data_in, verbose=(verbosity>2), turns=turns)
     ! Loop through each data level and interpolate to model grid.
     ! After interpolating, fill in points which will be needed to define the layers.
     do k=1,kd
@@ -840,6 +840,14 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id,  Time, conversion, G, t
 
       call fill_miss_2d(tr_outf, good2, fill2, tr_prev, G, smooth=.true., answers_2018=answers_2018)
 
+      ! now fill in missing values using "ICE-nine" algorithm.
+      tr_outf(:,:) = tr_out(:,:)
+      if (k==1) tr_prev(:,:) = tr_outf(:,:)
+      good2(:,:) = good(:,:)
+      fill2(:,:) = fill(:,:)
+
+      call fill_miss_2d(tr_outf, good2, fill2, tr_prev, G, smooth=.true., answers_2018=answers_2018)
+
 !     if (debug) then
 !       call hchksum(tr_outf, 'field from fill_miss_2d ', G%HI)
 !     endif
@@ -856,7 +864,7 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id,  Time, conversion, G, t
 
     enddo ! kd
   else
-      call time_interp_external(fms_id, Time, data_in, verbose=(verbosity>5), turns=turns)
+      call time_interp_external(fms_id, Time, data_in, verbose=(verbosity>2), turns=turns)
       do k=1,kd
         do j=js,je
           do i=is,ie
@@ -867,6 +875,7 @@ subroutine horiz_interp_and_extrap_tracer_fms_id(fms_id,  Time, conversion, G, t
         enddo
       enddo
   endif
+
 end subroutine horiz_interp_and_extrap_tracer_fms_id
 
 !> Create a 2d-mesh of grid coordinates from 1-d arrays.