Changeset 1175 for trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90

Timestamp:

2008-09-11T18:26:34+02:00 (16 years ago)

Author:

cetlod

Message:

update transport modules to take into account new trends organization, see ticket:248

File:

• trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90 (modified) (10 diffs)

trunk/NEMO/TOP_SRC/TRP/trcldf_lap.F90

@@ -1 +1 @@
 MODULE trcldf_lap
-   !!==============================================================================
+   !!======================================================================
    !!                       ***  MODULE  trcldf_lap  ***
    !! Ocean passive tracers:  horizontal component of the lateral tracer mixing trend
-   !!==============================================================================
+   !!======================================================================
+   !! History :  1.0  !  87-06  (P. Andrich, D. L Hostis)  Original code
+   !!                 !  91-11  (G. Madec)
+   !!                 !  95-02  (M. Levy)    passive tracers
+   !!                 !  95-11  (G. Madec)  suppress volumetric scale factors
+   !!                 !  96-01  (G. Madec)  statement function for e3
+   !!            8.5  !  02-06  (G. Madec)  F90: Free form and module
+   !!            9.0  !  04-03  (C. Ethe)   passive tracer
+   !!                 !  07-02  (C. Deltel)  Diagnose ML trends for passive tracers
+   !!----------------------------------------------------------------------
 #if defined key_top
-   !!----------------------------------------------------------------------
-   !!   'key_top'                                                TOP models
    !!----------------------------------------------------------------------
    !!   trc_ldf_lap  : update the tracer trend with the horizontal diffusion
@@ -14 +21 @@
    USE trp_trc                 ! ocean passive tracers variables
    USE prtctl_trc          ! Print control for debbuging
+   USE trdmld_trc
+   USE trdmld_trc_oce     
 
    IMPLICIT NONE
@@ -25 +34 @@
    !!----------------------------------------------------------------------
    !!   TOP 1.0 , LOCEAN-IPSL (2005) 
-   !! $Id$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/TRP/trcldf_lap.F90,v 1.10 2006/09/12 11:10:14 opalod Exp $ 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
    
@@ -53 +62 @@
       !! ** Action  : - Update tra arrays with the before iso-level 
       !!                harmonic mixing trend.
-      !!              - Save the trends in trtrd ('key_trc_diatrd')
-      !!
-      !! History :
-      !!   1.0  !  87-06  (P. Andrich, D. L Hostis)  Original code
-      !!        !  91-11  (G. Madec)
-      !!        !  95-02  (M. Levy)    passive tracers
-      !!        !  95-11  (G. Madec)  suppress volumetric scale factors
-      !!        !  96-01  (G. Madec)  statement function for e3
-      !!   8.5  !  02-06  (G. Madec)  F90: Free form and module
-      !!   9.0  !  04-03  (C. Ethe)   passive tracer
+      !!              - Save the trends ('key_trdmld_trc')
       !!----------------------------------------------------------------------
       USE oce_trc          , ztu => ua,  &  ! use ua as workspace
          &                   ztv => va      ! use va as workspace
 
-      !! * Arguments
       INTEGER, INTENT( in ) ::   kt       ! ocean time-step index
-      
+
       !! * Local save
       REAL(wp), DIMENSION(jpi,jpj), SAVE ::   &
          ze1ur, ze2vr, zbtr2              ! scale factor coefficients
       
-      !! * Local declarations
-      INTEGER ::   ji, jj, jk,jn         ! dummy loop indices
+      INTEGER ::   ji, jj, jk, jn         ! dummy loop indices
       REAL(wp) ::   &
-         zabe1, zabe2, zbtr              ! temporary scalars
+         zabe1, zabe2, zbtr               ! temporary scalars
 
       REAL(wp) ::   &
-         ztra, ztrax, ztray              ! workspace
+         ztra, ztrax, ztray               ! workspace
       CHARACTER (len=22) :: charout
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrtrd
       !!----------------------------------------------------------------------
       
@@ -92 +91 @@
          zbtr2(:,:) = 1. / ( e1t(:,:) * e2t(:,:) )
       ENDIF
- 
-      DO jn = 1, jptra
+
+      IF( l_trdtrc ) THEN
+!         STOP 'trcldf_lap: this was never validated, please comment this line to proceed...'
+      ENDIF
+
+      IF( l_trdtrc ) ALLOCATE( ztrtrd(jpi,jpj,jpk) )
+
+      !                                                          ! ===========
+      DO jn = 1, jptra                                           ! tracer loop
+         !                                                       ! ===========
+         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)
          
          !                                                  ! =============
@@ -99 +107 @@
             !                                               ! =============
             ! 1. First derivative (gradient)
-            ! -------------------
+            ! ------------------------------
             DO jj = 1, jpjm1
                DO ji = 1, fs_jpim1   ! vector opt.
@@ -116 +124 @@
 
             ! 2. Second derivative (divergence)
-            ! --------------------
+            ! ---------------------------------
             DO jj = 2, jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.
@@ -130 +138 @@
                   ! add it to the general tracer trends
                   tra(ji,jj,jk,jn) = tra(ji,jj,jk,jn) + ztrax + ztray
-
 #if defined key_trc_diatrd
                   ! save the horizontal diffusive trends
@@ -136 +143 @@
                   IF (luttrd(jn)) trtrd(ji,jj,jk,ikeep(jn),5) = ztray
 #endif
+
                END DO
             END DO
@@ -141 +149 @@
          END DO                                             !  End of slab  
          !                                                  ! =============
+         IF( l_trdtrc ) THEN
+            ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:)
+            IF (luttrd(jn)) CALL trd_mod_trc( ztrtrd, jn, jptrc_trd_ldf, kt )   ! trends diags
+         END IF
+         !                                                       ! ===========
+      END DO                                                     ! tracer loop
+      !                                                          ! ===========
 
-      END DO
+      IF( l_trdtrc ) DEALLOCATE( ztrtrd )
 
-     IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
+      IF( ln_ctl ) THEN    ! print mean trends (used for debugging)
          WRITE(charout, FMT="('ldf - lap')")
-         CALL prt_ctl_trc_info(charout)
-         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm,clinfo2='trd')
+         CALL prt_ctl_trc_info( charout )
+         CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' )
       ENDIF