/[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.41 by jmc, Mon Nov 7 15:53:31 2011 UTC revision 1.49 by jmc, Tue Feb 4 19:39:57 2014 UTC
# Line 76  C     sTitle     :: "salt" title Line 76  C     sTitle     :: "salt" title
76        CHARACTER*16   diagUnits        CHARACTER*16   diagUnits
77        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
78        CHARACTER*2    rUnit2c        CHARACTER*2    rUnit2c
79        CHARACTER*4    tUnit4c, sUnit4c        CHARACTER*4    tUnit4c
80          CHARACTER*5    sUnit5c
81        CHARACTER*(10) rTitle, eTitle, fTitle        CHARACTER*(10) rTitle, eTitle, fTitle
82        CHARACTER*(20) pTitle, sTitle        CHARACTER*(20) pTitle, sTitle
83    
# Line 102  C--------------------------------------- Line 103  C---------------------------------------
103        ENDIF        ENDIF
104        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
105          tUnit4c= 'K   '          tUnit4c= 'K   '
106          sUnit4c= 'g/kg'          sUnit5c= 'kg/kg'
107          sTitle = ' Specific Humidity  '          sTitle = ' Specific Humidity  '
108  #ifdef ALLOW_FIZHI          IF (useAIM) sUnit5c= 'g/kg '
109          IF (useFIZHI) sUnit4c= 'kg/kg'        ELSEIF ( eosType.EQ.'TEOS10' ) THEN
110  #endif /* ALLOW_FIZHI */          tUnit4c= 'degC'
111            sUnit5c= 'g/kg '
112    c       tTitle = 'Conservative Temp.   '
113            sTitle = ' Absolute Salinity  '
114        ELSE        ELSE
115          tUnit4c= 'degC'          tUnit4c= 'degC'
116          sUnit4c= 'psu '          sUnit5c= 'psu  '
117    c       tTitle = 'Potential Temperature'
118          sTitle = ' Salinity           '          sTitle = ' Salinity           '
119        ENDIF        ENDIF
120  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
# Line 169  c    I          diagName, diagCode, diag Line 174  c    I          diagName, diagCode, diag
174    
175        diagName  = 'SALT    '        diagName  = 'SALT    '
176        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
177        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
178        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
179        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
180       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 183  c    I          diagName, diagCode, diag Line 188  c    I          diagName, diagCode, diag
188    
189  c     diagName  = 'SSS     '  c     diagName  = 'SSS     '
190  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
191  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
192  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
193  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
194  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 191  c    I          diagName, diagCode, diag Line 196  c    I          diagName, diagCode, diag
196        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
197        diagName  = 'SALTanom'        diagName  = 'SALTanom'
198        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
199        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
200        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
201        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
202       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 229  c    I          diagName, diagCode, diag Line 234  c    I          diagName, diagCode, diag
234    
235        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
236        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
237        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
238        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
239        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
240       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 237  c    I          diagName, diagCode, diag Line 242  c    I          diagName, diagCode, diag
242        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
243        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
244        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
245        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
246        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
247        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
248       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 378  C-    use 'PhiVEL' as mate. Line 383  C-    use 'PhiVEL' as mate.
383        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
384        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
385       I                           //sTitle, myThid )       I                           //sTitle, myThid )
386        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
387        diagCode  = 'UUr     MR      '        diagCode  = 'UUr     MR      '
388        diagMate  = diagNum + 2        diagMate  = diagNum + 2
389        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 387  C-    use 'PhiVEL' as mate. Line 392  C-    use 'PhiVEL' as mate.
392        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
393        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
394       I                           //sTitle, myThid )       I                           //sTitle, myThid )
395        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
396        diagCode  = 'VVr     MR      '        diagCode  = 'VVr     MR      '
397        diagMate  = diagNum        diagMate  = diagNum
398        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 396  C-    use 'PhiVEL' as mate. Line 401  C-    use 'PhiVEL' as mate.
401        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
402        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
403       I                           //sTitle, myThid )       I                           //sTitle, myThid )
404        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
405        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
406        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
407       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 427  C-    use 'PhiVEL' as mate. Line 432  C-    use 'PhiVEL' as mate.
432        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
433        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
434       I                          //sTitle, myThid )       I                          //sTitle, myThid )
435        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
436        diagCode  = 'UUR     MR      '        diagCode  = 'UUR     MR      '
437        diagMate  = diagNum + 2        diagMate  = diagNum + 2
438        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 436  C-    use 'PhiVEL' as mate. Line 441  C-    use 'PhiVEL' as mate.
441        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
442        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
443       I                          //sTitle, myThid )       I                          //sTitle, myThid )
444        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
445        diagCode  = 'VVR     MR      '        diagCode  = 'VVR     MR      '
446        diagMate  = diagNum        diagMate  = diagNum
447        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 445  C-    use 'PhiVEL' as mate. Line 450  C-    use 'PhiVEL' as mate.
450        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
451        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
452       I                          //sTitle, myThid )       I                          //sTitle, myThid )
453        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
454        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
455        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
456       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 479  C---+----1----+----2----+----3----+----4 Line 484  C---+----1----+----2----+----3----+----4
484    
485        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
486        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
487        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6        '
488        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
489        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
490       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 501  C---+----1----+----2----+----3----+----4 Line 506  C---+----1----+----2----+----3----+----4
506       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
507    
508        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
509        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Density'
510          diagUnits = 'kg/m^2/s        '
511          diagCode  = 'WM      LR      '
512          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
513         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
514    
515          diagName  = 'WdRHO_P '
516          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
517          diagUnits = 'kg/m^2/s        '
518          diagCode  = 'WM      LR      '
519          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
521    
522          diagName  = 'WdRHOdP '
523          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
524        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
525        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
526        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 741  c     diagTitle = 'Free-Surface r-Positi Line 760  c     diagTitle = 'Free-Surface r-Positi
760    
761        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
762        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
763        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
764        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
765        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
766       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
767    
768          diagName  = 'MoistCor'
769          diagTitle = 'Heating correction due to moist thermodynamics'
770          diagUnits = 'W/m^2           '
771          diagCode  = 'SM      MR      '
772          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
773         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
774    
775        diagName  = 'AB_gU   '        diagName  = 'AB_gU   '
776        diagTitle = 'U momentum tendency from Adams-Bashforth'        diagTitle = 'U momentum tendency from Adams-Bashforth'
777        diagUnits = 'm/s^2           '        diagUnits = 'm/s^2           '
# Line 782  c     diagTitle = 'Free-Surface r-Positi Line 808  c     diagTitle = 'Free-Surface r-Positi
808        diagName  = 'AB_gS   '        diagName  = 'AB_gS   '
809        diagTitle = DIAGS_MK_TITLE(        diagTitle = DIAGS_MK_TITLE(
810       &            sTitle//'tendency from Adams-Bashforth', myThid )       &            sTitle//'tendency from Adams-Bashforth', myThid )
811        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
812          diagCode  = 'SMR     MR      '
813          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
814         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
815    
816          diagName  = 'gTinAB  '
817          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
818          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
819          diagCode  = 'SMR     MR      '
820          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
821         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
822    
823          diagName  = 'gSinAB  '
824          diagTitle = DIAGS_MK_TITLE(
825         &            sTitle//'tendency going in Adams-Bashforth', myThid )
826          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
827        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
828        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
829       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
830  #endif /* ALLOW_GENERIC_ADVDIFF */  #endif /* ALLOW_GENERIC_ADVDIFF */
831    
832    #ifdef ALLOW_EDDYPSI
833          diagName  = 'TAUXEDDY'
834          diagTitle = 'Zonal Eddy Stress'
835          diagUnits = 'N/m**2          '
836          diagCode  = 'UU P    MR      '
837          diagMate  = diagNum
838          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
839         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
840            
841          diagName  = 'TAUYEDDY'
842          diagTitle = 'Meridional Eddy Stress'
843          diagUnits = 'N/m**2          '
844          diagCode  = 'VV P    MR      '
845          diagMate  = diagNum
846          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
847         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
848    
849    #ifdef ALLOW_GMREDI
850          diagName  = 'UVELMEAN'
851          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
852          diagUnits = 'm/s             '
853          diagCode  = 'UUR     MR      '
854          diagMate  = diagNum + 2
855          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
856         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
857            
858          diagName  = 'VVELMEAN'
859          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
860          diagUnits = 'm/s             '
861          diagCode  = 'VVR     MR      '
862          diagMate  = diagNum
863          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
864         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
865    
866    #endif /* ALLOW_GMREDI */
867    
868    #endif /* ALLOW_EDDYPSI */
869    
870        RETURN        RETURN
871        END        END
Legend:
Removed from v.1.41  
changed lines
  Added in v.1.49
  ViewVC Help
Powered by ViewVC 1.1.22