Changeset 216 for trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

Timestamp:

2005-03-17T15:02:38+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE151 : New trends organization

File:

• trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/DYN/dynldf_bilap.F90

@@ -14 +14 @@
    USE ldfdyn_oce      ! ocean dynamics: lateral physics
    USE in_out_manager  ! I/O manager
-   USE trddyn_oce      ! dynamics trends diagnostics variables
+   USE trdmod          ! ocean dynamics trends 
+   USE trdmod_oce      ! ocean variables trends
    USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 
@@ -66 +67 @@
       !! ** Action : - Update (ua,va) with the before iso-level biharmonic
       !!               mixing trend.
-      !!             - Save in (utrd,vtrd) the trends ('key_diatrends')
+      !!             - Save in (ztdua,ztdva) the trends ('key_trddyn')
       !!
       !! History :
@@ -75 +76 @@
       !!        !  97-07  (G. Madec)  lbc calls
       !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
-      !!----------------------------------------------------------------------
+      !!   9.0  !  04-08  (C. Talandier) New trends organization
+      !!----------------------------------------------------------------------
+      !! * Modules used     
+      USE oce, ONLY :    ztdua => ta,      & ! use ta as 3D workspace   
+                         ztdva => sa         ! use sa as 3D workspace   
+
       !! * Arguments
       INTEGER, INTENT( in ) ::   kt           ! ocean time-step index
@@ -97 +103 @@
       zlu(:,:) = 0.e0
       zlv(:,:) = 0.e0
+
+      ! Save ua and va trends
+      IF( l_trddyn )   THEN
+         ztdua(:,:,:) = ua(:,:,:) 
+         ztdva(:,:,:) = va(:,:,:) 
+      ENDIF
       !                                                ! ===============
       DO jk = 1, jpkm1                                 ! Horizontal slab
@@ -208 +220 @@
                ua(ji,jj,jk) = ua(ji,jj,jk) + zua
                va(ji,jj,jk) = va(ji,jj,jk) + zva
-#if defined key_trddyn || defined key_trd_vor
-               ! save the horizontal biharmonic diffusive trends
-               utrd(ji,jj,jk,5) = zua
-               vtrd(ji,jj,jk,5) = zva
-#endif
             END DO
          END DO
@@ -219 +226 @@
       END DO                                           !   End of slab
       !                                                ! ===============
+      ! save the lateral diffusion trends for diagnostic
+      ! momentum trends
+      IF( l_trddyn )   THEN
+         ztdua(:,:,:) = ua(:,:,:) - ztdua(:,:,:)
+         ztdva(:,:,:) = va(:,:,:) - ztdva(:,:,:)
+
+         CALL trd_mod(ztdua, ztdva, jpdtdldf, 'DYN', kt)
+      ENDIF
 
       IF(l_ctl) THEN         ! print sum trends (used for debugging)