initialized all eap history variables to zero on land (#632)

CICE-Consortium · Sep 15, 2021 · 8c63ea0 · 8c63ea0
1 parent 554a302
commit 8c63ea0
Showing 1 changed file with 11 additions and 4 deletions.
diff --git 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