Changeset 216 for trunk/NEMO/OPA_SRC/DYN/dynldf_lap.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_lap.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/DYN/dynldf_lap.F90

@@ -15 +15 @@
    USE zdf_oce         ! ocean vertical physics
    USE in_out_manager  ! I/O manager
-   USE trddyn_oce      ! ocean dynamics trends
+   USE trdmod          ! ocean dynamics trends 
+   USE trdmod_oce      ! ocean variables trends
    USE ldfslp          ! iso-neutral slopes 
 
@@ -56 +57 @@
       !! ** Action : - Update (ua,va) with the before iso-level harmonic 
       !!               mixing trend.
-      !!             - Save in (utrd,vtrd) arrays the trends ('key_diatrends')
+      !!             - Save in (ztdua,ztdva) arrays the trends ('key_trddyn')
       !!
       !! History :
@@ -62 +63 @@
       !!        !  91-11 (G. Madec)
       !!        !  96-01 (G. Madec) statement function for e3 and ahm
-      !!   8.5  !  02-06  (G. Madec)  F90: Free form and module
+      !!   8.5  !  02-06 (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
+      INTEGER, INTENT( in ) ::   kt       ! ocean time-step index
 
       !! * Local declarations
-      INTEGER  ::   ji, jj, jk           ! dummy loop indices
+      INTEGER  ::   ji, jj, jk            ! dummy loop indices
       REAL(wp) ::   &
-         zua, zva, ze2u, ze1v            ! temporary scalars
+         zua, zva, ze2u, ze1v             ! temporary scalars
       !!----------------------------------------------------------------------
 
@@ -77 +83 @@
          IF(lwp) WRITE(numout,*) 'dyn_ldf : iso-level harmonic (laplacien) operator'
          IF(lwp) WRITE(numout,*) '~~~~~~~ '
+      ENDIF
+
+      ! Save ua and va trends
+      IF( l_trddyn )   THEN
+         ztdua(:,:,:) = ua(:,:,:) 
+         ztdva(:,:,:) = va(:,:,:) 
       ENDIF
 
@@ -107 +119 @@
                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 diffusive trends
-               utrd(ji,jj,jk,5) = zua
-               vtrd(ji,jj,jk,5) = zva
-#endif
             END DO
          END DO
@@ -117 +124 @@
       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)