1 
C $Header: /u/gcmpack/MITgcm/pkg/ggl90/ggl90_diagnostics_init.F,v 1.2 2010/01/03 18:48:53 jmc Exp $ 
2 
C $Name: $ 
3 

4 
#include "GGL90_OPTIONS.h" 
5 

6 
C+1+2+3+4+5+6+7+ 
7 
CBOP 0 
8 
C !ROUTINE: GGL90_DIAGNOSTICS_INIT 
9 

10 
C !INTERFACE: 
11 
SUBROUTINE GGL90_DIAGNOSTICS_INIT( myThid ) 
12 

13 
C !DESCRIPTION: 
14 
C Initialize list of all available diagnostics 
15 

16 
C !USES: 
17 
IMPLICIT NONE 
18 
#include "EEPARAMS.h" 
19 
#include "SIZE.h" 
20 
#include "GGL90.h" 
21 

22 
C !INPUT/OUTPUT PARAMETERS: 
23 
C myThid :: my Thread Id number 
24 
INTEGER myThid 
25 
CEOP 
26 

27 
#ifdef ALLOW_DIAGNOSTICS 
28 
C !LOCAL VARIABLES: 
29 
C === Local variables === 
30 
C diagNum :: diagnostics number in the (long) list of available diag. 
31 
C diagMate :: diag. mate number in the (long) list of available diag. 
32 
C diagName :: local short name (8c) of a diagnostics 
33 
C diagCode :: local parser field with characteristics of the diagnostics 
34 
C cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT 
35 
C diagUnits :: local string (16c): physical units of a diagnostic field 
36 
C diagTitle :: local string (80c): description of field in diagnostic 
37 
INTEGER diagNum 
38 
c INTEGER diagMate 
39 
CHARACTER*8 diagName 
40 
CHARACTER*16 diagCode 
41 
CHARACTER*16 diagUnits 
42 
CHARACTER*(80) diagTitle 
43 

44 
C+1+2+3+4+5+6+7+ 
45 

46 
diagName = 'GGL90TKE' 
47 
diagTitle = 'GGL90 subgrid turbulent kinetic energy' 
48 
diagUnits = 'm^2/s^2 ' 
49 
diagCode = 'SM LR ' 
50 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
51 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
52 

53 
diagName = 'GGL90Lmx' 
54 
diagTitle = 'Mixing length scale ' 
55 
diagUnits = 'm ' 
56 
diagCode = 'SM LR ' 
57 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
58 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
59 

60 
diagName = 'GGL90Prl' 
61 
diagTitle = 'Prandtl number used in GGL90' 
62 
diagUnits = '1 ' 
63 
diagCode = 'SM LR ' 
64 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
65 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
66 

67 
diagName = 'GGL90ArU' 
68 
diagTitle = 'GGL90 eddy viscosity at Upoint' 
69 
diagUnits = 'm^2/s ' 
70 
diagCode = 'SM LR ' 
71 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
72 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
73 

74 
diagName = 'GGL90ArV' 
75 
diagTitle = 'GGL90 eddy viscosity at Vpoint' 
76 
diagUnits = 'm^2/s ' 
77 
diagCode = 'SM LR ' 
78 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
79 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
80 

81 
diagName = 'GGL90Kr' 
82 
diagTitle = 'GGL90 diffusion coefficient for temperature' 
83 
diagUnits = 'm^2/s ' 
84 
diagCode = 'SM LR ' 
85 
CALL DIAGNOSTICS_ADDTOLIST( diagNum, 
86 
I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) 
87 

88 
#endif /* ALLOW_DIAGNOSTICS */ 
89 

90 
RETURN 
91 
END 