/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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

revision 1.5 by dimitri, Fri Dec 17 19:17:57 2004 UTC revision 1.11 by dimitri, Thu May 5 16:41:54 2005 UTC
# Line 24  C       \begin{tabular}[h]{|c|c|}\hline Line 24  C       \begin{tabular}[h]{|c|c|}\hline
24  C         \textbf{Positions}  &  \textbf{Characters}  C         \textbf{Positions}  &  \textbf{Characters}
25  C         &  \textbf{Meanings} \\\hline  C         &  \textbf{Meanings} \\\hline
26  C         parse(1)  &  S  &  scalar \\  C         parse(1)  &  S  &  scalar \\
27  C                   &  U  &  U-vector component  \\  C                   &  U  &  vector component in X direction \\
28  C                   &  V  &  V-vector component  \\  C                   &  V  &  vector component in Y direction \\
29    C                   &  W  &  vector component in vertical direction \\
30  C         parse(2)  &  U  &  C-grid U-Point  \\  C         parse(2)  &  U  &  C-grid U-Point  \\
31  C                   &  V  &  C-grid V-Point  \\  C                   &  V  &  C-grid V-Point  \\
32  C                   &  M  &  C-grid Mass Point  \\  C                   &  M  &  C-grid Mass Point  \\
# Line 63  C     !LOCAL VARIABLES: Line 64  C     !LOCAL VARIABLES:
64        CHARACTER*16  diagCode        CHARACTER*16  diagCode
65        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
66        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
67          CHARACTER*2   rUnit2c
68          CHARACTER*4   tUnit4c, sUnit4c
69    
70          CHARACTER*(16) DIAGS_MK_UNITS
71          EXTERNAL DIAGS_MK_UNITS
72    
73  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
74  C     For each output variable,  C     For each output variable,
# Line 70  C     specify Name (cdiag, 8c), Descript Line 76  C     specify Name (cdiag, 8c), Descript
76  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)  C         and Type/Parms (location on C grid, 2D/3D, ...) (gdiag, 16c)
77  C----------------------------------------------------------------------  C----------------------------------------------------------------------
78    
79          IF ( usingPCoords ) THEN
80            rUnit2c= 'Pa'
81          ELSE
82            rUnit2c= 'm '
83          ENDIF
84          IF ( fluidIsAir ) THEN
85            tUnit4c= 'K   '
86            sUnit4c= 'g/kg'
87          ELSE
88            tUnit4c= 'degC'
89            sUnit4c= 'psu '
90          ENDIF
91    
92    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
93    C-    state variables of the main code (and related quadratic var):
94    
95          diagName  = 'ETAN    '
96          diagTitle = 'Perturbation of Surface (pressure, height) (Pa,m)'
97          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
98          diagCode  = 'SM      M1      '
99          CALL DIAGNOSTICS_ADD2LIST( diagNum,
100         I          diagName, diagCode, diagUnits, diagTitle, myThid )
101    
102          diagName  = 'ETANSQ  '
103          diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)'
104          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
105          diagCode  = 'SM      M1      '
106          CALL DIAGNOSTICS_ADD2LIST( diagNum,
107         I          diagName, diagCode, diagUnits, diagTitle, myThid )
108    
109          diagName  = 'DETADT2 '
110          diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'
111          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
112          diagCode  = 'SM      M1      '
113          CALL DIAGNOSTICS_ADD2LIST( diagNum,
114         I          diagName, diagCode, diagUnits, diagTitle, myThid )
115    
116          diagName  = 'THETA   '
117          diagTitle = 'Potential Temperature (degC,K)'
118          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
119          diagCode  = 'SM      MR      '
120          CALL DIAGNOSTICS_ADD2LIST( diagNum,
121         I          diagName, diagCode, diagUnits, diagTitle, myThid )
122    
123          diagName  = 'SALT    '
124          diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
125          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
126          diagCode  = 'SM      MR      '
127          CALL DIAGNOSTICS_ADD2LIST( diagNum,
128         I          diagName, diagCode, diagUnits, diagTitle, myThid )
129    
130          diagName  = 'SALTanom'
131          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
132          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
133          diagCode  = 'SM      MR      '
134          CALL DIAGNOSTICS_ADD2LIST( diagNum,
135         I          diagName, diagCode, diagUnits, diagTitle, myThid )
136    
137          diagName  = 'UVEL    '
138          diagTitle = 'Zonal Component of Velocity (m/s)'
139          diagUnits = 'm/s             '
140          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
141          CALL DIAGNOSTICS_ADD2LIST( diagNum,
142         I          diagName, diagCode, diagUnits, diagTitle, myThid )
143    
144          diagName  = 'VVEL    '
145          diagTitle = 'Meridional Component of Velocity (m/s)'
146          diagUnits = 'm/s             '
147          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
148          CALL DIAGNOSTICS_ADD2LIST( diagNum,
149         I          diagName, diagCode, diagUnits, diagTitle, myThid )
150    
151          diagName  = 'WVEL    '
152          diagTitle = 'Vertical Component of Velocity (r_units/s)'
153          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
154          diagCode  = 'WM      LR      '
155          CALL DIAGNOSTICS_ADD2LIST( diagNum,
156         I          diagName, diagCode, diagUnits, diagTitle, myThid )
157    
158          diagName  = 'THETASQ '
159          diagTitle = 'Square of Potential Temperature (K^2)'
160          diagUnits = DIAGS_MK_UNITS( tUnit4c//'^2', myThid )
161          diagCode  = 'SM      MR      '
162          CALL DIAGNOSTICS_ADD2LIST( diagNum,
163         I          diagName, diagCode, diagUnits, diagTitle, myThid )
164    
165          diagName  = 'SALTSQ  '
166          diagTitle = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
167          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
168          diagCode  = 'SM      MR      '
169          CALL DIAGNOSTICS_ADD2LIST( diagNum,
170         I          diagName, diagCode, diagUnits, diagTitle, myThid )
171    
172          diagName  = 'SALTSQan'
173          diagTitle = 'Square of Salt anomaly (=SALTSQ-35^2 (g^2/kg^2)'
174          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
175          diagCode  = 'SM      MR      '
176          CALL DIAGNOSTICS_ADD2LIST( diagNum,
177         I          diagName, diagCode, diagUnits, diagTitle, myThid )
178    
179          diagName  = 'UVELSQ  '
180          diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
181          diagUnits = 'm^2/s^2         '
182          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
183          CALL DIAGNOSTICS_ADD2LIST( diagNum,
184         I          diagName, diagCode, diagUnits, diagTitle, myThid )
185    
186          diagName  = 'VVELSQ  '
187          diagTitle = 'Square of Meridional Comp of Velocity (m^2/s^2)'
188          diagUnits = 'm^2/s^2         '
189          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
190          CALL DIAGNOSTICS_ADD2LIST( diagNum,
191         I          diagName, diagCode, diagUnits, diagTitle, myThid )
192    
193          diagName  = 'WVELSQ  '
194          diagTitle = 'Square of Vertical Comp of Velocity (m^2/s^2)'
195          diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2/s^2', myThid )
196          diagCode  = 'WM      LR      '
197          CALL DIAGNOSTICS_ADD2LIST( diagNum,
198         I          diagName, diagCode, diagUnits, diagTitle, myThid )
199    
200          diagName  = 'UV_VEL_C'
201          diagTitle ='Product of horizontal Comp of velocity (cell center)'
202          diagUnits = 'm^2/s^2         '
203          WRITE(diagCode,'(A,I3.3,A)') 'UM   ', diagNum+1 ,'MR      '
204          CALL DIAGNOSTICS_ADD2LIST( diagNum,
205         I          diagName, diagCode, diagUnits, diagTitle, myThid )
206    
207          diagName  = 'UV_VEL_Z'
208          diagTitle = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
209          diagUnits = 'm^2/s^2         '
210          WRITE(diagCode,'(A,I3.3,A)') 'UZ   ', diagNum+1 ,'MR      '
211          CALL DIAGNOSTICS_ADD2LIST( diagNum,
212         I          diagName, diagCode, diagUnits, diagTitle, myThid )
213    
214          diagName  = 'WU_VEL  '
215          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
216          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
217          diagCode  = 'WU      LR      '
218          CALL DIAGNOSTICS_ADD2LIST( diagNum,
219         I          diagName, diagCode, diagUnits, diagTitle, myThid )
220    
221          diagName  = 'WV_VEL  '
222          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
223          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
224          diagCode  = 'WV      LR      '
225          CALL DIAGNOSTICS_ADD2LIST( diagNum,
226         I          diagName, diagCode, diagUnits, diagTitle, myThid )
227    
228          diagName  = 'UVELMASS'
229          diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
230          diagUnits = 'm/s             '
231          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
232          CALL DIAGNOSTICS_ADD2LIST( diagNum,
233         I          diagName, diagCode, diagUnits, diagTitle, myThid )
234    
235          diagName  = 'VVELMASS'
236          diagTitle = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
237          diagUnits = 'm/s             '
238          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
239          CALL DIAGNOSTICS_ADD2LIST( diagNum,
240         I          diagName, diagCode, diagUnits, diagTitle, myThid )
241    
242          diagName  = 'WVELMASS'
243          diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
244          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
245          diagCode  = 'WM      LR      '
246          CALL DIAGNOSTICS_ADD2LIST( diagNum,
247         I          diagName, diagCode, diagUnits, diagTitle, myThid )
248    
249          diagName  = 'UTHMASS '
250          diagTitle = 'Zonal Mass-Weight Transp of Pot Temp (K.m/s)'
251          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
252          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
253          CALL DIAGNOSTICS_ADD2LIST( diagNum,
254         I          diagName, diagCode, diagUnits, diagTitle, myThid )
255    
256          diagName  = 'VTHMASS '
257          diagTitle = 'Meridional Mass-Weight Transp of Pot Temp (K.m/s)'
258          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
259          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
260          CALL DIAGNOSTICS_ADD2LIST( diagNum,
261         I          diagName, diagCode, diagUnits, diagTitle, myThid )
262    
263          diagName  = 'WTHMASS '
264          diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
265          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
266          diagCode  = 'WM      LR      '
267          CALL DIAGNOSTICS_ADD2LIST( diagNum,
268         I          diagName, diagCode, diagUnits, diagTitle, myThid )
269    
270          diagName  = 'USLTMASS'
271          diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'
272          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
273          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
274          CALL DIAGNOSTICS_ADD2LIST( diagNum,
275         I          diagName, diagCode, diagUnits, diagTitle, myThid )
276    
277          diagName  = 'VSLTMASS'
278          diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'
279          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
280          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
281          CALL DIAGNOSTICS_ADD2LIST( diagNum,
282         I          diagName, diagCode, diagUnits, diagTitle, myThid )
283    
284          diagName  = 'WSLTMASS'
285          diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
286          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
287          diagCode  = 'WM      LR      '
288          CALL DIAGNOSTICS_ADD2LIST( diagNum,
289         I          diagName, diagCode, diagUnits, diagTitle, myThid )
290    
291          diagName  = 'UVELTH  '
292          diagTitle = 'Zonal Transp of Pot Temp (K.m/s)'
293          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
294          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
295          CALL DIAGNOSTICS_ADD2LIST( diagNum,
296         I          diagName, diagCode, diagUnits, diagTitle, myThid )
297    
298          diagName  = 'VVELTH  '
299          diagTitle = 'Meridional Transp of Pot Temp (K.m/s)'
300          diagUnits = DIAGS_MK_UNITS( tUnit4c//'.m/s', myThid )
301          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
302          CALL DIAGNOSTICS_ADD2LIST( diagNum,
303         I          diagName, diagCode, diagUnits, diagTitle, myThid )
304    
305          diagName  = 'WVELTH  '
306          diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
307          diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
308          diagCode  = 'WM      LR      '
309          CALL DIAGNOSTICS_ADD2LIST( diagNum,
310         I          diagName, diagCode, diagUnits, diagTitle, myThid )
311    
312          diagName  = 'UVELSLT '
313          diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'
314          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
315          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
316          CALL DIAGNOSTICS_ADD2LIST( diagNum,
317         I          diagName, diagCode, diagUnits, diagTitle, myThid )
318    
319          diagName  = 'VVELSLT '
320          diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'
321          diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
322          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
323          CALL DIAGNOSTICS_ADD2LIST( diagNum,
324         I          diagName, diagCode, diagUnits, diagTitle, myThid )
325    
326          diagName  = 'WVELSLT '
327          diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
328          diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
329          diagCode  = 'WM      LR      '
330          CALL DIAGNOSTICS_ADD2LIST( diagNum,
331         I          diagName, diagCode, diagUnits, diagTitle, myThid )
332    
333    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
334    
335          diagName  = 'RHOAnoma'
336          diagTitle = 'Density Anomaly (=Rho-rhoConst)'
337          diagUnits = 'kg/m^3          '
338          diagCode  = 'SM      MR      '
339          CALL DIAGNOSTICS_ADD2LIST( diagNum,
340         I          diagName, diagCode, diagUnits, diagTitle, myThid )
341    
342          diagName  = 'RHOANOSQ'
343          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
344          diagUnits = 'kg^2/m^6          '
345          diagCode  = 'SM      MR      '
346          CALL DIAGNOSTICS_ADD2LIST( diagNum,
347         I          diagName, diagCode, diagUnits, diagTitle, myThid )
348    
349          diagName  = 'URHOMASS'
350          diagTitle = 'Zonal Transport of Density'
351          diagUnits = 'kg/m^2/s        '
352          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
353          CALL DIAGNOSTICS_ADD2LIST( diagNum,
354         I          diagName, diagCode, diagUnits, diagTitle, myThid )
355    
356          diagName  = 'VRHOMASS'
357          diagTitle = 'Meridional Transport of Density'
358          diagUnits = 'kg/m^2/s        '
359          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
360          CALL DIAGNOSTICS_ADD2LIST( diagNum,
361         I          diagName, diagCode, diagUnits, diagTitle, myThid )
362    
363          diagName  = 'WRHOMASS'
364          diagTitle = 'Vertical Transport of Potential Density'
365          diagUnits = 'kg/m^2/s        '
366          diagCode  = 'WM      LR      '
367          CALL DIAGNOSTICS_ADD2LIST( diagNum,
368         I          diagName, diagCode, diagUnits, diagTitle, myThid )
369    
370        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
371        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
372        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
373        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
374        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
375       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
376    
377        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
378        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
379        diagUnits = 'm2/s2           '        diagUnits = 'm^2/s^2         '
380        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
381        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
382       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
383    
384        diagName  = 'PHIBOTSQ'        diagName  = 'PHIBOTSQ'
385        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'        diagTitle = 'Square of ocean bottom pressure / top. geo-Potential'
386        diagUnits = 'm4/s4           '        diagUnits = 'm^4/s^4         '
387        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
388        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
389       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
390    
391        diagName  = 'DRHODR  '        diagName  = 'DRHODR  '
392        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)       '        diagTitle = 'Stratification: d.Sigma/dr (kg/m3/r_unit)'
393        diagUnits = 'kg/m4           '        diagUnits = 'kg/m^4          '
394          IF ( usingPCoords ) diagUnits = 's^2/m^2         '
395        diagCode  = 'SM      LR      '        diagCode  = 'SM      LR      '
396        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
397       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
398    
399        diagName  = 'DETADT2 '  C--   Momentum terms:
400        diagTitle = 'Square of Eta (Surf.P,SSH) Tendency ((r_unit/s)^2)'        diagName  = 'VISCA4  '
401        diagUnits = 'm2 or Pa2 /sec2 '        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '
402        diagCode  = 'SM      M1      '        diagUnits = 'm^4/s           '
403          diagCode  = 'SM      MR      '
404          CALL DIAGNOSTICS_ADD2LIST( diagNum,
405         I          diagName, diagCode, diagUnits, diagTitle, myThid )
406    
407          diagName  = 'VISCAH  '
408          diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '
409          diagUnits = 'm^2/s           '
410          diagCode  = 'SM      MR      '
411        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
412       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
413    
414    C--   surface fluxes:
415        diagName  = 'TAUX    '        diagName  = 'TAUX    '
416        diagTitle = 'zonal surface wind stress, >0 increases uVel    '        diagTitle = 'zonal surface wind stress, >0 increases uVel    '
417        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
418        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
419        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
420       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
421    
422        diagName  = 'TAUY    '        diagName  = 'TAUY    '
423        diagTitle = 'meridional surf. wind stress, >0 increases vVel '        diagTitle = 'meridional surf. wind stress, >0 increases vVel '
424        diagUnits = 'N/m2            '        diagUnits = 'N/m^2           '
425        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
426        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
427       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
428    
429        diagName  = 'TFLUX   '        diagName  = 'TFLUX   '
430        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'net surface heat flux, >0 increases theta       '
431        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
432        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
433        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
434       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
435    
436        diagName  = 'TRELAX  '        diagName  = 'TRELAX  '
437        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'surface temperature relaxation, >0 increases theta'
438        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
439        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
440        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
441       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
442    
443        diagName  = 'TICE    '        diagName  = 'TICE    '
444        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'
445        diagUnits = 'W/m2            '        diagUnits = 'W/m^2           '
446        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
447        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
448       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
449    
450        diagName  = 'SFLUX   '        diagName  = 'SFLUX   '
451        diagTitle = 'net surface salt flux, >0 increases salt        '        diagTitle = 'net surface salt flux, >0 increases salt        '
452        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
453        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
454        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
455       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
456    
457        diagName  = 'SRELAX  '        diagName  = 'SRELAX  '
458        diagTitle = 'surface salinity relaxation, >0 increases salt  '        diagTitle = 'surface salinity relaxation, >0 increases salt  '
459        diagUnits = 'g/m2/s          '        diagUnits = 'g/m^2/s         '
460        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
461        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
462       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22