diff --git a/src/core_ocean/shared/mpas_ocn_diagnostics.F b/src/core_ocean/shared/mpas_ocn_diagnostics.F
index 1a2f03590d..de1ad3674e 100644
--- a/src/core_ocean/shared/mpas_ocn_diagnostics.F
+++ b/src/core_ocean/shared/mpas_ocn_diagnostics.F
@@ -1859,7 +1859,6 @@ subroutine ocn_reconstruct_gm_vectors(diagnosticsPool, meshPool) !{{{
          transportVelocityMeridional, normalGMBolusVelocity, GMBolusVelocityX, GMBolusVelocityY, GMBolusVelocityZ, &
          GMBolusVelocityZonal, GMBolusVelocityMeridional, &
          gmStreamFuncTopOfEdge, GMStreamFuncX, GMStreamFuncY, GMStreamFuncZ, GMStreamFuncZonal, GMStreamFuncMeridional
-      real(kind=RKIND), dimension(:,:,:), pointer :: slopeTriadUp, slopeTriadDown
 
       call mpas_timer_start('reconstruct gm vecs')
 
@@ -1884,12 +1883,7 @@ subroutine ocn_reconstruct_gm_vectors(diagnosticsPool, meshPool) !{{{
       call mpas_pool_get_array(diagnosticsPool, 'GMStreamFuncZonal', GMStreamFuncZonal)
       call mpas_pool_get_array(diagnosticsPool, 'GMStreamFuncMeridional', GMStreamFuncMeridional)
 
-      call mpas_pool_get_array(diagnosticsPool, 'slopeTriadUp', slopeTriadUp)
-      call mpas_pool_get_array(diagnosticsPool, 'slopeTriadDown', slopeTriadDown)
-
-      ! mrp temporary change to see slope
-      !call mpas_reconstruct(meshPool, normalTransportVelocity,          &
-      call mpas_reconstruct(meshPool, slopeTriadUp(:,1,:),          &
+      call mpas_reconstruct(meshPool, normalTransportVelocity,          &
                        transportVelocityX,            &
                        transportVelocityY,            &
                        transportVelocityZ,            &
@@ -1897,9 +1891,7 @@ subroutine ocn_reconstruct_gm_vectors(diagnosticsPool, meshPool) !{{{
                        transportVelocityMeridional    &
                       )
 
-      ! mrp temporary change to see slope
-      !call mpas_reconstruct(meshPool, normalGMBolusVelocity,          &
-      call mpas_reconstruct(meshPool, slopeTriadUp(:,2,:),          &
+      call mpas_reconstruct(meshPool, normalGMBolusVelocity,          &
                        GMBolusVelocityX,            &
                        GMBolusVelocityY,            &
                        GMBolusVelocityZ,            &
@@ -1907,9 +1899,7 @@ subroutine ocn_reconstruct_gm_vectors(diagnosticsPool, meshPool) !{{{
                        GMBolusVelocityMeridional    &
                       )
 
-      ! mrp temporary change to see slope
-      !call mpas_reconstruct(meshPool, gmStreamFuncTopOfEdge,          &
-      call mpas_reconstruct(meshPool, slopeTriadUp(:,1,:),          &
+      call mpas_reconstruct(meshPool, gmStreamFuncTopOfEdge,          &
                       GMStreamFuncX,            &
                       GMStreamFuncY,            &
                       GMStreamFuncZ,            &
diff --git a/src/core_ocean/shared/mpas_ocn_equation_of_state_linear.F b/src/core_ocean/shared/mpas_ocn_equation_of_state_linear.F
index ab2e7a8660..d8b2a4f555 100644
--- a/src/core_ocean/shared/mpas_ocn_equation_of_state_linear.F
+++ b/src/core_ocean/shared/mpas_ocn_equation_of_state_linear.F
@@ -299,9 +299,6 @@ subroutine ocn_equation_of_state_linear_density_exp(           &
 
       ! copy some intent(in) into local work space
 
-            ! mrp printing debug loop. Remove later
-            !print *, 'EOS. displacementType=',trim(displacementType)
-            !    print *, 'k,density(k,iCell),thermalExpansionCoeff(k,iCell),salineContractionCoeff(k,iCell)'
       !*** if surfaceDisplaced, then compute density at all levels
       !*** based on input surface values
       !*** return fixed values of expansion coefficients
@@ -320,18 +317,11 @@ subroutine ocn_equation_of_state_linear_density_exp(           &
                                (tracersSurfaceLayerValue(indexS,iCell) - &
                                 ocnEqStateLinearSref)
 
-            !mrp replace this ********** temporary
-!            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha/ &
-!                                             density(k,iCell)
-!
-!            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta/ &
-!                                              density(k,iCell)
-            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha
-            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta
-            ! mrp printing debug loop. Remove later
-            !if (iCell==6) then
-            !    print *, k,density(k,iCell),thermalExpansionCoeff(k,iCell),salineContractionCoeff(k,iCell)
-            !endif
+            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha/ &
+                                             density(k,iCell)
+
+            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta/ &
+                                              density(k,iCell)
          end do
          end do
          !$omp end do
@@ -350,24 +340,16 @@ subroutine ocn_equation_of_state_linear_density_exp(           &
                                (tracers(indexS,k,iCell) - &
                                 ocnEqStateLinearSref)
 
-            !mrp replace this ********** temporary
-!            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha/ &
-!                                             density(k,iCell)
-!
-!            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta/ &
-!                                              density(k,iCell)
-            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha
-            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta
-            ! mrp printing debug loop. Remove later
-            !if (iCell==6) then
-            !    print *, k,density(k,iCell),thermalExpansionCoeff(k,iCell),salineContractionCoeff(k,iCell)
-            !endif
+            thermalExpansionCoeff(k,iCell) = ocnEqStateLinearAlpha/ &
+                                             density(k,iCell)
+
+            salineContractionCoeff(k,iCell) = ocnEqStateLinearBeta/ &
+                                              density(k,iCell)
          end do
          end do
          !$omp end do
 
       endif
-      
 
       !-----------------------------------------------------------------
 
diff --git a/src/core_ocean/shared/mpas_ocn_gm.F b/src/core_ocean/shared/mpas_ocn_gm.F
index 6dedaabe33..8060e455db 100644
--- a/src/core_ocean/shared/mpas_ocn_gm.F
+++ b/src/core_ocean/shared/mpas_ocn_gm.F
@@ -508,25 +508,10 @@ subroutine ocn_GM_compute_Bolus_velocity(statePool, diagnosticsPool, &
          ! Normalize k33
          do k = 2, maxLevelCell(iCell)
             k33(k,iCell) = k33(k,iCell) / k33Norm(k)
-                !print '(a,2i5,1000e9.2)', 'iCell,k,Up',iCell,k,slopeTriadUp(k,:,iCell)
-                !print '(a,2i5,1000e9.2)', 'iCell,k,Dn',iCell,k,slopeTriadDown(k,:,iCell)
-if (k33Norm(k)<1e-15) then
-    print *, 'WARNING: k33Norm(k)<1e-15'
-                !print '(a,100i9)', 'k,iCell: ',k,iCell
-                !print '(a,100e9.2)', 'k33: ',k33(k,iCell)
-                !print '(a,100e9.2)', 'k33Norm: ',k33Norm(k)
-                !print '(a,100f5.2)', 'slopeTaper: ',slopeTaper(k)
-                !print '(a,1000e9.2)', 'slopeTriadUp: ',slopeTriadUp(k,:,iCell)
-                !print '(a,1000e9.2)', 'slopeTriadDown: ',slopeTriadDown(k,:,iCell)
-endif
          end do
       end do ! iCell
       !$omp end do
 
-!      print *, 'in gm: max/min, up/down',&
-!      maxval(slopeTriadUp),minval(slopeTriadUp), &
-!      maxval(slopeTriadDown),minval(slopeTriadDown)
-
       deallocate(dzTop)
       deallocate(dTdzTop)
       deallocate(dSdzTop)