/[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.16 by molod, Wed Jun 15 12:59:10 2005 UTC revision 1.20 by molod, Tue Aug 16 21:29:38 2005 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 129  C-    state variables of the main code ( Line 134  C-    state variables of the main code (
134       I          diagName, diagCode, diagUnits, diagTitle, myThid )       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  = 'RELHUM  '
144          diagTitle = 'Relative Humidity               '
145          diagUnits = 'percent         '
146          diagCode  = 'SM      MR      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150        diagName  = 'SSS     '        diagName  = 'SSS     '
151        diagTitle = 'Sea Surface Salinity (g/kg)'        diagTitle = 'Sea Surface Salinity '
152        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153        diagCode  = 'SM      M1      '        diagCode  = 'SM      M1      '
154        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
155       I          diagName, diagCode, diagUnits, diagTitle, myThid )       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)'
159        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
160        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
161        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 297  C-    state variables of the main code ( Line 309  C-    state variables of the main code (
309       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
310    
311        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
312        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal Mass-Weight Transp of Salt '
313        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
314        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
315        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
316       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
317    
318        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
319        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional Mass-Weight Transp of Salt'
320        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )
321        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
322        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
323       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
325        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
326        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Salt '
327        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
328        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
329        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 339  C-    state variables of the main code ( Line 351  C-    state variables of the main code (
351       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
352    
353        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
354        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal Transp of Salt '
355        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
356        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
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  = 'VVELSLT '        diagName  = 'VVELSLT '
361        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional Transp of Salt '
362        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )
363        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
364        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
365       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
367        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
368        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Transp of Salt '
369        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
370        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
371        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 497  C--   surface fluxes: Line 509  C--   surface fluxes:
509        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
510       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
511    
512    #ifdef SHORTWAVE_HEATING
513          diagName  = 'SWFLUX  '
514          diagTitle = 'net upward SW radiation, >0 increases theta     '
515          diagUnits = 'W/m^2           '
516          diagCode  = 'SM      U1      '
517          CALL DIAGNOSTICS_ADD2LIST( diagNum,
518         I          diagName, diagCode, diagUnits, diagTitle, myThid )
519    #endif
520    
521  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
522        diagName  = 'PRESSURE'        diagName  = 'PRESSURE'
523        diagTitle = 'Atmospheric Pressure (Pa)     '        diagTitle = 'Atmospheric Pressure (Pa)     '
# Line 505  C---+----1----+----2----+----3----+----4 Line 526  C---+----1----+----2----+----3----+----4
526        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
527       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529          diagName  = 'RSURF   '
530          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
531          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
532          diagCode  = 'SM      M1      '
533          CALL DIAGNOSTICS_ADD2LIST( diagNum,
534         I          diagName, diagCode, diagUnits, diagTitle, myThid )
535    
536    
537        RETURN        RETURN
538        END        END
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22