Changeset 2236 for branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

Timestamp:

2010-10-12T20:49:32+02:00 (14 years ago)

Author:

cetlod

Message:

First guess of NEMO_v3.3

File:

• branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynzdf_imp.F90 (modified) (4 diffs)

branches/DEV_r2106_LOCEAN2010/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

@@ -20 +20 @@
    USE phycst          ! physical constants
    USE in_out_manager  ! I/O manager
+#if defined key_zdfgls
+   USE zdfbfr, ONLY : bfrua, bfrva, wbotu, wbotv ! bottom stresses
+   USE lbclnk         ! ocean lateral boundary conditions (or mpp link)
+#endif
 
    IMPLICIT NONE
@@ -65 +69 @@
       REAL(wp) ::   zzwi, zzws, zrhs       ! temporary scalars
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwi        ! 3D workspace
+#if defined key_zdfgls
+      INTEGER :: ikbu, ikbv, ikbum1, ikbvm1
+      REAL(wp) :: zcbcu, zcbcv
+#endif
       !!----------------------------------------------------------------------
 
@@ -155 +163 @@
          END DO
       END DO
-      !
-      DO jk = 1, jpkm1        !==  Normalization to obtain the general momentum trend ua  ==
+
+#if defined key_zdfgls
+      ! Save bottom stress for next time step
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ikbu   = MIN( mbathy(ji+1,jj  ), mbathy(ji,jj) )
+            ikbum1 = MAX( ikbu-1, 1 )
+            wbotu(ji,jj) = bfrua(ji,jj) * ua(ji,jj,ikbum1) * umask(ji,jj,ikbum1)
+         END DO
+      END DO
+      CALL lbc_lnk( wbotu(:,:), 'U', -1. )
+#endif
+
+      ! Normalization to obtain the general momentum trend ua
+      DO jk = 1, jpkm1
          DO jj = 2, jpjm1   
             DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -243 +264 @@
          END DO
       END DO
-      !
-      DO jk = 1, jpkm1     !==  Normalization to obtain the general momentum trend va  ==
+
+#if defined key_zdfgls
+      ! Save bottom stress for next time step
+      DO jj = 2, jpjm1
+         DO ji = fs_2, fs_jpim1   ! vector opt.
+            ikbv   = MIN( mbathy(ji,jj+1), mbathy(ji,jj) )
+            ikbvm1 = MAX( ikbv-1, 1 )
+            wbotv(ji,jj) = bfrva(ji,jj) * va(ji,jj,ikbvm1) * vmask(ji,jj,ikbvm1)
+         END DO
+      END DO
+      CALL lbc_lnk( wbotv(:,:), 'V', -1. )
+#endif
+
+      ! Normalization to obtain the general momentum trend va
+      DO jk = 1, jpkm1
          DO jj = 2, jpjm1   
             DO ji = fs_2, fs_jpim1   ! vector opt.