Changeset 15055 for NEMO/trunk/src/OCE/ISF/isfhdiv.F90

Timestamp:

2021-06-25T00:36:40+02:00 (3 years ago)

Author:

clem

Message:

nn_hls=2: make sure ww is defined on the halos since it is used in several places. So we need to ensure that hdiv is also defined on the halos

File:

• NEMO/trunk/src/OCE/ISF/isfhdiv.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/ISF/isfhdiv.F90

@@ -97 +97 @@
       !
       ! compute integrated divergence correction
-      DO_2D( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
+      DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
          zhdiv(ji,jj) = 0.5_wp * ( pfwf(ji,jj) + pfwf_b(ji,jj) ) * r1_rho0 / phtbl(ji,jj)
       END_2D
       !
       ! update divergence at each level affected by ice shelf top boundary layer
-      DO_2D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls )
+      DO_2D_OVR( nn_hls, nn_hls, nn_hls, nn_hls )
          ikt = ktop(ji,jj)
          ikb = kbot(ji,jj)
@@ -136 +136 @@
       !!----------------------------------------------------------------------
       !
-      DO_3D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpk )
+      DO_3D_OVR( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk )
          phdiv(ji,jj,jk) =  phdiv(ji,jj,jk) + pqvol(ji,jj,jk) * r1_e1e2t(ji,jj)   &
             &                             / e3t(ji,jj,jk,Kmm)