/[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.46 by m_bates, Wed Sep 18 01:41:05 2013 UTC revision 1.50 by jmc, Sat Jul 19 16:46:25 2014 UTC
# Line 132  C-    free-surface (eTitle) and fixed-bo Line 132  C-    free-surface (eTitle) and fixed-bo
132  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
134    
   
135        diagName  = 'ETAN    '        diagName  = 'ETAN    '
136        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
137  c     IF ( fluidIsWater .AND. usingZCoords )  c     IF ( fluidIsWater .AND. usingZCoords )
# Line 465  C-    use 'PhiVEL' as mate. Line 464  C-    use 'PhiVEL' as mate.
464       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
467        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Merid. Mass-Weight Transp of '
468       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
469        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
470        diagCode  = 'VVr     MR      '        diagCode  = 'VVr     MR      '
# Line 484  C---+----1----+----2----+----3----+----4 Line 483  C---+----1----+----2----+----3----+----4
483    
484        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
485        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
486        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6        '
487        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
488        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
489       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 767  c     diagTitle = 'Free-Surface r-Positi Line 766  c     diagTitle = 'Free-Surface r-Positi
766    
767        diagName  = 'MoistCor'        diagName  = 'MoistCor'
768        diagTitle = 'Heating correction due to moist thermodynamics'        diagTitle = 'Heating correction due to moist thermodynamics'
769        diagUnits = 'W/s^2           '        diagUnits = 'W/m^2           '
770        diagCode  = 'SM      MR      '        diagCode  = 'SM      MR      '
771        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
772       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
773    
774        diagName  = 'AB_gU   '  #ifdef ALLOW_GENERIC_ADVDIFF
775        diagTitle = 'U momentum tendency from Adams-Bashforth'        diagName  = 'gT_Forc '
776        diagUnits = 'm/s^2           '        diagTitle = 'Potential Temp. forcing tendency'
777        diagCode  = 'UUR     MR      '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
778        diagMate  = diagNum + 2        diagCode  = 'SMR     MR      '
       CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
   
       diagName  = 'AB_gV   '  
       diagTitle = 'V momentum tendency from Adams-Bashforth'  
       diagUnits = 'm/s^2           '  
       diagCode  = 'VVR     MR      '  
       diagMate  = diagNum  
779        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
780       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
781    
782  #ifdef ALLOW_NONHYDROSTATIC        diagName  = 'gS_Forc '
783        diagName  = 'AB_gW   '        diagTitle = DIAGS_MK_TITLE(
784        diagTitle = 'W momentum tendency from Adams-Bashforth'       &            sTitle//'forcing tendency', myThid )
785        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
786        diagCode  = 'WM      LR      '        diagCode  = 'SMR     MR      '
787        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
788       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
 #endif /* ALLOW_NONHYDROSTATIC */  
789    
 #ifdef ALLOW_GENERIC_ADVDIFF  
790        diagName  = 'AB_gT   '        diagName  = 'AB_gT   '
791        diagTitle = 'Potential Temp. tendency from Adams-Bashforth'        diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
792        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
# Line 829  c     diagTitle = 'Free-Surface r-Positi Line 818  c     diagTitle = 'Free-Surface r-Positi
818       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
819  #endif /* ALLOW_GENERIC_ADVDIFF */  #endif /* ALLOW_GENERIC_ADVDIFF */
820    
821  #ifdef ALLOW_DIAGNOSTICS        diagName  = 'AB_gU   '
822          diagName  = 'TAUXEDDY'        diagTitle = 'U momentum tendency from Adams-Bashforth'
823          diagTitle = 'Zonal Eddy Stress'        diagUnits = 'm/s^2           '
824          diagUnits = 'N/m**2          '        diagCode  = 'UUR     MR      '
825          diagCode  = 'UU P    MR      '        diagMate  = diagNum + 2
826          diagMate  = diagNum        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
827          CALL DIAGNOSTICS_ADDTOLIST( diagNum,       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
828       I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
829                  diagName  = 'AB_gV   '
830          diagName  = 'TAUYEDDY'        diagTitle = 'V momentum tendency from Adams-Bashforth'
831          diagTitle = 'Meridional Eddy Stress'        diagUnits = 'm/s^2           '
832          diagUnits = 'N/m**2          '        diagCode  = 'VVR     MR      '
833          diagCode  = 'VV P    MR      '        diagMate  = diagNum
834          diagMate  = diagNum        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
835          CALL DIAGNOSTICS_ADDTOLIST( diagNum,       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
836       I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
837  #endif /* ALLOW_DIAGNOSTICS */  #ifdef ALLOW_NONHYDROSTATIC
838          diagName  = 'AB_gW   '
839          diagTitle = 'W momentum tendency from Adams-Bashforth'
840          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
841          diagCode  = 'WM      LR      '
842          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
843         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
844    #endif /* ALLOW_NONHYDROSTATIC */
845    
846    #ifdef ALLOW_EDDYPSI
847          diagName  = 'TAUXEDDY'
848          diagTitle = 'Zonal Eddy Stress'
849          diagUnits = 'N/m**2          '
850          diagCode  = 'UU P    MR      '
851          diagMate  = diagNum
852          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
853         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
854    
855          diagName  = 'TAUYEDDY'
856          diagTitle = 'Meridional Eddy Stress'
857          diagUnits = 'N/m**2          '
858          diagCode  = 'VV P    MR      '
859          diagMate  = diagNum
860          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
861         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
862    
863    # ifdef ALLOW_GMREDI
864          diagName  = 'UVELMEAN'
865          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
866          diagUnits = 'm/s             '
867          diagCode  = 'UUR     MR      '
868          diagMate  = diagNum + 2
869          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872          diagName  = 'VVELMEAN'
873          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
874          diagUnits = 'm/s             '
875          diagCode  = 'VVR     MR      '
876          diagMate  = diagNum
877          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
878         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
879    # endif /* ALLOW_GMREDI */
880    #endif /* ALLOW_EDDYPSI */
881    
882        RETURN        RETURN
883        END        END
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22