/[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.18 by heimbach, Fri Jul 29 15:50:53 2005 UTC revision 1.25 by molod, Mon Feb 27 17:55:14 2006 UTC
# Line 85  C--------------------------------------- Line 85  C---------------------------------------
85        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
86          tUnit4c= 'K   '          tUnit4c= 'K   '
87          sUnit4c= 'g/kg'          sUnit4c= 'g/kg'
88    #ifdef ALLOW_FIZHI
89          IF(useFIZHI)THEN
90            sUnit4c= 'kg/kg'
91          ENDIF
92    #endif /* ALLOW_FIZHI */
93        ELSE        ELSE
94          tUnit4c= 'degC'          tUnit4c= 'degC'
95          sUnit4c= 'psu '          sUnit4c= 'psu '
# Line 101  C-    state variables of the main code ( Line 106  C-    state variables of the main code (
106       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
107    
108        diagName  = 'ETANSQ  '        diagName  = 'ETANSQ  '
109        diagTitle = 'Square of Perturbation of Sfc (Pa^2,m^2)'        diagTitle = 'Square of Perturbation of Sfc position (r_unit^2)'
110        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'^2', myThid )
111        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
112        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 121  C-    state variables of the main code ( Line 126  C-    state variables of the main code (
126        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
127       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
128    
129        diagName  = 'SST     '  c     diagName  = 'SST     '
130        diagTitle = 'Sea Surface Temperature (degC,K)'  c     diagTitle = 'Sea Surface Temperature (degC,K)'
131        diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
132        diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
133        CALL DIAGNOSTICS_ADD2LIST( diagNum,  c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
134       I          diagName, diagCode, diagUnits, diagTitle, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
135    
136        diagName  = 'SALT    '        diagName  = 'SALT    '
137        diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)'        diagTitle = 'Salt or Water Vapor mixing ratio'
138        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
139        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
140        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
141       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
142    
143        diagName  = 'SSS     '        diagName  = 'RELHUM  '
144        diagTitle = 'Sea Surface Salinity (g/kg)'        diagTitle = 'Relative Humidity               '
145        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = 'percent         '
146        diagCode  = 'SM      M1      '        diagCode  = 'SM      MR      '
147        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
148       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150    c     diagName  = 'SSS     '
151    c     diagTitle = 'Sea Surface Salinity '
152    c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153    c     diagCode  = 'SM      M1      '
154    c     CALL DIAGNOSTICS_ADD2LIST( diagNum,
155    c    I          diagName, diagCode, diagUnits, diagTitle, myThid )
156    
157        diagName  = 'SALTanom'        diagName  = 'SALTanom'
158        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
159        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
# Line 163  C-    state variables of the main code ( Line 175  C-    state variables of the main code (
175        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
176       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
       diagName  = 'UVEL_k2 '  
       diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'  
       diagUnits = 'm/s             '  
       WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
   
       diagName  = 'VVEL_k2 '  
       diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'  
       diagUnits = 'm/s             '  
       WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
   
178        diagName  = 'WVEL    '        diagName  = 'WVEL    '
179        diagTitle = 'Vertical Component of Velocity (r_units/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
180        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
# Line 297  C-    state variables of the main code ( Line 295  C-    state variables of the main code (
295       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
296    
297        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
298        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal Mass-Weight Transp of Salt '
299        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
300        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
301        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
302       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
303    
304        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
305        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional Mass-Weight Transp of Salt'
306        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
307        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
308        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
309       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
310    
311        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
312        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Salt '
313        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
314        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 339  C-    state variables of the main code ( Line 337  C-    state variables of the main code (
337       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
338    
339        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
340        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal Transp of Salt '
341        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
342        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
343        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
344       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
345    
346        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
347        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional Transp of Salt '
348        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
349        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
350        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
351       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
354        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Transp of Salt '
355        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
356        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
357        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
358       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
359    
360          diagName  = 'UVELPHI '
361          diagTitle = 'Zonal Transp of Geopotential Height (m2/s2.m/s)'
362          diagUnits = 'm^3/s^3         '
363          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
364          CALL DIAGNOSTICS_ADD2LIST( diagNum,
365         I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367          diagName  = 'VVELPHI '
368          diagTitle = 'Meridional Transp of Geopotential Height (m2/s2.m/s)'
369          diagUnits = 'm^3/s^3         '
370          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
371          CALL DIAGNOSTICS_ADD2LIST( diagNum,
372         I          diagName, diagCode, diagUnits, diagTitle, myThid )
373    
374  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
375    
376        diagName  = 'RHOAnoma'        diagName  = 'RHOAnoma'
# Line 432  C---+----1----+----2----+----3----+----4 Line 444  C---+----1----+----2----+----3----+----4
444        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
445       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
446    
447  C--   Momentum terms:        diagName  = 'CONVADJ '
448        diagName  = 'VISCA4  '        diagTitle = 'Convective Adjustment Index [0-1] '
449        diagTitle = 'Biharmonic Viscosity Coefficient in (m4/s)        '        diagUnits = 'fraction        '
450        diagUnits = 'm^4/s           '        diagCode  = 'SM      LR      '
       diagCode  = 'SM      MR      '  
       CALL DIAGNOSTICS_ADD2LIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, myThid )  
   
       diagName  = 'VISCAH  '  
       diagTitle = 'Harmonic Viscosity Coefficient in (m2/s)        '  
       diagUnits = 'm^2/s           '  
       diagCode  = 'SM      MR      '  
451        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
452       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
453    
# Line 521  C---+----1----+----2----+----3----+----4 Line 525  C---+----1----+----2----+----3----+----4
525        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
526       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
527    
528          diagName  = 'TOTUTEND'
529          diagTitle = 'Tendency of Zonal Component of Velocity (m/s^2)'
530          diagUnits = 'm/s^2           '
531          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
532          CALL DIAGNOSTICS_ADD2LIST( diagNum,
533         I          diagName, diagCode, diagUnits, diagTitle, myThid )
534    
535          diagName  = 'TOTVTEND'
536          diagTitle = 'Tendency of Meridional Component of Velocity (m/s^2)'
537          diagUnits = 'm/s^2           '
538          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
539          CALL DIAGNOSTICS_ADD2LIST( diagNum,
540         I          diagName, diagCode, diagUnits, diagTitle, myThid )
541    
542          diagName  = 'TOTTTEND'
543          diagTitle = 'Tendency of Potential Temperature (degC,K)'
544          diagUnits = 'deg K/sec'
545          diagCode  = 'SM      MR      '
546          CALL DIAGNOSTICS_ADD2LIST( diagNum,
547         I          diagName, diagCode, diagUnits, diagTitle, myThid )
548    
549          diagName  = 'TOTSTEND'
550          diagTitle = 'Tendency of Salt or Water Vapor mixing ratio'
551          diagUnits = 'kg/kg-sec'
552          diagCode  = 'SM      MR      '
553          CALL DIAGNOSTICS_ADD2LIST( diagNum,
554         I          diagName, diagCode, diagUnits, diagTitle, myThid )
555    
556    
557        RETURN        RETURN
558        END        END
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22