/[MITgcm]/MITgcm_contrib/submesoscale/code/gmredi_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm_contrib/submesoscale/code/gmredi_diagnostics_init.F

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

revision 1.2 by dimitri, Fri May 30 22:13:42 2008 UTC revision 1.4 by zhc, Fri Mar 12 18:31:00 2010 UTC
# Line 1  Line 1 
   
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
# Line 35  CEOP Line 34  CEOP
34  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
35  C     === Local variables ===  C     === Local variables ===
36  C     diagNum   :: diagnostics number in the (long) list of available diag.  C     diagNum   :: diagnostics number in the (long) list of available diag.
37    C     diagMate  :: diag. mate number in the (long) list of available diag.
38  C     diagName  :: local short name (8c) of a diagnostics  C     diagName  :: local short name (8c) of a diagnostics
39  C     diagCode  :: local parser field with characteristics of the diagnostics  C     diagCode  :: local parser field with characteristics of the diagnostics
40  C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT  C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
41  C     diagUnits :: local string (16c): physical units of a diagnostic field  C     diagUnits :: local string (16c): physical units of a diagnostic field
42  C     diagTitle :: local string (80c): description of field in diagnostic  C     diagTitle :: local string (80c): description of field in diagnostic
43        INTEGER       diagNum        INTEGER       diagNum
44          INTEGER       diagMate
45        CHARACTER*8   diagName        CHARACTER*8   diagName
46        CHARACTER*16  diagCode        CHARACTER*16  diagCode
47        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
# Line 51  C---+----1----+----2----+----3----+----4 Line 52  C---+----1----+----2----+----3----+----4
52  c     IF ( useDiagnotics ) THEN  c     IF ( useDiagnotics ) THEN
53    
54          diagName  = 'GM_VisbK'          diagName  = 'GM_VisbK'
55          diagTitle =          diagTitle =
56       &     'Mixing coefficient from Visbeck etal parameterization'       &     'Mixing coefficient from Visbeck etal parameterization'
57          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
58          diagCode  = 'SM P    M1      '          diagCode  = 'SM P    M1      '
59          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
60       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
61    
62          diagName  = 'GM_hTrsL'          diagName  = 'GM_hTrsL'
63          diagTitle = 'Base depth (>0) of the Transition Layer'          diagTitle = 'Base depth (>0) of the Transition Layer'
64          diagUnits = 'm               '          diagUnits = 'm               '
65          diagCode  = 'SM P    M1      '          diagCode  = 'SM P    M1      '
66          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
67       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
68    
69          diagName  = 'GM_baseS'          diagName  = 'GM_baseS'
70          diagTitle = 'Slope at the base of the Transition Layer'          diagTitle = 'Slope at the base of the Transition Layer'
71          diagUnits = '1               '          diagUnits = '1               '
72          diagCode  = 'SM P    M1      '          diagCode  = 'SM P    M1      '
73          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
74       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
75    
76          diagName  = 'GM_rLamb'          diagName  = 'GM_rLamb'
77          diagTitle =          diagTitle =
78       &    'Slope vertical gradient at Trans. Layer Base (=recip.Lambda)'       &    'Slope vertical gradient at Trans. Layer Base (=recip.Lambda)'
79          diagUnits = '1/m             '          diagUnits = '1/m             '
80          diagCode  = 'SM P    M1      '          diagCode  = 'SM P    M1      '
81          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
82       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
83    
84          diagName  = 'GM_Kux  '          diagName  = 'GM_Kux  '
85          diagTitle = 'K_11 element (U.point, X.dir) of GM-Redi tensor'          diagTitle = 'K_11 element (U.point, X.dir) of GM-Redi tensor'
86          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
87          WRITE(diagCode,'(A,I3.3,A)') 'UU P ',diagNum+2,'MR      '          diagCode  = 'UU P    MR      '
88          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
89       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
90         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
91    
92          diagName  = 'GM_Kvy  '          diagName  = 'GM_Kvy  '
93          diagTitle = 'K_22 element (V.point, Y.dir) of GM-Redi tensor'          diagTitle = 'K_22 element (V.point, Y.dir) of GM-Redi tensor'
94          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
95          WRITE(diagCode,'(A,I3.3,A)') 'VV P ', diagNum ,'MR      '          diagCode  = 'VV P    MR      '
96          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
97       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
98         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
99    
100          diagName  = 'GM_Kuz  '          diagName  = 'GM_Kuz  '
101          diagTitle = 'K_13 element (U.point, Z.dir) of GM-Redi tensor'          diagTitle = 'K_13 element (U.point, Z.dir) of GM-Redi tensor'
102          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
103          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '          diagCode  = 'UU      MR      '
104          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
105       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
106         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
107    
108          diagName  = 'GM_Kvz  '          diagName  = 'GM_Kvz  '
109          diagTitle = 'K_23 element (V.point, Z.dir) of GM-Redi tensor'          diagTitle = 'K_23 element (V.point, Z.dir) of GM-Redi tensor'
110          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
111          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '          diagCode  = 'VV      MR      '
112          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
113       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
114         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
115    
116          diagName  = 'GM_Kwx  '          diagName  = 'GM_Kwx  '
117          diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'          diagTitle = 'K_31 element (W.point, X.dir) of GM-Redi tensor'
118          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
119          WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'LR      '          diagCode  = 'UM      LR      '
120          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
121       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
122         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
123    
124          diagName  = 'GM_Kwy  '          diagName  = 'GM_Kwy  '
125          diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'          diagTitle = 'K_32 element (W.point, Y.dir) of GM-Redi tensor'
126          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
127          WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum ,'LR      '          diagCode  = 'VM      LR      '
128          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
129       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
130         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
131    
132          diagName  = 'GM_Kwz  '          diagName  = 'GM_Kwz  '
133          diagTitle = 'K_33 element (W.point, Z.dir) of GM-Redi tensor'          diagTitle = 'K_33 element (W.point, Z.dir) of GM-Redi tensor'
134          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
135          diagCode  = 'WM P    LR      '          diagCode  = 'WM P    LR      '
136          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
137       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
138    
139          diagName  = 'GM_PsiX '          diagName  = 'GM_PsiX '
140          diagTitle = 'GM Bolus transport stream-function : X component'          diagTitle = 'GM Bolus transport stream-function : X component'
141          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
142          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'LR      '          diagCode  = 'UU      LR      '
143          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
144       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
145         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
146    
147          diagName  = 'GM_PsiY '          diagName  = 'GM_PsiY '
148          diagTitle = 'GM Bolus transport stream-function : Y component'          diagTitle = 'GM Bolus transport stream-function : Y component'
149          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
150          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'LR      '          diagCode  = 'VV      LR      '
151          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
152       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
153         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
154    
155          diagName  = 'GM_KuzTz'          diagName  = 'GM_KuzTz'
156          diagTitle = 'Redi Off-diagonal Temperature flux: X component'          diagTitle = 'Redi Off-diagonal Temperature flux: X component'
157          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
158          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '          diagCode  = 'UU      MR      '
159          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
160       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
161         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
162    
163          diagName  = 'GM_KvzTz'          diagName  = 'GM_KvzTz'
164          diagTitle = 'Redi Off-diagonal Temperature flux: Y component'          diagTitle = 'Redi Off-diagonal Temperature flux: Y component'
165          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
166          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '          diagCode  = 'VV      MR      '
167          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
168       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
169         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
170    
171          diagName  = 'GM_ubT  '          diagName  = 'GM_ubT  '
172          diagTitle = 'Zonal Mass-Weight Bolus Transp of Pot Temp'          diagTitle = 'Zonal Mass-Weight Bolus Transp of Pot Temp'
173          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
174          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '          diagCode  = 'UUr     MR      '
175          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
176       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
177         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
178    
179          diagName  = 'GM_vbT  '          diagName  = 'GM_vbT  '
180          diagTitle = 'Meridional Mass-Weight Bolus Transp of Pot Temp'          diagTitle = 'Meridional Mass-Weight Bolus Transp of Pot Temp'
181          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
182          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '          diagCode  = 'VVr     MR      '
183          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
184       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
185         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
186    
187  #ifdef GM_SUBMESO  #ifdef GM_SUBMESO
188  CBFK These are diagnosed whether or not Bolus Advection is used...  CBFK These are diagnosed whether or not Bolus Advection is used...
# Line 180  CBFK (i.e., GM diagnostic=GM+Submeso) Line 193  CBFK (i.e., GM diagnostic=GM+Submeso)
193          diagName  = 'SM_PsiX '          diagName  = 'SM_PsiX '
194          diagTitle = 'Submeso Bolus transport quasi-streamfunction : X'          diagTitle = 'Submeso Bolus transport quasi-streamfunction : X'
195          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
196          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'LR      '          diagCode  = 'UU      LR      '
197          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
198       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
199         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
200    
201          diagName  = 'SM_PsiY '          diagName  = 'SM_PsiY '
202          diagTitle = 'Submeso Bolus transport quasi-streamfunction : Y'          diagTitle = 'Submeso Bolus transport quasi-streamfunction : Y'
203          diagUnits = 'm^2/s           '          diagUnits = 'm^2/s           '
204          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'LR      '          diagCode  = 'VV      LR      '
205          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
206       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
207         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
208    
209          diagName  = 'SM_ubT  '          diagName  = 'SM_ubT  '
210          diagTitle = 'Zon. Submeso Mass-Weight Bolus Transp of Pot Temp'          diagTitle = 'Zon. Submeso Mass-Weight Bolus Transp of Pot Temp'
211          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
212          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '          diagCode  = 'UUr     MR      '
213          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
214       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
215         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
216    
217          diagName  = 'SM_vbT  '          diagName  = 'SM_vbT  '
218          diagTitle = 'Mer. Submeso Mass-Weight Bolus Transp of Pot Temp'          diagTitle = 'Mer. Submeso Mass-Weight Bolus Transp of Pot Temp'
219          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
220          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '          diagCode  = 'VVr     MR      '
221          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
222       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
223         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
224    
225          diagName  = 'SM_wbT  '          diagName  = 'SM_wbT  '
226          diagTitle = 'Rvel Submeso Mass-Weight Bolus Transp of Pot Temp'          diagTitle = 'Rvel Submeso Mass-Weight Bolus Transp of Pot Temp'
227          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
228          WRITE(diagCode,'(A,I3.3,A)') 'WM   ',diagNum,'LR      '          diagCode  = 'WM      LR      '
229          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
230       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
231         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
232    
233          diagName  = 'SM_KuzTz'          diagName  = 'SM_KuzTz'
234          diagTitle = 'Zon. Submeso Mass-Weight Kappa Transp of Pot Temp'          diagTitle = 'Zon. Submeso Mass-Weight Kappa Transp of Pot Temp'
235          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
236          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '          diagCode  = 'UU      MR      '
237          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
238       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
240    
241          diagName  = 'SM_KvzTz'          diagName  = 'SM_KvzTz'
242          diagTitle = 'Mer. Submeso Mass-Weight Kappa Transp of Pot Temp'          diagTitle = 'Mer. Submeso Mass-Weight Kappa Transp of Pot Temp'
243          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
244          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '          diagCode  = 'VV      MR      '
245          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
246       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
247         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
248    
249          diagName  = 'SM_KrddT'          diagName  = 'SM_KrddT'
250          diagTitle = 'Rvel Submeso Mass-Weight Kappa Transp of Pot Temp'          diagTitle = 'Rvel Submeso Mass-Weight Kappa Transp of Pot Temp'
251          diagUnits = 'degC.m^3/s      '          diagUnits = 'degC.m^3/s      '
252          WRITE(diagCode,'(A,I3.3,A)') 'WM   ',diagNum,'LR      '          diagCode  = 'WM      LR      '
253          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
254       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
255         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
256  #endif  #endif
257    
258  #ifdef ALLOW_TAU_EDDY  #ifdef ALLOW_EDDYPSI
259          diagName  = 'GMEdTauX'          diagName  = 'GMEdTauX'
260          diagTitle = 'eddy-induced stress X-comp. estimated from Kwx'          diagTitle = 'eddy-induced stress X-comp. estimated from Kwx'
261          diagUnits = 'N/m^2           '          diagUnits = 'N/m^2           '
262          WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'LR      '          diagCode  = 'UM      LR      '
263          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum + 2
264       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
265         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
266    
267          diagName  = 'GMEdTauY'          diagName  = 'GMEdTauY'
268          diagTitle = 'eddy-induced stress Y-comp. estimated from Kwy'          diagTitle = 'eddy-induced stress Y-comp. estimated from Kwy'
269          diagUnits = 'N/m^2           '          diagUnits = 'N/m^2           '
270          WRITE(diagCode,'(A,I3.3,A)') 'VM   ',diagNum+2,'LR      '          diagCode  = 'VM      LR      '
271          CALL DIAGNOSTICS_ADD2LIST( diagNum,          diagMate  = diagNum
272       I          diagName, diagCode, diagUnits, diagTitle, myThid )          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
273         I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
274  #endif  #endif
275    
276  c     ENDIF  c     ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22