Changeset 3318 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

Timestamp:

2012-02-25T16:50:01+01:00 (12 years ago)

Author:

gm

Message:

Ediag branche: #927 split TRA/DYN trd computation

File:

• branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90 (modified) (5 diffs)

branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

@@ -4 +4 @@
    !! Ocean passive tracers:  surface boundary condition
    !!======================================================================
-   !! History :  8.2  !  1998-10  (G. Madec, G. Roullet, M. Imbard)  Original code
-   !!            8.2  !  2001-02  (D. Ludicone)  sea ice and free surface
-   !!            8.5  !  2002-06  (G. Madec)  F90: Free form and module
-   !!            9.0  !  2004-03  (C. Ethe)  adapted for passive tracers
-   !!                 !  2006-08  (C. Deltel) Diagnose ML trends for passive tracers
+   !! History :  1.0  !  2004-03  (C. Ethe)  adapt trasbc to passive tracers
+   !!            2.0  !  2006-08  (C. Deltel) Diagnose ML trends for passive tracers
    !!==============================================================================
 #if defined key_top
@@ -19 +16 @@
    USE trc             ! ocean  passive tracers variables
    USE prtctl_trc      ! Print control for debbuging
-   USE trdmod_oce
-   USE trdtra
+   USE trd_oce         ! trends: ocean variables
+   USE trdtra          ! trends: tracer manager
 
    IMPLICIT NONE
@@ -44 +41 @@
       !!       trend of tracer equations.
       !!
-      !! ** Method :
-      !!      * concentration/dilution effect:
+      !! ** Method :    concentration/dilution effect:
       !!            The surface freshwater flux modify the ocean volume
       !!         and thus the concentration of a tracer as :
@@ -55 +51 @@
       !! ** Action  : - Update the 1st level of tra with the trend associated
       !!                with the tracer surface boundary condition 
-      !!
       !!----------------------------------------------------------------------
+      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index
       !
-      INTEGER, INTENT( in ) ::   kt          ! ocean time-step index
-      !
-      INTEGER  ::   ji, jj, jn           ! dummy loop indices
+      INTEGER  ::   ji, jj, jn     ! dummy loop indices
       REAL(wp) ::   zsrau, zse3t   ! temporary scalars
-      CHARACTER (len=22) :: charout
-      REAL(wp), POINTER, DIMENSION(:,:  ) :: zemps
-      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrtrd
+      CHARACTER (len=22) ::   charout
+      REAL(wp), POINTER, DIMENSION(:,:  ) ::   zemps    ! 2D workspace
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztrtrd   ! 3D workspace
       !!---------------------------------------------------------------------
       !
@@ -108 +102 @@
          IF( l_trdtrc ) THEN
             ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:)
-            CALL trd_tra( kt, 'TRC', jn, jptra_trd_nsr, ztrtrd )
+            CALL trd_tra( kt, 'TRC', jn, jptra_nsr, ztrtrd )
          END IF
          !                                                       ! ===========