Changeset 15104 for NEMO/trunk/src/OCE/TRD/trdken.F90

Timestamp:

2021-07-07T16:36:00+02:00 (3 years ago)

Author:

clem

Message:

nn_hls=2: repair some loops.

File:

• NEMO/trunk/src/OCE/TRD/trdken.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/TRD/trdken.F90

@@ -100 +100 @@
       !
       zke(:,:,jpk) = 0._wp
-      zke(1,:, : ) = 0._wp
-      zke(:,1, : ) = 0._wp
-      DO_3D( 0, 1, 0, 1, 1, jpkm1 )
+      zke(1:nn_hls,:, : ) = 0._wp
+      zke(:,1:nn_hls, : ) = 0._wp
+      DO_3D( 0, nn_hls, 0, nn_hls, 1, jpkm1 )
          zke(ji,jj,jk) = 0.5_wp * rho0 *( uu(ji  ,jj,jk,Kmm) * putrd(ji  ,jj,jk) * bu(ji  ,jj,jk)  &
             &                           + uu(ji-1,jj,jk,Kmm) * putrd(ji-1,jj,jk) * bu(ji-1,jj,jk)  &
@@ -122 +122 @@
                            z2dx(:,:) = uu(:,:,1,Kmm) * ( utau_b(:,:) + utau(:,:) ) * e1e2u(:,:) * umask(:,:,1)
                            z2dy(:,:) = vv(:,:,1,Kmm) * ( vtau_b(:,:) + vtau(:,:) ) * e1e2v(:,:) * vmask(:,:,1)
-                           zke2d(1,:) = 0._wp   ;   zke2d(:,1) = 0._wp
-                           DO_2D( 0, 1, 0, 1 )
+                           zke2d(1:nn_hls,:) = 0._wp   ;   zke2d(:,1:nn_hls) = 0._wp
+                           DO_2D( 0, nn_hls, 0, nn_hls )
                               zke2d(ji,jj) = r1_rho0 * 0.5_wp * (   z2dx(ji,jj) + z2dx(ji-1,jj)   &
                               &                                   + z2dy(ji,jj) + z2dy(ji,jj-1)   ) * r1_bt(ji,jj,1)
@@ -219 +219 @@
 
       ! conv value on T-point
-      DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+      DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 )
          zcoef = 0.5_wp / e3t(ji,jj,jk,Kmm)
          pconv(ji,jj,jk) = zcoef * ( zconv(ji,jj,jk) + zconv(ji,jj,jk+1) ) * tmask(ji,jj,jk)