Changeset 1120

Timestamp:

2008-06-20T17:21:25+02:00 (16 years ago)

Author:

cetlod

Message:

update trc_trp.F90 module, see ticket 144

File:

• trunk/NEMO/TOP_SRC/TRP/trp_trc.F90 (modified) (2 diffs)

trunk/NEMO/TOP_SRC/TRP/trp_trc.F90

@@ -14 +14 @@
    !!   'key_top'                                                TOP models
    !!----------------------------------------------------------------------
-
-   !! passive tracers number
-   USE par_trc, ONLY :   jptra    =>   jptra    !: number of passive tracers
+   USE par_oce
+   USE par_trc
 
    !! passive tracers fields 
@@ -26 +25 @@
    !! time step
    USE trc    , ONLY :   ndttrc   =>   ndttrc   !: frequency of step on passive tracers (NAMELIST)
+   USE trc    , ONLY :   ctrcnm   =>   ctrcnm   !: 
+   USE trc    , ONLY :   numnat   =>   numnat   !: namelist
+   USE trc    , ONLY :   nittrc000 =>  nittrc000
 
-   !! non-centered advection scheme (smolarkiewicz)
-   USE trc    , ONLY :   rtrn     =>   rtrn     !: value for truncation (NAMELIST)
-   USE trc    , ONLY :   ctrcnm   =>   ctrcnm   !: value for truncation (NAMELIST)
+   !! isopycnal scheme for passive tracers
+   !! ------------------------------------  
+   REAL(wp), PUBLIC ::   ahtrb0     !: background diffusivity coefficient for passive tracer (m2/s)
+   REAL(wp), PUBLIC ::   trcrat     !: ratio between passive and active tracer coeff for diffusion
+   REAL(wp), PUBLIC ::   ahtrc0     !: horizontal eddy diffusivity for passive tracers (m2/s)
+   REAL(wp), PUBLIC ::   aeivtr0    !: eddy induced velocity coefficient (m2/s)
 
-# if defined key_trc_diaadd
-   USE par_trc, ONLY :   jpdia2d  =>  jpdia2d   !: number of 2D passive tracers diag
-   USE par_trc, ONLY :   jpdia3d  =>  jpdia3d   !: number of 3D passive tracers diag
-   USE trc    , ONLY :   trc2d    =>   trc2d    !: additional 2D variable for ouputs
-   USE trc    , ONLY :   trc3d    =>   trc3d    !: additional 3D variable for ouputs
+   !! interpolated gradient
+   !!--------------------------------------------------  
+   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jptra) ::   gtru   !: horizontal gradient at u-points at bottom ocean level
+   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jptra) ::   gtrv   !: horizontal gradient at v-points at bottom ocean level
+
+   
+# if defined key_trcldf_eiv && defined key_diaeiv
+   !! The three component of the eddy induced velocity
+   !! --------------------------------------------------
+   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jpk) ::   u_trc_eiv   !: u-eiv (m/s)
+   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jpk) ::   v_trc_eiv   !: v-eiv (m/s)
+   REAL(wp), PUBLIC, DIMENSION (jpi,jpj,jpk) ::   w_trc_eiv   !: w-eiv (m/s)
 # endif
+
+# if defined key_trc_diatrd
+   
+   !!  non conservative trends (biological, ...)
+   !! --------------------------------------------------
+   LOGICAL, PUBLIC, DIMENSION (jptra) ::   luttrd   !: large trends diagnostic to write or not (namelist)
+   
+   !!  Advection-diffusion trends
+   REAL(wp), PUBLIC, DIMENSION(:,:,:,:,:), ALLOCATABLE ::   trtrd   !: trends of the tracer equations
+   
+   INTEGER, PUBLIC, DIMENSION(jptra) :: ikeep ! indice of tracer for which dyn trends are stored
+   INTEGER, PUBLIC                   :: nkeep ! number of tracers for which dyn trends are stored 
+   !                                          ! (used to allocate trtrd buffer)
+
+   !! netcdf files and index common
+   !! --------------------------------------------------
+   INTEGER , PUBLIC ::   nwritetrd   !: frequency of additional arrays outputs(namelist)
+   
+# endif 
 
 #else