Fix diag logic

iangrooms · May 24, 2024 · 8f1c468 · 8f1c468
1 parent bc7aa4c
commit 8f1c468
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/parameterizations/stochastic/MOM_stochastics.F90 b/src/parameterizations/stochastic/MOM_stochastics.F90
@@ -179,9 +179,9 @@ subroutine stochastics_init(dt, grid, GV, CS, param_file, diag, Time)
 
   ! Initialize the "taper" fields. These fields multiply the components of the stochastic
   ! velocity increment in such a way as to smoothly taper them to zero at land boundaries.
-  if (CS%do_skeb) then
-    ALLOC_(CS%taperCu(grid%isdB:grid%iedB,grid%jsd:grid%jed))
-    ALLOC_(CS%taperCv(grid%isd:grid%ied,grid%jsdB:grid%jedB))
+  if ((CS%do_skeb) .or. (CS%id_skeb_taperu > 0) .or. (CS%id_skeb_taperv > 0)) then
+    ALLOC_(CS%taperCu(grid%IsdB:grid%IedB,grid%jsd:grid%jed))
+    ALLOC_(CS%taperCv(grid%isd:grid%ied,grid%JsdB:grid%JedB))
     ! Initialize taper from land mask
     do j=grid%jsd,grid%jed ; do I=grid%isdB,grid%iedB
       CS%taperCu(I,j) = grid%mask2dCu(I,j)