Changeset 911 for trunk/NEMO/OPA_SRC/DYN/divcur.F90

Timestamp:

2008-04-28T11:31:32+02:00 (16 years ago)

Author:

ctlod

Message:

Implementation of the BDY package, see ticket: #126

File:

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

trunk/NEMO/OPA_SRC/DYN/divcur.F90

@@ -14 +14 @@
    USE in_out_manager  ! I/O manager
    USE obc_oce         ! ocean lateral open boundary condition
+   USE bdy_oce        ! Unstructured open boundaries variables
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 
@@ -77 +78 @@
       !!   8.2  !  00-03  (G. Madec)  no slip accurate
       !!   9.0  !  03-08  (G. Madec)  merged of cur and div, free form, F90
+      !!        !  05-01  (J. Chanut, A. Sellar) unstructured open boundaries
       !!----------------------------------------------------------------------
       !! * Arguments
@@ -134 +136 @@
          ENDIF
 #endif
+#endif         
+#if defined key_bdy || key_bdy_tides
+         ! unstructured open boundaries (div must be zero behind the open boundary)
+         DO jj = 1, jpj
+           DO ji = 1, jpi
+             hdivn(ji,jj,jk)=hdivn(ji,jj,jk)*bdytmask(ji,jj)
+           END DO
+         END DO
 #endif         
 #if defined key_agrif
@@ -291 +301 @@
       !!   8.1  !  97-08  (J.M. Molines)  Open boundaries
       !!   9.0  !  02-09  (G. Madec, E. Durand)  Free form, F90
+      !!        !  05-01  (J. Chanut) Unstructured open boundaries
       !!----------------------------------------------------------------------
       !! * Arguments
@@ -344 +355 @@
 #endif
 #endif         
+#if defined key_bdy || key_bdy_tides
+         ! unstructured open boundaries (div must be zero behind the open boundary)
+         DO jj = 1, jpj
+           DO ji = 1, jpi
+             hdivn(ji,jj,jk)=hdivn(ji,jj,jk)*bdytmask(ji,jj)
+           END DO
+         END DO
+#endif        
 #if defined key_agrif
          if ( .NOT. AGRIF_Root() ) then