/[MITgcm]/MITgcm/pkg/layers/layers_diagnostics_init.F
ViewVC logotype

Annotation of /MITgcm/pkg/layers/layers_diagnostics_init.F

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

Revision 1.4 - (hide annotations) (download)
Wed Oct 17 18:49:15 2012 UTC (11 years, 8 months ago) by rpa
Branch: MAIN
Changes since 1.3: +68 -1 lines 
Layers variable names updated and new diagnostics added
1 rpa 1.4 C $Header: /u/gcmpack/MITgcm/pkg/layers/layers_diagnostics_init.F,v 1.3 2012/09/19 22:46:56 gforget Exp $
2 gforget 1.1 C $Name: $
3    
4     #include "LAYERS_OPTIONS.h"
5    
6     C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7    
8     SUBROUTINE LAYERS_DIAGNOSTICS_INIT( myThid )
9    
10     C ===================================================================
11     C Initialize LAYERS variables that are kept fixed during the run.
12     C ===================================================================
13    
14     IMPLICIT NONE
15     #include "EEPARAMS.h"
16     #include "SIZE.h"
17     #include "PARAMS.h"
18     #include "GRID.h"
19     #ifdef ALLOW_LAYERS
20     # include "LAYERS_SIZE.h"
21     # include "LAYERS.h"
22     #endif
23     #ifdef ALLOW_DIAGNOSTICS
24     # include "DIAGNOSTICS_SIZE.h"
25     # include "DIAGNOSTICS.h"
26     #endif
27    
28     C INPUT/OUTPUT PARAMETERS:
29     C myThid :: my Thread Id number
30     INTEGER myThid
31    
32 gforget 1.2 INTEGER iLa
33 gforget 1.1 INTEGER diagNum
34     INTEGER diagMate
35     CHARACTER*8 diagName
36     CHARACTER*16 diagCode
37     CHARACTER*16 diagUnits
38     CHARACTER*(80) diagTitle
39    
40     #ifdef ALLOW_LAYERS
41     #ifdef ALLOW_DIAGNOSTICS
42     IF ( useDiagnostics ) THEN
43    
44     diagNum=ndiagt
45 gforget 1.2
46     DO iLa=1,layers_maxNum
47     IF ( layers_num(iLa).NE.0 ) THEN
48    
49 rpa 1.4 C --- UH, UV
50 gforget 1.1 #ifdef LAYERS_UFLUX
51 gforget 1.2 WRITE(diagName,'(A4,I1,A3)') 'LaUH',iLa,layers_name(iLa)
52 gforget 1.1 diagTitle = 'Layer Integrated U (UH in m2/s)'
53     diagUnits = 'm2/s '
54     diagCode = 'UU MX '
55 gforget 1.3 #ifdef LAYERS_VFLUX
56 gforget 1.1 diagMate = diagNum + 2
57 gforget 1.3 #else
58     diagMate = 0
59     #endif
60 gforget 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
61     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
62     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
63     #endif
64    
65     #ifdef LAYERS_VFLUX
66 gforget 1.2 WRITE(diagName,'(A4,I1,A3)') 'LaVH',iLa,layers_name(iLa)
67 gforget 1.1 diagTitle = 'Layer Integrated V (VH m2/s)'
68     diagUnits = 'm2/s '
69     diagCode = 'VV MX '
70 gforget 1.3 #ifdef LAYERS_UFLUX
71 gforget 1.1 diagMate = diagNum
72 gforget 1.3 #else
73     diagMate = 0
74     #endif
75 gforget 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
76     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
77     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
78     #endif
79    
80 rpa 1.4 C --- Thickness
81 gforget 1.1 #ifdef LAYERS_THICKNESS
82     #ifdef LAYERS_UFLUX
83 gforget 1.2 WRITE(diagName,'(A4,I1,A3)') 'LaHw',iLa,layers_name(iLa)
84 gforget 1.1 diagTitle = 'Layer Thickness at U points (m)'
85     diagUnits = 'm '
86     diagCode = 'UU MX '
87 gforget 1.3 #ifdef LAYERS_VFLUX
88 gforget 1.1 diagMate = diagNum + 2
89 gforget 1.3 #else
90     diagMate = 0
91     #endif
92 gforget 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
93     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
94     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
95     #endif
96    
97     #ifdef LAYERS_VFLUX
98 gforget 1.2 WRITE(diagName,'(A4,I1,A3)') 'LaHs',iLa,layers_name(iLa)
99 gforget 1.1 diagTitle = 'Layer Thickness at V points (m)'
100     diagUnits = 'm '
101     diagCode = 'VV MX '
102 gforget 1.3 #ifdef LAYERS_UFLUX
103 gforget 1.1 diagMate = diagNum
104 gforget 1.3 #else
105     diagMate = 0
106     #endif
107 gforget 1.1 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
108     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
109     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
110     #endif
111 rpa 1.4
112     C --- PI
113     #ifdef LAYERS_THICKNESS
114     #ifdef LAYERS_UFLUX
115     WRITE(diagName,'(A4,I1,A3)') 'LaPw',iLa,layers_name(iLa)
116     diagTitle = 'Layer Probability at U points (m)'
117     diagUnits = ' '
118     diagCode = 'UU MX '
119     #ifdef LAYERS_VFLUX
120     diagMate = diagNum + 2
121     #else
122     diagMate = 0
123     #endif
124     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
125     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
126     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
127     #endif
128    
129     #ifdef LAYERS_VFLUX
130     WRITE(diagName,'(A4,I1,A3)') 'LaPs',iLa,layers_name(iLa)
131     diagTitle = 'Layer Probability at V points (m)'
132     diagUnits = ' '
133     diagCode = 'VV MX '
134     #ifdef LAYERS_UFLUX
135     diagMate = diagNum
136     #else
137     diagMate = 0
138     #endif
139     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
140     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
141     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
142     #endif
143    
144     C --- U, V
145     #ifdef LAYERS_THICKNESS
146     #ifdef LAYERS_UFLUX
147     WRITE(diagName,'(A4,I1,A3)') 'LaUa',iLa,layers_name(iLa)
148     diagTitle = 'Layer-averaged U velocity (non-weighted) (m/s)'
149     diagUnits = 'm/s '
150     diagCode = 'UU MX '
151     #ifdef LAYERS_VFLUX
152     diagMate = diagNum + 2
153     #else
154     diagMate = 0
155     #endif
156     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
157     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
158     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
159     #endif
160    
161     #ifdef LAYERS_VFLUX
162     WRITE(diagName,'(A4,I1,A3)') 'LaVa',iLa,layers_name(iLa)
163     diagTitle = 'Layer-averaged V velocity (non-weighted) (m/s)'
164     diagUnits = 'm/s '
165     diagCode = 'VV MX '
166     #ifdef LAYERS_UFLUX
167     diagMate = diagNum
168     #else
169     diagMate = 0
170     #endif
171     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
172     I diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
173     CALL DIAGNOSTICS_SETKLEV( diagName, Nlayers , myThid )
174     #endif
175    
176 gforget 1.1 #endif
177    
178     ENDIF
179 gforget 1.2 ENDDO
180    
181     ENDIF
182 gforget 1.1 #endif /* ALLOW_DIAGNOSTICS */
183     #endif /* ALLOW_LAYERS */
184    
185     RETURN
186     END
  ViewVC Help
Powered by ViewVC 1.1.22