Changeset 15071 for NEMO/trunk/src/OCE/ZDF/zdftke.F90

Timestamp:

2021-07-02T15:12:08+02:00 (3 years ago)

Author:

clem

Message:

make tke and gls work with tiling in debug mode. But it needs to be checked over

File:

• NEMO/trunk/src/OCE/ZDF/zdftke.F90 (modified) (2 diffs)

NEMO/trunk/src/OCE/ZDF/zdftke.F90

@@ -428 +428 @@
       !
       !                          !* Matrix inversion from level 2 (tke prescribed at level 1)
-      DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 )                ! First recurrence : Dk = Dk - Lk * Uk-1 / Dk-1
+      DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 )                ! First recurrence : Dk = Dk - Lk * Uk-1 / Dk-1
          zdiag(ji,jj,jk) = zdiag(ji,jj,jk) - zd_lw(ji,jj,jk) * zd_up(ji,jj,jk-1) / zdiag(ji,jj,jk-1)
       END_3D
@@ -435 +435 @@
 !         zd_lw(ji,jj,2) = en(ji,jj,2) - zd_lw(ji,jj,2) * en(ji,jj,1)    ! Surface boudary conditions on tke
 !      END_2D
-      DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 )
+      DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 2, jpkm1 )
          zd_lw(ji,jj,jk) = en(ji,jj,jk) - zd_lw(ji,jj,jk) / zdiag(ji,jj,jk-1) *zd_lw(ji,jj,jk-1)
       END_3D