diff --git a/cicecore/cicedynB/dynamics/ice_dyn_eap.F90 b/cicecore/cicedynB/dynamics/ice_dyn_eap.F90 index 9c52bb888..9ecc79305 100644 --- a/cicecore/cicedynB/dynamics/ice_dyn_eap.F90 +++ b/cicecore/cicedynB/dynamics/ice_dyn_eap.F90 @@ -326,11 +326,18 @@ subroutine eap (dt) do j = 1, ny_block do i = 1, nx_block if (icetmask(i,j,iblk)==0) then + if (tmask(i,j,iblk)) then ! structure tensor - a11_1(i,j,iblk) = p5 - a11_2(i,j,iblk) = p5 - a11_3(i,j,iblk) = p5 - a11_4(i,j,iblk) = p5 + a11_1(i,j,iblk) = p5 + a11_2(i,j,iblk) = p5 + a11_3(i,j,iblk) = p5 + a11_4(i,j,iblk) = p5 + else + a11_1(i,j,iblk) = c0 + a11_2(i,j,iblk) = c0 + a11_3(i,j,iblk) = c0 + a11_4(i,j,iblk) = c0 + endif a12_1(i,j,iblk) = c0 a12_2(i,j,iblk) = c0 a12_3(i,j,iblk) = c0