Changeset 15055 for NEMO/trunk/src/OCE/DYN/divhor.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/DYN/divhor.F90 (modified) (3 diffs)

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

@@ -79 +79 @@
       ENDIF
       !
-      DO_3D_OVR( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpkm1 )                                    !==  Horizontal divergence  ==!
+      DO_3D_OVR( 0, 0, 0, 0, 1, jpkm1 )                                          !==  Horizontal divergence  ==!
          ! round brackets added to fix the order of floating point operations
          ! needed to ensure halo 1 - halo 2 compatibility
@@ -90 +90 @@
             &             )  * r1_e1e2t(ji,jj) / e3t(ji,jj,jk,Kmm)
       END_3D
+      CALL lbc_lnk( 'divhor', hdiv, 'T', 1.0_wp )   !   (no sign change)
+      !                                             ! needed for ww
       !
-      IF( ln_rnf )   CALL sbc_rnf_div( hdiv, Kmm )                     !==  runoffs    ==!   (update hdiv field)
+      IF( ln_rnf )   CALL sbc_rnf_div( hdiv, Kmm )                               !==  runoffs    ==!   (update hdiv field)
       !
 #if defined key_asminc 
@@ -99 +101 @@
       IF( ln_isf )                      CALL isf_hdiv( kt, Kmm, hdiv )           !==  ice shelf         ==!   (update hdiv field)
       !
-      IF (nn_hls==1) CALL lbc_lnk( 'divhor', hdiv, 'T', 1.0_wp )   !   (no sign change)
-      !
       IF( ln_timing )   CALL timing_stop('div_hor')
       !