/[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.10 by jmc, Tue Feb 15 20:57:38 2005 UTC revision 1.18 by heimbach, Fri Jul 29 15:50:53 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 360  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 447  C--   surface fluxes: Line 497  C--   surface fluxes:
497        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
498       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500    #ifdef SHORTWAVE_HEATING
501          diagName  = 'SWFLUX  '
502          diagTitle = 'net upward SW radiation, >0 increases theta     '
503          diagUnits = 'W/m^2           '
504          diagCode  = 'SM      U1      '
505          CALL DIAGNOSTICS_ADD2LIST( diagNum,
506         I          diagName, diagCode, diagUnits, diagTitle, myThid )
507    #endif
508    
509  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
510          diagName  = 'PRESSURE'
511          diagTitle = 'Atmospheric Pressure (Pa)     '
512          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
513          diagCode  = 'SM      MR      '
514          CALL DIAGNOSTICS_ADD2LIST( diagNum,
515         I          diagName, diagCode, diagUnits, diagTitle, myThid )
516    
517          diagName  = 'RSURF   '
518          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
519          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
520          diagCode  = 'SM      M1      '
521          CALL DIAGNOSTICS_ADD2LIST( diagNum,
522         I          diagName, diagCode, diagUnits, diagTitle, myThid )
523    
524    
525        RETURN        RETURN
526        END        END
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.22