Changeset 14834 for NEMO/trunk/src/OCE/DYN/divhor.F90

Timestamp:

2021-05-11T11:24:44+02:00 (3 years ago)

Author:

hadcv

Message:

#2600: Merge in dev_r14273_HPC-02_Daley_Tiling

File:

• NEMO/trunk/src/OCE/DYN/divhor.F90 (modified) (4 diffs)

NEMO/trunk/src/OCE/DYN/divhor.F90

@@ -64 +64 @@
       !
       INTEGER  ::   ji, jj, jk    ! dummy loop indices
-      REAL(wp) ::   zraur, zdep   ! local scalars
-      REAL(wp), DIMENSION(jpi,jpj) :: ztmp
       !!----------------------------------------------------------------------
       !
@@ -71 +69 @@
       !
       IF( kt == nit000 ) THEN
-         IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*) 'div_hor : horizontal velocity divergence '
-         IF(lwp) WRITE(numout,*) '~~~~~~~   '
-         hdiv(:,:,:) = 0._wp    ! initialize hdiv for the halos at the first time step
+         IF( .NOT. l_istiled .OR. ntile == 1 )  THEN                       ! Do only on the first tile
+            IF(lwp) WRITE(numout,*)
+            IF(lwp) WRITE(numout,*) 'div_hor : horizontal velocity divergence '
+            IF(lwp) WRITE(numout,*) '~~~~~~~   '
+         ENDIF
+         DO_3D_OVR( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpk )
+            hdiv(ji,jj,jk) = 0._wp    ! initialize hdiv for the halos at the first time step
+         END_3D
       ENDIF
       !
-      DO_3D( 0, 0, 0, 0, 1, jpkm1 )                                    !==  Horizontal divergence  ==!
+      DO_3D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpkm1 )                                    !==  Horizontal divergence  ==!
          ! round brackets added to fix the order of floating point operations
          ! needed to ensure halo 1 - halo 2 compatibility
@@ -84 +86 @@
             &                )                                                             & ! bracket for halo 1 - halo 2 compatibility
             &              + ( e1v(ji,jj  ) * e3v(ji,jj  ,jk,Kmm) * vv(ji,jj  ,jk,Kmm)     &
-            &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     & 
+            &                - e1v(ji,jj-1) * e3v(ji,jj-1,jk,Kmm) * vv(ji,jj-1,jk,Kmm)     &
             &                )                                                             & ! bracket for halo 1 - halo 2 compatibility
             &             )  * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm)
@@ -95 +97 @@
       ! 
 #endif
-      !
       IF( ln_isf )                      CALL isf_hdiv( kt, Kmm, hdiv )           !==  ice shelf         ==!   (update hdiv field)
       !
-      CALL lbc_lnk( 'divhor', hdiv, 'T', 1.0_wp )   !   (no sign change)
+      IF (nn_hls==1) CALL lbc_lnk( 'divhor', hdiv, 'T', 1.0_wp )   !   (no sign change)
       !
       IF( ln_timing )   CALL timing_stop('div_hor')