Changeset 2587 for branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynzdf_exp_tam.F90

Timestamp:

2011-02-15T12:58:59+01:00 (13 years ago)

Author:

vidard

Message:

refer to ticket #798

File:

• branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynzdf_exp_tam.F90 (modified) (3 diffs)

branches/TAM_V3_0/NEMOTAM/OPATAM_SRC/DYN/dynzdf_exp_tam.F90

@@ -20 +20 @@
    !! * Modules used
 
-   USE in_out_manager  ! I/O manager
+   USE par_kind      , ONLY: & ! Precision variables
+      & wp
+   USE par_oce       , ONLY: & ! Ocean space and time domain variables
+      & jpi,                 &
+      & jpj,                 & 
+      & jpk,                 & 
+      & jpim1,               &
+      & jpjm1,               & 
+      & jpkm1
+   USE oce_tam       , ONLY: & ! ocean dynamics and tracers
+      & ub_tl,               &
+      & vb_tl,               &
+      & ua_tl,               &
+      & va_tl,               &
+      & ub_ad,               &
+      & vb_ad,               &
+      & ua_ad,               &
+      & va_ad
+   USE zdf_oce       , ONLY: & ! ocean vertical physics
+      & avmu,                &
+      & avmv,                &
+      & n_zdfexp
+   USE dom_oce       , ONLY: & ! ocean space and time domain
+#if defined key_zco
+      & e3t_0,               &
+      & e3w_0,               &
+#else
+      & e3u,                 &
+      & e3v,                 &
+      & e3uw,                &
+      & e3vw,                &
+#endif
+      & umask,               &
+      & vmask
+   USE phycst        , ONLY: & ! physical constants
+      & rau0
+   USE in_out_manager, ONLY: & ! I/O manager
+      & nit000,              &
+      & nitend,              &
+      & numout,              &
+      & lwp,                 &
+      & ctl_stop
    IMPLICIT NONE
    PRIVATE
@@ -59 +100 @@
       !! * Local declarations
       INTEGER ::   ji, jj, jk, jl                              ! dummy loop indices
-      REAL(wp) ::   zrau0r, zlavmr, zua, zva                   ! temporary scalars
-      REAL(wp), DIMENSION(jpi,jpk) ::   zwx, zwy, zwz, zww     ! temporary workspace arrays
+      REAL(wp) ::   zrau0r, zlavmr, zuatl, zvatl                   ! temporary scalars
+      REAL(wp), DIMENSION(jpi,jpk) ::   zwxtl, zwytl, zwztl, zwwtl     ! temporary workspace arrays
       !!----------------------------------------------------------------------
 
@@ -68 +109 @@
          IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~~~ '
       ENDIF
-
-      CALL ctl_stop ('dyn_zdf_exp_tan not available yet')
+      ! Local constant initialization
+      ! -----------------------------
+      zrau0r = 1. / rau0                                   ! inverse of the reference density
+      zlavmr = 1. / float( n_zdfexp )                      ! inverse of the number of sub time step
+
+      !                                                ! ===============
+      DO jj = 2, jpjm1                                 !  Vertical slab
+         !                                             ! ===============
+
+         ! Surface boundary condition
+         DO ji = 2, jpim1
+            zwytl(ji,1) = 0.0_wp
+            zwwtl(ji,1) = 0.0_wp
+         END DO  
+
+         ! Initialization of x, z and contingently trends array
+         DO jk = 1, jpk
+            DO ji = 2, jpim1
+               zwxtl(ji,jk) = ub_tl(ji,jj,jk)
+               zwztl(ji,jk) = vb_tl(ji,jj,jk)
+            END DO  
+         END DO  
+
+         ! Time splitting loop
+         DO jl = 1, n_zdfexp
+
+            ! First vertical derivative
+            DO jk = 2, jpk
+               DO ji = 2, jpim1
+                  zwytl(ji,jk) = avmu(ji,jj,jk) * ( zwxtl(ji,jk-1) - zwxtl(ji,jk) ) / fse3uw(ji,jj,jk) 
+                  zwwtl(ji,jk) = avmv(ji,jj,jk) * ( zwztl(ji,jk-1) - zwztl(ji,jk) ) / fse3vw(ji,jj,jk)
+               END DO  
+            END DO  
+
+            ! Second vertical derivative and trend estimation at kt+l*rdt/n_zdfexp
+            DO jk = 1, jpkm1
+               DO ji = 2, jpim1
+                  zuatl = zlavmr*( zwytl(ji,jk) - zwytl(ji,jk+1) ) / fse3u(ji,jj,jk)
+                  zvatl = zlavmr*( zwwtl(ji,jk) - zwwtl(ji,jk+1) ) / fse3v(ji,jj,jk)
+                  ua_tl(ji,jj,jk) = ua_tl(ji,jj,jk) + zuatl
+                  va_tl(ji,jj,jk) = va_tl(ji,jj,jk) + zvatl
+
+                  zwxtl(ji,jk) = zwxtl(ji,jk) + p2dt*zuatl*umask(ji,jj,jk)
+                  zwztl(ji,jk) = zwztl(ji,jk) + p2dt*zvatl*vmask(ji,jj,jk)
+               END DO  
+            END DO  
+
+         END DO  
+
+         !                                             ! ===============
+      END DO                                           !   End of slab
+      !                                                ! ===============
+
+
    END SUBROUTINE dyn_zdf_exp_tan
    SUBROUTINE dyn_zdf_exp_adj( kt, p2dt )