/[MITgcm]/MITgcm/pkg/kpp/kpp_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/kpp/kpp_diagnostics_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.6 - (hide annotations) (download)
Fri Sep 18 11:40:22 2009 UTC (14 years, 9 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint62, checkpoint61v, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.5: +17 -1 lines 
add double diffusive contributions as a hack, turned off by default
for now, and the code can be excluded with a CPP-flag
EXCLUDE_KPP_DOUBLEDIFF just as EXCLUDE_KPP_SHEAR_MIX
1 mlosch 1.6 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_diagnostics_init.F,v 1.5 2008/10/30 18:51:35 dfer Exp $
2 dimitri 1.1 C $Name: $
3    
4     #include "KPP_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: KPP_DIAGNOSTICS_INIT
8     C !INTERFACE:
9     SUBROUTINE KPP_DIAGNOSTICS_INIT( myThid )
10    
11     C !DESCRIPTION: \bv
12     C *==========================================================*
13     C | SUBROUTINE KPP_DIAGNOSTICS_INIT
14     C | o Routine to initialize list of all available diagnostics
15     C | for KPP package
16     C *==========================================================*
17     C \ev
18     C !USES:
19     IMPLICIT NONE
20    
21     C === Global variables ===
22     #include "EEPARAMS.h"
23    
24     C !INPUT/OUTPUT PARAMETERS:
25     C === Routine arguments ===
26     C myThid :: my Thread Id number
27     INTEGER myThid
28     CEOP
29    
30     #ifdef ALLOW_DIAGNOSTICS
31     C !LOCAL VARIABLES:
32     C === Local variables ===
33     C diagNum :: diagnostics number in the (long) list of available diag.
34     C diagName :: local short name (8c) of a diagnostics
35     C diagCode :: local parser field with characteristics of the diagnostics
36     C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
37     C diagUnits :: local string (16c): physical units of a diagnostic field
38     C diagTitle :: local string (80c): description of field in diagnostic
39     INTEGER diagNum
40     CHARACTER*8 diagName
41     CHARACTER*16 diagCode
42     CHARACTER*16 diagUnits
43     CHARACTER*(80) diagTitle
44    
45     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
46    
47     diagName = 'KPPviscA'
48     diagTitle = 'KPP vertical eddy viscosity coefficient'
49     diagUnits = 'm^2/s '
50     diagCode = 'SM P UR '
51     CALL DIAGNOSTICS_ADD2LIST( diagNum,
52     I diagName, diagCode, diagUnits, diagTitle, myThid )
53    
54     diagName = 'KPPdiffS'
55     diagTitle = 'Vertical diffusion coefficient for salt & tracers'
56     diagUnits = 'm^2/s '
57     diagCode = 'SM P UR '
58     CALL DIAGNOSTICS_ADD2LIST( diagNum,
59     I diagName, diagCode, diagUnits, diagTitle, myThid )
60    
61     diagName = 'KPPdiffT'
62     diagTitle = 'Vertical diffusion coefficient for heat'
63     diagUnits = 'm^2/s '
64     diagCode = 'SM P UR '
65     CALL DIAGNOSTICS_ADD2LIST( diagNum,
66     I diagName, diagCode, diagUnits, diagTitle, myThid )
67    
68     diagName = 'KPPghat '
69     diagTitle = 'Nonlocal transport coefficient'
70     diagUnits = 's/m^2 '
71     diagCode = 'SM P LR '
72     CALL DIAGNOSTICS_ADD2LIST( diagNum,
73     I diagName, diagCode, diagUnits, diagTitle, myThid )
74    
75     diagName = 'KPPhbl '
76 dimitri 1.2 diagTitle = 'KPP boundary layer depth, bulk Ri criterion'
77 dimitri 1.1 diagUnits = 'm '
78     diagCode = 'SM P 1 '
79     CALL DIAGNOSTICS_ADD2LIST( diagNum,
80     I diagName, diagCode, diagUnits, diagTitle, myThid )
81    
82     diagName = 'KPPfrac '
83     diagTitle = 'Short-wave flux fraction penetrating mixing layer'
84     diagUnits = ' '
85     diagCode = 'SM P 1 '
86     CALL DIAGNOSTICS_ADD2LIST( diagNum,
87     I diagName, diagCode, diagUnits, diagTitle, myThid )
88    
89 atn 1.3 diagName = 'KPPdbsfc'
90     diagTitle = 'Buoyancy difference with respect to surface'
91     diagUnits = 'm/s^2 '
92     diagCode = 'SM P UR '
93     CALL DIAGNOSTICS_ADD2LIST( diagNum,
94     I diagName, diagCode, diagUnits, diagTitle, myThid )
95    
96     diagName = 'KPPbfsfc'
97     diagTitle = 'Bo+radiation absorbed to d=hbf*hbl + plume'
98     diagUnits = 'm^2/s^3 '
99     diagCode = 'SM P UR '
100     CALL DIAGNOSTICS_ADD2LIST( diagNum,
101     I diagName, diagCode, diagUnits, diagTitle, myThid )
102    
103 dfer 1.4 diagName = 'KPPRi '
104     diagTitle = 'Bulk Richardson number'
105     diagUnits = 'non-dimensional '
106     diagCode = 'SM P UR '
107     CALL DIAGNOSTICS_ADD2LIST( diagNum,
108     I diagName, diagCode, diagUnits, diagTitle, myThid )
109    
110 atn 1.3 diagName = 'KPPbo '
111     diagTitle = 'Surface turbulent buoyancy forcing'
112     diagUnits = 'm^2/s^3 '
113     diagCode = 'SM P 1 '
114     CALL DIAGNOSTICS_ADD2LIST( diagNum,
115     I diagName, diagCode, diagUnits, diagTitle, myThid )
116    
117     diagName = 'KPPbosol'
118     diagTitle = 'surface radiative buoyancy forcing'
119     diagUnits = 'm^2/s^3 '
120     diagCode = 'SM P 1 '
121     CALL DIAGNOSTICS_ADD2LIST( diagNum,
122     I diagName, diagCode, diagUnits, diagTitle, myThid )
123    
124     diagName = 'KPPdbloc'
125     diagTitle = 'Local delta buoyancy across interfaces'
126     diagUnits = 'm/s^2 '
127     diagCode = 'SM P UR '
128     CALL DIAGNOSTICS_ADD2LIST( diagNum,
129     I diagName, diagCode, diagUnits, diagTitle, myThid )
130    
131 mlosch 1.6 #ifndef EXCLUDE_KPP_DOUBLEDIFF
132     diagName = 'KPPnuddt'
133     diagTitle = 'Vertical double diffusion coefficient for heat'
134     diagUnits = 'm^2/s '
135     diagCode = 'SM P UR '
136     CALL DIAGNOSTICS_ADD2LIST( diagNum,
137     I diagName, diagCode, diagUnits, diagTitle, myThid )
138    
139     diagName = 'KPPnudds'
140     diagTitle = 'Vertical double diffusion coefficient for salt'
141     diagUnits = 'm^2/s '
142     diagCode = 'SM P UR '
143     CALL DIAGNOSTICS_ADD2LIST( diagNum,
144     I diagName, diagCode, diagUnits, diagTitle, myThid )
145     #endif /* ndef EXCLUDE_KPP_DOUBLEDIFF */
146    
147 atn 1.3 #ifdef ALLOW_SALT_PLUME
148     diagName = 'KPPpfrac'
149     diagTitle = 'Salt plume flux fraction penetrating mixing layer'
150     diagUnits = ' '
151     diagCode = 'SM P 1 '
152     CALL DIAGNOSTICS_ADD2LIST( diagNum,
153     I diagName, diagCode, diagUnits, diagTitle, myThid )
154    
155     diagName = 'KPPboplm'
156     diagTitle = 'Surface haline buoyancy forcing'
157     diagUnits = 'm^2/s^3 '
158     diagCode = 'SM P 1 '
159     CALL DIAGNOSTICS_ADD2LIST( diagNum,
160     I diagName, diagCode, diagUnits, diagTitle, myThid )
161     #endif /* ALLOW_SALT_PLUME */
162    
163 dimitri 1.1 #endif /* ALLOW_DIAGNOSTICS */
164    
165     RETURN
166     END
  ViewVC Help
Powered by ViewVC 1.1.22