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

trunk/NEMO/OPA_SRC/DYN/dynspg_fsc.F90

@@ -16 +16 @@
    USE oce             ! ocean dynamics and tracers 
    USE dom_oce         ! ocean space and time domain 
-   USE trdtra_oce      ! ocean active tracer   trend
-   USE trddyn_oce      ! ocean active dynamics trend
+   USE trdmod          ! ocean dynamics trends 
+   USE trdmod_oce      ! ocean variables trends
    USE zdf_oce         ! ocean vertical physics
    USE in_out_manager  ! I/O manager
    USE phycst          ! physical constants
    USE ocesbc          ! ocean surface boundary condition
-   USE flxrnf          ! ???
-   USE sol_oce         ! solver variables
+   USE flxrnf          ! ocean runoffs
+   USE sol_oce         ! ocean elliptic solver
    USE solpcg          ! preconditionned conjugate gradient solver
    USE solsor          ! Successive Over-relaxation solver
@@ -30 +30 @@
    USE obcdyn          ! ocean open boundary condition (obc_dyn routines)
    USE obcvol          ! ocean open boundary condition (obc_vol routines)
-   USE lib_mpp         ! ???
-   USE lbclnk          ! ???
-   USE cla_dynspg      ! ???
+   USE lib_mpp         ! distributed memory computing library
+   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
+   USE cla_dynspg      ! cross land advection
 
    IMPLICIT NONE
@@ -94 +94 @@
       !!
       !! ** Action : - Update (ua,va) with the surf. pressure gradient trend
-      !!             - Save the trends in (utrd,vtrd) ('key_diatrends')
+      !!             - Save the trends in (ztdua,ztdva) ('key_trddyn')
       !!
       !! References :
@@ -101 +101 @@
       !! History :
       !!        !  98-05 (G. Roullet)  Original code
-      !!        !  98-10  (G. Madec, M. Imbard)  release 8.2
-      !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
-      !!        !  02-11  (C. Talandier, A-M Treguier) Open boundaries
+      !!        !  98-10 (G. Madec, M. Imbard)  release 8.2
+      !!   8.5  !  02-08 (G. Madec)  F90: Free form and module
+      !!        !  02-11 (C. Talandier, A-M Treguier) Open boundaries
+      !!   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
@@ -155 +160 @@
             spgu(ji,jj) = zspgu
             spgv(ji,jj) = zspgv
-#if defined key_trddyn || defined key_trd_vor
-            ! save the surface pressure gradient trends 
-            utrd2(ji,jj,1) = zspgu
-            vtrd2(ji,jj,1) = zspgv
-#endif
          END DO 
       END DO 
@@ -170 +170 @@
                ua(ji,jj,jk) = ua(ji,jj,jk) + spgu(ji,jj)
                va(ji,jj,jk) = va(ji,jj,jk) + spgv(ji,jj)
-#if defined key_trddyn || defined key_trd_vor
-               ! save the surface pressure gradient trend for diagnostics
-               utrd(ji,jj,jk,8) = spgu(ji,jj)
-               vtrd(ji,jj,jk,8) = spgv(ji,jj) 
-#endif
-            END DO
-         END DO
-      END DO
+            END DO
+         END DO
+      END DO
+
+      ! Save the surface pressure gradient trend for diagnostics
+      IF( l_trddyn )   THEN
+         DO jk = 1, jpkm1
+            ztdua(:,:,jk) = spgu(:,:) 
+            ztdva(:,:,jk) = spgv(:,:) 
+         END DO
+      ENDIF
       
       ! 1. Evaluate the masked next velocity
@@ -310 +313 @@
             spgv(ji,jj) = z2dt * ztdgv
 #endif
-#if defined key_trddyn || defined key_trd_vor
-            ! save the transport divergence gradient trends
-            utrd2(ji,jj,2) = utrd2(ji,jj,2) + ztdgu
-            vtrd2(ji,jj,2) = vtrd2(ji,jj,2) + ztdgv
-#endif
          END DO
       END DO
@@ -328 +326 @@
                ua(ji,jj,jk) = (ua(ji,jj,jk) + spgu(ji,jj)) * umask(ji,jj,jk)
                va(ji,jj,jk) = (va(ji,jj,jk) + spgv(ji,jj)) * vmask(ji,jj,jk)
-#if defined key_trddyn || defined key_trd_vor
-               ! save the surface pressure gradient trend for diagnostics
-               utrd(ji,jj,jk,8) = utrd(ji,jj,jk,8) + spgu(ji,jj)/z2dt
-               vtrd(ji,jj,jk,8) = vtrd(ji,jj,jk,8) + spgv(ji,jj)/z2dt
-#endif
-            END DO
-         END DO
-      END DO
+            END DO
+         END DO
+      END DO
+
+      ! save the surface pressure gradient trends for diagnostic
+      ! momentum trends
+      IF( l_trddyn )   THEN
+         DO jk = 1, jpkm1
+            ztdua(:,:,jk) = ztdua(:,:,jk) + spgu(:,:)/z2dt 
+            ztdva(:,:,jk) = ztdva(:,:,jk) + spgv(:,:)/z2dt 
+         END DO
+
+         CALL trd_mod(ztdua, ztdva, jpdtdspg, 'DYN', kt)
+      ENDIF
 
       IF(l_ctl) THEN         ! print sum trends (used for debugging)