fix debug failure when grid_ice=C

* compiling in debug mode using -init=snan,arrays requires initialization of variables
DeniseWorthen · Jun 20, 2022 · 26498db · 26498db
1 parent 2d5487a
commit 26498db
Showing 1 changed file with 36 additions and 5 deletions.
diff --git a/cicecore/cicedynB/dynamics/ice_dyn_shared.F90 b/cicecore/cicedynB/dynamics/ice_dyn_shared.F90
@@ -197,7 +197,11 @@ subroutine init_dyn (dt)
           stressm_1, stressm_2, stressm_3, stressm_4, &
           stress12_1, stress12_2, stress12_3, stress12_4, &
           stresspT, stressmT, stress12T, &
-          stresspU, stressmU, stress12U
+          stresspU, stressmU, stress12U, &
+          taubxE, taubyE, strairxE, strairyE, strocnxE, strocnyE, &
+          strtltxE, strtltyE, strintxE, strintyE, iceEmask, fmE, TbE, &
+          taubxN, taubyN, strairxN, strairyN, strocnxN, strocnyN, &
+          strtltxN, strtltyN, strintxN, strintyN, iceNmask, fmN, TbN
       use ice_state, only: uvel, vvel, uvelE, vvelE, uvelN, vvelN, divu, shear
       use ice_grid, only: ULAT, NLAT, ELAT, tarea
 
@@ -244,10 +248,37 @@ subroutine init_dyn (dt)
          uvel(i,j,iblk) = c0    ! m/s
          vvel(i,j,iblk) = c0    ! m/s
          if (grid_ice == 'CD' .or. grid_ice == 'C') then ! extra velocity variables
-            uvelE(i,j,iblk) = c0
-            vvelE(i,j,iblk) = c0
-            uvelN(i,j,iblk) = c0
-            vvelN(i,j,iblk) = c0
+            uvelE  (i,j,iblk) = c0
+            vvelE  (i,j,iblk) = c0
+           taubxE  (i,j,iblk) = c0
+           taubyE  (i,j,iblk) = c0
+           strairxE(i,j,iblk) = c0
+           strairyE(i,j,iblk) = c0
+           strocnxE(i,j,iblk) = c0
+           strocnyE(i,j,iblk) = c0
+           strtltxE(i,j,iblk) = c0
+           strtltyE(i,j,iblk) = c0
+           strintxE(i,j,iblk) = c0
+           strintyE(i,j,iblk) = c0
+           iceEmask(i,j,iblk) = c0
+           fmE     (i,j,iblk) = c0
+           TbE     (i,j,iblk) = c0
+
+            uvelN  (i,j,iblk) = c0
+            vvelN  (i,j,iblk) = c0
+           taubxN  (i,j,iblk) = c0
+           taubyN  (i,j,iblk) = c0
+           strairxN(i,j,iblk) = c0
+           strairyN(i,j,iblk) = c0
+           strocnxN(i,j,iblk) = c0
+           strocnyN(i,j,iblk) = c0
+           strtltxN(i,j,iblk) = c0
+           strtltyN(i,j,iblk) = c0
+           strintxN(i,j,iblk) = c0
+           strintyN(i,j,iblk) = c0
+           iceNmask(i,j,iblk) = c0
+           fmN     (i,j,iblk) = c0
+           TbN     (i,j,iblk) = c0
          endif
 
          ! strain rates