Changeset 3318 for branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/TOP_SRC/TRP/trdmld_trc.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/trdmld_trc.F90 (modified) (5 diffs)

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

@@ -4 +4 @@
    !! Ocean diagnostics:  mixed layer passive tracer trends 
    !!======================================================================
-   !! History :  9.0  !  06-08  (C. Deltel)  Original code (from trdmld.F90)
-   !!                 !  07-04  (C. Deltel)  Bug fix : add trcrad trends
-   !!                 !  07-06  (C. Deltel)  key_gyre : do not call lbc_lnk
+   !! History :  1.0  !  2006-08  (C. Deltel)  Original code (from trdmld.F90)
+   !!             -   !  2007-04  (C. Deltel)  Bug fix : add trcrad trends
+   !!             -   !  2007-06  (C. Deltel)  key_gyre : do not call lbc_lnk
    !!----------------------------------------------------------------------
 #if   defined key_top && ( defined key_trdmld_trc   ||   defined key_esopa )
@@ -23 +23 @@
    USE zdfddm  , ONLY : avs  !: salinity vertical diffusivity coeff. at w-point
 # endif
+   USE trdmld_trc_rst    ! restart for diagnosing the ML trends
    USE trcnam_trp        ! passive tracers transport namelist variables
    USE trdmod_trc_oce    ! definition of main arrays used for trends computations
-   USE in_out_manager    ! I/O manager
    USE dianam            ! build the name of file (routine)
    USE ldfslp            ! iso-neutral slopes 
+   USE in_out_manager    ! I/O manager
    USE ioipsl            ! NetCDF library
    USE lbclnk            ! ocean lateral boundary conditions (or mpp link)
    USE lib_mpp           ! MPP library
-   USE trdmld_trc_rst    ! restart for diagnosing the ML trends
    USE prtctl            ! print control
    USE sms_pisces        ! PISCES bio-model
    USE sms_lobster       ! LOBSTER bio-model
-
+   USE wrk_nemo          ! Memory Allocation
+   
    IMPLICIT NONE
    PRIVATE
@@ -46 +47 @@
    PUBLIC trd_mld_bio_zint
 
-   CHARACTER (LEN=40) ::  clhstnam                                ! name of the trends NetCDF file
-   INTEGER ::   nmoymltrd
+   CHARACTER (LEN=40) ::   clhstnam                                ! name of the trends NetCDF file
+   INTEGER            ::   nmoymltrd
    INTEGER, ALLOCATABLE, SAVE, DIMENSION(:) ::   ndextrd1
    INTEGER, DIMENSION(jptra) ::   nidtrd, nh_t
@@ -61 +62 @@
 
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:,:) ::  ztmltrd2   !
-#if defined key_lobster
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::  ztmltrdbio2  ! only needed for mean diagnostics in trd_mld_bio()
-#endif
+   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:)   ::  ztmltrdbio2  ! only needed for mean diagnostics in trd_mld_bio()
 
    !! * Substitutions
 #  include "top_substitute.h90"
+#  include "zdfddm_substitute.h90"
    !!----------------------------------------------------------------------
    !! NEMO/TOP 3.3 , NEMO Consortium (2010)
@@ -367 +367 @@
       !!       In IV), the appropriate trends are written in the trends NetCDF file.
       !!
-      !! References :
-      !!       - Vialard & al.
-      !!       - See NEMO documentation (in preparation)
-      !!----------------------------------------------------------------------
-      !
+      !! References: Vialard et al., 2001, JPO.
+      !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt   ! ocean time-step index
       !
-      INTEGER ::   ji, jj, jk, jl, ik, it, itmod, jn
-      REAL(wp) ::   zavt, zfn, zfn2
-      !
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltot             ! d(trc)/dt over the anlysis window (incl. Asselin)
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlres             ! residual = dh/dt entrainment term
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlatf             ! for storage only
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlrad             ! for storage only (for trb<0 corr in trcrad)
-      !
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltot2            ! -+
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlres2            !  | working arrays to diagnose the trends
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltrdm2           !  | associated with the time meaned ML
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlatf2            !  | passive tracers
-      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlrad2            !  | (-> for trb<0 corr in trcrad)
+      INTEGER ::   ji, jj, jk, jl      ! dummy loop indices
+      INTEGER ::   ik, it, itmod, jn   ! local integers
+      REAL(wp)::   zavt, zfn, zfn2    ! local scalars
+      !
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltot   ! d(trc)/dt over the anlysis window (incl. Asselin)
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlres   ! residual = dh/dt entrainment term
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlatf   ! for storage only
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlrad   ! for storage only (for trb<0 corr in trcrad)
+      !
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltot2   ! -+
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlres2   !  | working arrays to diagnose the trends
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmltrdm2  !  | associated with the time meaned ML
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlatf2   !  | passive tracers
+      REAL(wp), POINTER, DIMENSION(:,:,:) ::   ztmlrad2   !  | (-> for trb<0 corr in trcrad)
       !
       CHARACTER (LEN= 5) ::   clvar
 #if defined key_dimgout
-      INTEGER ::   iyear,imon,iday
+      INTEGER ::   iyear, imon, iday
       CHARACTER(LEN=80) ::   cltext, clmode
 #endif