diff --git a/src/ALE/coord_scalar.F90 b/src/ALE/coord_scalar.F90
index 105aebd005..e235a536ec 100644
--- a/src/ALE/coord_scalar.F90
+++ b/src/ALE/coord_scalar.F90
@@ -158,9 +158,6 @@ subroutine build_scalar_column(CS, nz, depth, h, T, S, eqn_of_state, z_interface
       densities(k) = T(mapping(k))
     enddo
 
-    print *, "xTmp before sort", xTmp
-    print *, "h_nv before sort", h_nv
-
     ! Sort densities and get scalar array for sorting
     if ( CS%needs_sorting ) then
       call sort_scalar_k_1d(count_nonzero_layers, densities, ksort)
@@ -177,17 +174,10 @@ subroutine build_scalar_column(CS, nz, depth, h, T, S, eqn_of_state, z_interface
       xTmp(k+1) = xTmp(k) + h_nv(k)
     enddo
 
-    print *, "xTmp after sort", xTmp
-    print *, "h_nv after sort", h_nv
-
     ! Based on source column density profile, interpolate to generate a new grid
     call build_and_interpolate_grid(CS%interp_CS, densities, count_nonzero_layers, &
                                     h_nv, xTmp, CS%target_density, CS%nk, h_new, &
                                     x1, h_neglect, h_neglect_edge)
-    
-    print *, "densities", densities
-    print *, "ksort", ksort
-    print *, "h_new", h_new
 
     ! Inflate vanished layers
     call old_inflate_layers_1d(CS%min_thickness, CS%nk, h_new)
diff --git a/src/framework/MOM_diag_remap.F90 b/src/framework/MOM_diag_remap.F90
index 66116a331c..75ab36e4e5 100644
--- a/src/framework/MOM_diag_remap.F90
+++ b/src/framework/MOM_diag_remap.F90
@@ -458,7 +458,7 @@ subroutine diag_remap_do_remap(remap_cs, G, GV, h, staggered_in_x, staggered_in_
         if ( needs_sorting ) then
           ! Get sorting indices
           ksort(:) = remap_cs%ksort3d(i,j,:)
-          ! Sort field to be remapped
+          ! Sort source thicknesses, and field to be remapped
           do k=1, nz_src
             h_src(k) = h(i,j,ksort(k))
             field1d(k) = field(i,j,ksort(k))