/[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.9 by jmc, Thu Feb 10 15:01:34 2005 UTC revision 1.19 by molod, Tue Aug 16 20:36:10 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 84  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 120  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     '
130          diagTitle = 'Sea Surface Temperature (degC,K)'
131          diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid )
132          diagCode  = 'SM      M1      '
133          CALL DIAGNOSTICS_ADD2LIST( diagNum,
134         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 )
139          diagCode  = 'SM      MR      '
140          CALL DIAGNOSTICS_ADD2LIST( diagNum,
141         I          diagName, diagCode, diagUnits, diagTitle, myThid )
142    
143          diagName  = 'SSS     '
144          diagTitle = 'Sea Surface Salinity '
145          diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
146          diagCode  = 'SM      M1      '
147          CALL DIAGNOSTICS_ADD2LIST( diagNum,
148         I          diagName, diagCode, diagUnits, diagTitle, myThid )
149    
150          diagName  = 'SALTanom'
151          diagTitle = 'Salt anomaly (=SALT-35)'
152        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )
153        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
154        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
# Line 141  C-    state variables of the main code ( Line 168  C-    state variables of the main code (
168        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
169       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
170    
171          diagName  = 'UVEL_k2 '
172          diagTitle = 'Zonal Component of Velocity at level 2 (m/s)'
173          diagUnits = 'm/s             '
174          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
175          CALL DIAGNOSTICS_ADD2LIST( diagNum,
176         I          diagName, diagCode, diagUnits, diagTitle, myThid )
177    
178          diagName  = 'VVEL_k2 '
179          diagTitle = 'Meridional Component of Velocity at level 2 (m/s)'
180          diagUnits = 'm/s             '
181          WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'M1      '
182          CALL DIAGNOSTICS_ADD2LIST( diagNum,
183         I          diagName, diagCode, diagUnits, diagTitle, myThid )
184    
185        diagName  = 'WVEL    '        diagName  = 'WVEL    '
186        diagTitle = 'Vertical Component of Velocity (r_units/s)'        diagTitle = 'Vertical Component of Velocity (r_units/s)'
187        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
# Line 162  C-    state variables of the main code ( Line 203  C-    state variables of the main code (
203        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
204       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
205    
206          diagName  = 'SALTSQan'
207          diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
208          diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )
209          diagCode  = 'SM      MR      '
210          CALL DIAGNOSTICS_ADD2LIST( diagNum,
211         I          diagName, diagCode, diagUnits, diagTitle, myThid )
212    
213        diagName  = 'UVELSQ  '        diagName  = 'UVELSQ  '
214        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'        diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)'
215        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 228  C-    state variables of the main code ( Line 276  C-    state variables of the main code (
276        diagName  = 'WVELMASS'        diagName  = 'WVELMASS'
277        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'        diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
278        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid )
279        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
280        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
281       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
282    
# Line 249  C-    state variables of the main code ( Line 297  C-    state variables of the main code (
297        diagName  = 'WTHMASS '        diagName  = 'WTHMASS '
298        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)'
299        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
300        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
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  = 'USLTMASS'        diagName  = 'USLTMASS'
305        diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal 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)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'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  = 'VSLTMASS'        diagName  = 'VSLTMASS'
312        diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional 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)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'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  = 'WSLTMASS'        diagName  = 'WSLTMASS'
319        diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Mass-Weight Transp of Salt '
320        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
321        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
322        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
323       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
324    
# Line 291  C-    state variables of the main code ( Line 339  C-    state variables of the main code (
339        diagName  = 'WVELTH  '        diagName  = 'WVELTH  '
340        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'        diagTitle = 'Vertical Transp of Pot Temp (K.m/s)'
341        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid )
342        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
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  = 'UVELSLT '        diagName  = 'UVELSLT '
347        diagTitle = 'Zonal Transp of Salt (g/kg.m/s)'        diagTitle = 'Zonal 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)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'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  = 'VVELSLT '        diagName  = 'VVELSLT '
354        diagTitle = 'Meridional Transp of Salt (g/kg.m/s)'        diagTitle = 'Meridional 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)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'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  = 'WVELSLT '        diagName  = 'WVELSLT '
361        diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'        diagTitle = 'Vertical Transp of Salt '
362        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )
363        diagCode  = 'WM      MR      '        diagCode  = 'WM      LR      '
364        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
365       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
366    
# Line 325  C---+----1----+----2----+----3----+----4 Line 373  C---+----1----+----2----+----3----+----4
373        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
374       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
375    
376          diagName  = 'RHOANOSQ'
377          diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
378          diagUnits = 'kg^2/m^6          '
379          diagCode  = 'SM      MR      '
380          CALL DIAGNOSTICS_ADD2LIST( diagNum,
381         I          diagName, diagCode, diagUnits, diagTitle, myThid )
382    
383          diagName  = 'URHOMASS'
384          diagTitle = 'Zonal Transport of Density'
385          diagUnits = 'kg/m^2/s        '
386          WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
387          CALL DIAGNOSTICS_ADD2LIST( diagNum,
388         I          diagName, diagCode, diagUnits, diagTitle, myThid )
389    
390          diagName  = 'VRHOMASS'
391          diagTitle = 'Meridional Transport of Density'
392          diagUnits = 'kg/m^2/s        '
393          WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
394          CALL DIAGNOSTICS_ADD2LIST( diagNum,
395         I          diagName, diagCode, diagUnits, diagTitle, myThid )
396    
397          diagName  = 'WRHOMASS'
398          diagTitle = 'Vertical Transport of Potential Density'
399          diagUnits = 'kg/m^2/s        '
400          diagCode  = 'WM      LR      '
401          CALL DIAGNOSTICS_ADD2LIST( diagNum,
402         I          diagName, diagCode, diagUnits, diagTitle, myThid )
403    
404        diagName  = 'PHIHYD  '        diagName  = 'PHIHYD  '
405        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'        diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential'
406        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 332  C---+----1----+----2----+----3----+----4 Line 408  C---+----1----+----2----+----3----+----4
408        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
409       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
410    
411          diagName  = 'PHIHYDSQ'
412          diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential'
413          diagUnits = 'm^4/s^4         '
414          diagCode  = 'SM      MR      '
415          CALL DIAGNOSTICS_ADD2LIST( diagNum,
416         I          diagName, diagCode, diagUnits, diagTitle, myThid )
417    
418        diagName  = 'PHIBOT  '        diagName  = 'PHIBOT  '
419        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'        diagTitle = 'ocean bottom pressure / top. atmos geo-Potential'
420        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
# Line 419  C--   surface fluxes: Line 502  C--   surface fluxes:
502        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
503       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
504    
505    #ifdef SHORTWAVE_HEATING
506          diagName  = 'SWFLUX  '
507          diagTitle = 'net upward SW radiation, >0 increases theta     '
508          diagUnits = 'W/m^2           '
509          diagCode  = 'SM      U1      '
510          CALL DIAGNOSTICS_ADD2LIST( diagNum,
511         I          diagName, diagCode, diagUnits, diagTitle, myThid )
512    #endif
513    
514  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
515          diagName  = 'PRESSURE'
516          diagTitle = 'Atmospheric Pressure (Pa)     '
517          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
518          diagCode  = 'SM      MR      '
519          CALL DIAGNOSTICS_ADD2LIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, myThid )
521    
522          diagName  = 'RSURF   '
523          diagTitle = 'Surface (Pressure, Height) (Pa,m)'
524          diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )
525          diagCode  = 'SM      M1      '
526          CALL DIAGNOSTICS_ADD2LIST( diagNum,
527         I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529    
530        RETURN        RETURN
531        END        END
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22