/[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.1 - (hide annotations) (download)
Fri Apr 29 18:47:02 2005 UTC (19 years, 1 month ago) by dimitri
Branch: MAIN
o added diagnostics to pkg/kpp, including computation of mixed layer
  depth based on a temperature/density criterion
o updated verification/natl_box to test the new pkg/kpp diagnostics

1 dimitri 1.1 C $Header: $
2     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     diagTitle = 'KPP mixing layer depth based on bulk Ri criterion'
77     diagUnits = 'm '
78     diagCode = 'SM P 1 '
79     CALL DIAGNOSTICS_ADD2LIST( diagNum,
80     I diagName, diagCode, diagUnits, diagTitle, myThid )
81    
82     diagName = 'KPPmld '
83     diagTitle = 'Mixed layer depth for dT=.8degC density criterion'
84     diagUnits = 'm '
85     diagCode = 'SM P 1 '
86     CALL DIAGNOSTICS_ADD2LIST( diagNum,
87     I diagName, diagCode, diagUnits, diagTitle, myThid )
88    
89     diagName = 'KPPfrac '
90     diagTitle = 'Short-wave flux fraction penetrating mixing layer'
91     diagUnits = ' '
92     diagCode = 'SM P 1 '
93     CALL DIAGNOSTICS_ADD2LIST( diagNum,
94     I diagName, diagCode, diagUnits, diagTitle, myThid )
95    
96     #endif /* ALLOW_DIAGNOSTICS */
97    
98     RETURN
99     END

  ViewVC Help
Powered by ViewVC 1.1.22