Changeset 64

Timestamp:

12/27/12 21:55:06 (11 years ago)

Author:

smasson

Message:

add back missing part of NEMO 3.4: div must be zero behind the open boundary

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90 (modified) (2 diffs)

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/divcur.F90

@@ -27 +27 @@
    USE sbc_oce, ONLY : ln_rnf   ! surface boundary condition: ocean
    USE sbcrnf          ! river runoff 
+   USE obc_oce         ! ocean lateral open boundary condition
    USE cla             ! cross land advection             (cla_div routine)
    USE in_out_manager  ! I/O manager
@@ -302 +303 @@
          END DO  
 
+#if defined key_obc
+         IF( Agrif_Root() ) THEN
+            ! open boundaries (div must be zero behind the open boundary)
+            !  mpp remark: The zeroing of hdivn can probably be extended to 1->jpi/jpj for the correct row/column
+            IF( lp_obc_east  )   hdivn(nie0p1:nie1p1,nje0  :nje1  ,jk) = 0.e0      ! east
+            IF( lp_obc_west  )   hdivn(niw0  :niw1  ,njw0  :njw1  ,jk) = 0.e0      ! west
+            IF( lp_obc_north )   hdivn(nin0  :nin1  ,njn0p1:njn1p1,jk) = 0.e0      ! north
+            IF( lp_obc_south )   hdivn(nis0  :nis1  ,njs0  :njs1  ,jk) = 0.e0      ! south
+         ENDIF
+#endif 
+
          IF( .NOT. AGRIF_Root() ) THEN
             IF ((nbondi ==  1).OR.(nbondi == 2)) hdivn(nlci-1 , :     ,jk) = 0.e0      ! east