diff --git a/cicecore/cicedynB/dynamics/ice_dyn_evp.F90 b/cicecore/cicedynB/dynamics/ice_dyn_evp.F90
index e411bb35c..f7650e119 100644
--- a/cicecore/cicedynB/dynamics/ice_dyn_evp.F90
+++ b/cicecore/cicedynB/dynamics/ice_dyn_evp.F90
@@ -1254,7 +1254,7 @@ subroutine stress_T   (nx_block,   ny_block,   &
          etax2T   , & ! etax2  = 2*eta  (shear viscous coeff)
          stresspT , & ! sigma11+sigma22
          stressmT , & ! sigma11-sigma22
-         stress12T   ! sigma12
+         stress12T    ! sigma12
 
       real (kind=dbl_kind), dimension (nx_block,ny_block), intent(inout) :: &
          shear    , & ! strain rate II component (1/s)
@@ -1340,7 +1340,7 @@ subroutine stress_T   (nx_block,   ny_block,   &
                               dxT,        dyT,        &
                               tarear    ,             &
                               shear     , divu      , &
-                              rdg_conv  , rdg_shear )
+                              rdg_conv  , rdg_shear   )
 
       endif
 
@@ -1354,7 +1354,7 @@ end subroutine stress_T
 ! Nov 2021      
 
       subroutine stress_U   (nx_block,   ny_block,  & 
-                             ksub,       icellu,    & !!!!new 
+                             ksub,       icellu,    &  
                              indxui,     indxuj,    &
                              uvelE,      vvelE,     &
                              uvelN,      vvelN,     &
@@ -1417,7 +1417,7 @@ subroutine stress_U   (nx_block,   ny_block,  &
 
       character(len=*), parameter :: subname = '(stress_U)'
 
-      do ij = 1, icellu !!!!WATCHOUT!!!!
+      do ij = 1, icellu
          i = indxui(ij)
          j = indxuj(ij)