fill_kd_epbl Subroutine

   subroutine fill_kd_epbl(state_handle, buf)
      !! EPBL interface diffusivity (m^2/s).  Buffer is layer-shaped
      !! (nx, ny, nz); we emit the value at the BOTTOM interface of
      !! each layer (kd_int(:, :, k) convention), losing only the
      !! identically-zero surface interface.
      class(*), intent(in) :: state_handle
      real(wp), intent(inout) :: buf(:, :, :)
      select type (state => state_handle)
      class is (ocean_state_t)
         call fill_kd_epbl_impl(state%epbl%kd_int, buf)
      end select
   end subroutine fill_kd_epbl