New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trdicp_oce.F90 in branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD – NEMO

source: branches/2012/dev_r3309_LOCEAN12_Ediag/NEMOGCM/NEMO/OPA_SRC/TRD/trdicp_oce.F90 @ 3316

Last change on this file since 3316 was 3316, checked in by gm, 12 years ago

Ediag branche: #927 add 3D output for dyn & tracer trends
  • Property svn:keywords set to Id
File size: 4.9 KB
Line 
1MODULE trdicp_oce
2   !!======================================================================
3   !!                   ***  MODULE trdicp_oce  ***
4   !! Ocean trends :   set tracer and momentum trend variables
5   !!======================================================================
6   !! History :  1.0  !  2004-08 (C. Talandier) New trends organization
7   !!            3.5  !  2012-02 (G. Madec)  add 3D tracer zdf trend output using iom
8   !!----------------------------------------------------------------------
9   USE par_oce                 ! ocean parameters
10
11   IMPLICIT NONE
12   PUBLIC
13
14#if  defined key_trdtra   &&   defined key_trddyn    ||   defined key_esopa
15   LOGICAL, PARAMETER ::   lk_trdtra = .TRUE.    !: tracers  trend flag
16   LOGICAL, PARAMETER ::   lk_trddyn = .TRUE.    !: momentum trend flag
17#elif  defined key_trdtra
18   LOGICAL, PARAMETER ::   lk_trdtra = .TRUE.    !: tracers  trend flag
19   LOGICAL, PARAMETER ::   lk_trddyn = .FALSE.   !: momentum trend flag
20#elif  defined key_trddyn
21   LOGICAL, PARAMETER ::   lk_trdtra = .FALSE.   !: tracers  trend flag
22   LOGICAL, PARAMETER ::   lk_trddyn = .TRUE.    !: momentum trend flag
23#else
24   LOGICAL, PARAMETER ::   lk_trdtra = .FALSE.   !: tracers  trend flag
25   LOGICAL, PARAMETER ::   lk_trddyn = .FALSE.   !: momentum trend flag
26#endif
27
28   !                                        !!! => tracer trends indexes <=
29   INTEGER, PARAMETER ::   jpicpt_xad =  1   !: x- horizontal advection
30   INTEGER, PARAMETER ::   jpicpt_yad =  2   !: y- horizontal advection
31   INTEGER, PARAMETER ::   jpicpt_zad =  3   !: z- vertical   advection
32   INTEGER, PARAMETER ::   jpicpt_ldf =  4   !: lateral       diffusion
33   INTEGER, PARAMETER ::   jpicpt_zdf =  5   !: vertical diffusion (Kz)
34   INTEGER, PARAMETER ::   jpicpt_bbc =  6   !: Bottom Boundary Condition (geoth. flux)
35   INTEGER, PARAMETER ::   jpicpt_bbl =  7   !: Bottom Boundary Layer (diffusive/convective)
36   INTEGER, PARAMETER ::   jpicpt_npc =  8   !: static instability mixing
37   INTEGER, PARAMETER ::   jpicpt_dmp =  9   !: damping
38   INTEGER, PARAMETER ::   jpicpt_qsr = 10   !: penetrative solar radiation
39   INTEGER, PARAMETER ::   jpicpt_nsr = 11   !: non solar radiation
40   INTEGER, PARAMETER ::   jpicpt_zl1 = 12   !: first level vertical flux
41
42   !                                        !!! => Total tracer trends indexes <=
43   INTEGER, PARAMETER ::   jptot_tra  = 12   !: change it when adding/removing one indice above
44   
45   !                                        !!! => dynamic trends indexes <=
46   INTEGER, PARAMETER ::   jpicpd_hpg =  1   !: hydrostatic pressure gradient
47   INTEGER, PARAMETER ::   jpicpd_keg =  2   !: kinetic energy gradient
48   INTEGER, PARAMETER ::   jpicpd_rvo =  3   !: relative vorticity
49   INTEGER, PARAMETER ::   jpicpd_pvo =  4   !: planetary vorticity
50   INTEGER, PARAMETER ::   jpicpd_ldf =  5   !: lateral diffusion
51   INTEGER, PARAMETER ::   jpicpd_had =  6   !: horizontal advection
52   INTEGER, PARAMETER ::   jpicpd_zad =  7   !: vertical advection
53   INTEGER, PARAMETER ::   jpicpd_zdf =  8   !: vertical diffusion
54   INTEGER, PARAMETER ::   jpicpd_spg =  9   !: surface pressure gradient
55   INTEGER, PARAMETER ::   jpicpd_dat = 10   !: damping term
56   INTEGER, PARAMETER ::   jpicpd_swf = 11   !: surface wind forcing
57   INTEGER, PARAMETER ::   jpicpd_bfr = 12   !: bottom friction
58
59   !                                        !!! => Total dynamic trends indexes <=
60   INTEGER, PARAMETER ::   jptot_dyn  = 12   !: change it when adding/removing one indice above
61   
62#if   defined key_trdtra   ||   defined key_trddyn   ||   defined key_esopa
63
64   !! Variables used for diagnostics
65   !!---------------------------------------------------------------------
66   REAL(wp) ::   tvolt        !: volume of the whole ocean computed at t-points
67   REAL(wp) ::   tvolu        !: volume of the whole ocean computed at u-points
68   REAL(wp) ::   tvolv        !: volume of the whole ocean computed at v-points
69
70   !! Tracers trends diagnostics variables
71   !!---------------------------------------------------------------------
72   REAL(wp), DIMENSION(jptot_tra) ::   tmo, smo         !: tracers trends average
73   REAL(wp), DIMENSION(jptot_tra) ::   t2, s2           !: tracers square trends average
74   
75   !! Momentum trends diagnostics variables
76   !!---------------------------------------------------------------------
77   REAL(wp), DIMENSION(jptot_dyn) ::   umo, vmo         !: momentum trends average
78   REAL(wp), DIMENSION(jptot_dyn) ::   hke              !: momentum square trends average
79   REAL(wp) ::   rpktrd   !: potential to kinetic energy conversion
80   REAL(wp) ::   peke     !: conversion potential energy - kinetic energy trend
81
82#endif
83   !!----------------------------------------------------------------------
84   !! NEMO/OPA 3.3 , NEMO Consortium (2010)
85   !! $Id$
86   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
87   !!======================================================================
88END MODULE trdicp_oce
Note: See TracBrowser for help on using the repository browser.