/[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.8 by jmc, Fri Jan 28 01:06:12 2005 UTC revision 1.17 by molod, Mon Jul 11 18:34:17 2005 UTC
# Line 44  C                   &  1  &  levels = 1 Line 44  C                   &  1  &  levels = 1
44  C                   &  R  &  levels = Nr  \\  C                   &  R  &  levels = Nr  \\
45  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\  C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
46  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\  C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
47    C                   &  G  &  levels = Ground_level Number \\
48    C                   &  I  &  levels = sea-Ice_level Number \\
49  C       \end{tabular}  C       \end{tabular}
50  C     \end{center}  C     \end{center}
51    
# Line 51  C     !USES: Line 53  C     !USES:
53        IMPLICIT NONE        IMPLICIT NONE
54  #include "SIZE.h"  #include "SIZE.h"
55  #include "EEPARAMS.h"  #include "EEPARAMS.h"
 c #include "EESUPPORT.h"  
56  #include "PARAMS.h"  #include "PARAMS.h"
57    
58  C     !INPUT PARAMETERS:  C     !INPUT PARAMETERS:
# Line 120  C-    state variables of the main code ( Line 121  C-    state variables of the main code (
121        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
122       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
123    
124          diagName  = 'SST     '
125          diagTitle = 'Sea Surface Temperature (degC,K)'
126          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
127          diagCode  = 'SM      M1      '
128          CALL DIAGNOSTICS_ADD2LIST( diagNum,
129         I          diagName, diagCode, diagUnits, diagTitle, myThid )
130    
131        diagName  = 'SALT    '        diagName  = 'SALT    '
132        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'
133        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
# Line 127  C-    state variables of the main code ( Line 135  C-    state variables of the main code (
135        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
136       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
137    
138          diagName  = 'SSS     '
139          diagTitle = 'Sea Surface Salinity (g/kg)'
140          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
141          diagCode  = 'SM      M1      '
142          CALL DIAGNOSTICS_ADD2LIST( diagNum,
143         I          diagName, diagCode, diagUnits, diagTitle, myThid )
144    
145          diagName  = 'SALTanom'
146          diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
147          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
148          diagCode  = 'SM      MR      '
149          CALL DIAGNOSTICS_ADD2LIST( diagNum,
150         I          diagName, diagCode, diagUnits, diagTitle, myThid )
151    
152        diagName  = 'UVEL    '        diagName  = 'UVEL    '
153        diagTitle = 'Zonal Component of Velocity (m/s)'        diagTitle = 'Zonal Component of Velocity (m/s)'
154        diagUnits = 'm/s             '        diagUnits = 'm/s             '
# Line 141  C-    state variables of the main code ( Line 163  C-    state variables of the main code (
163        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
164       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
165    
166          diagName  = 'UVEL_k2 '
167          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
168          diagUnits = 'm/s             '
169          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
170          CALL DIAGNOSTICS_ADD2LIST( diagNum,
171         I          diagName, diagCode, diagUnits, diagTitle, myThid )
172    
173          diagName  = 'VVEL_k2 '
174          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
175          diagUnits = 'm/s             '
176          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
177          CALL DIAGNOSTICS_ADD2LIST( diagNum,
178         I          diagName, diagCode, diagUnits, diagTitle, myThid )
179    
180        diagName  = 'WVEL    '        diagName  = 'WVEL    '
181        diagTitle = 'Vertical Component of Velocity (r_units/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
182        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
# Line 162  C-    state variables of the main code ( Line 198  C-    state variables of the main code (
198        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
199       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
200    
201          diagName  = 'SALTSQan'
202          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
203          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
204          diagCode  = 'SM      MR      '
205          CALL DIAGNOSTICS_ADD2LIST( diagNum,
206         I          diagName, diagCode, diagUnits, diagTitle, myThid )
207    
208        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
209        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
210        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 197  C-    state variables of the main code ( Line 240  C-    state variables of the main code (
240        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
241       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
242    
243          diagName  = 'WU_VEL  '
244          diagTitle = 'Vertical Transport of Zonal Momentum (m^2/s^2)'
245          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
246          diagCode  = 'WU      LR      '
247          CALL DIAGNOSTICS_ADD2LIST( diagNum,
248         I          diagName, diagCode, diagUnits, diagTitle, myThid )
249    
250          diagName  = 'WV_VEL  '
251          diagTitle ='Vertical Transport of Meridional Momentum (m^2/s^2)'
252          diagUnits = DIAGS_MK_UNITS( 'm.'//rUnit2c//'/s^2', myThid )
253          diagCode  = 'WV      LR      '
254          CALL DIAGNOSTICS_ADD2LIST( diagNum,
255         I          diagName, diagCode, diagUnits, diagTitle, myThid )
256    
257        diagName  = 'UVELMASS'        diagName  = 'UVELMASS'
258        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
# Line 215  C-    state variables of the main code ( Line 271  C-    state variables of the main code (
271        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
272        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
273        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
274        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
275        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
276       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
277    
# Line 236  C-    state variables of the main code ( Line 292  C-    state variables of the main code (
292        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
293        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
294        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
295        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
296        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
297       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
298    
# Line 257  C-    state variables of the main code ( Line 313  C-    state variables of the main code (
313        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
314        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
315        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
316        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
317        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
318       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
319    
# Line 278  C-    state variables of the main code ( Line 334  C-    state variables of the main code (
334        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
335        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
336        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
337        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
338        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
339       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
340    
# Line 299  C-    state variables of the main code ( Line 355  C-    state variables of the main code (
355        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
356        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
357        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
358        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
359        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
360       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
361    
# Line 312  C---+----1----+----2----+----3----+----4 Line 368  C---+----1----+----2----+----3----+----4
368        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
369       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
370    
371          diagName  = 'RHOANOSQ'
372          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
373          diagUnits = 'kg^2/m^6          '
374          diagCode  = 'SM      MR      '
375          CALL DIAGNOSTICS_ADD2LIST( diagNum,
376         I          diagName, diagCode, diagUnits, diagTitle, myThid )
377    
378          diagName  = 'URHOMASS'
379          diagTitle = 'Zonal Transport of Density'
380          diagUnits = 'kg/m^2/s        '
381          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
382          CALL DIAGNOSTICS_ADD2LIST( diagNum,
383         I          diagName, diagCode, diagUnits, diagTitle, myThid )
384    
385          diagName  = 'VRHOMASS'
386          diagTitle = 'Meridional Transport of Density'
387          diagUnits = 'kg/m^2/s        '
388          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
389          CALL DIAGNOSTICS_ADD2LIST( diagNum,
390         I          diagName, diagCode, diagUnits, diagTitle, myThid )
391    
392          diagName  = 'WRHOMASS'
393          diagTitle = 'Vertical Transport of Potential Density'
394          diagUnits = 'kg/m^2/s        '
395          diagCode  = 'WM      LR      '
396          CALL DIAGNOSTICS_ADD2LIST( diagNum,
397         I          diagName, diagCode, diagUnits, diagTitle, myThid )
398    
399        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
400        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
401        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 319  C---+----1----+----2----+----3----+----4 Line 403  C---+----1----+----2----+----3----+----4
403        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
404       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406          diagName  = 'PHIHYDSQ'
407          diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'
408          diagUnits = 'm^4/s^4         '
409          diagCode  = 'SM      MR      '
410          CALL DIAGNOSTICS_ADD2LIST( diagNum,
411         I          diagName, diagCode, diagUnits, diagTitle, myThid )
412    
413        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
414        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
415        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 407  C--   surface fluxes: Line 498  C--   surface fluxes:
498       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
501          diagName  = 'PRESSURE'
502          diagTitle = 'Atmospheric Pressure (Pa)     '
503          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
504          diagCode  = 'SM      MR      '
505          CALL DIAGNOSTICS_ADD2LIST( diagNum,
506         I          diagName, diagCode, diagUnits, diagTitle, myThid )
507    
508          diagName  = 'RSURF   '
509          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
510          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
511          diagCode  = 'SM      M1      '
512          CALL DIAGNOSTICS_ADD2LIST( diagNum,
513         I          diagName, diagCode, diagUnits, diagTitle, myThid )
514    
515    
516        RETURN        RETURN
517        END        END
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.22