Changeset 2132 for branches/devmercator2010_1/NEMO/OPA_SRC/DYN/dynzdf_imp.F90

Timestamp:

2010-09-29T17:31:36+02:00 (14 years ago)

Author:

cbricaud

Message:

add change from DEV_r1784_GLS

File:

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

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

@@ -24 +24 @@
    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
@@ -75 +79 @@
       REAL(wp) ::   zzwi                               ! temporary scalars
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwi        ! temporary workspace arrays
+#if defined key_zdfgls
+      INTEGER :: ikbu, ikbv, ikbum1, ikbvm1
+      REAL(wp) :: zcbcu, zcbcv
+#endif
       !!----------------------------------------------------------------------
 
@@ -173 +181 @@
          END DO
       END DO
+
+#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
@@ -272 +292 @@
       END DO
 
+#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