diff --git a/cicecore/cicedyn/dynamics/ice_dyn_evp.F90 b/cicecore/cicedyn/dynamics/ice_dyn_evp.F90
index cf111cccf..a24c8f57d 100644
--- a/cicecore/cicedyn/dynamics/ice_dyn_evp.F90
+++ b/cicecore/cicedyn/dynamics/ice_dyn_evp.F90
@@ -763,12 +763,14 @@ subroutine evp (dt)
                                                  icellE    (iblk),                 &
                                                  indxEi  (:,iblk), indxEj(:,iblk), &
                                                  vice  (:,:,iblk), aice(:,:,iblk), &
-                                                 hwater(:,:,iblk), TbE (:,:,iblk))
+                                                 hwater(:,:,iblk), TbE (:,:,iblk), &
+                                                 grid_location='E')
                   call seabed_stress_factor_LKD (nx_block        , ny_block,       &
                                                  icellN    (iblk),                 &
                                                  indxNi  (:,iblk), indxNj(:,iblk), &
                                                  vice  (:,:,iblk), aice(:,:,iblk), &
-                                                 hwater(:,:,iblk), TbN (:,:,iblk))
+                                                 hwater(:,:,iblk), TbN (:,:,iblk), &
+                                                 grid_location='N')
                enddo
                !$OMP END PARALLEL DO
 
diff --git a/cicecore/cicedyn/dynamics/ice_dyn_shared.F90 b/cicecore/cicedyn/dynamics/ice_dyn_shared.F90
index 50f1aae6e..69e552730 100644
--- a/cicecore/cicedyn/dynamics/ice_dyn_shared.F90
+++ b/cicecore/cicedyn/dynamics/ice_dyn_shared.F90
@@ -1289,7 +1289,7 @@ subroutine seabed_stress_factor_LKD (nx_block, ny_block,         &
       real (kind=dbl_kind), dimension (nx_block,ny_block), intent(inout) :: &
          TbU           ! seabed stress factor at 'grid_location' (N/m^2)
 
-      character(len=*), optional, intent(inout) :: &
+      character(len=*), optional, intent(in) :: &
          grid_location ! grid location (U, E, N), U assumed if not present
 
       real (kind=dbl_kind) :: &