/[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.49 by jmc, Tue Feb 4 19:39:57 2014 UTC revision 1.51 by jmc, Tue Jan 20 20:50:15 2015 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 772  c     diagTitle = 'Free-Surface r-Positi Line 771  c     diagTitle = 'Free-Surface r-Positi
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_FRICTION_HEATING
775        diagTitle = 'U momentum tendency from Adams-Bashforth'        diagName  = 'HeatDiss'
776        diagUnits = 'm/s^2           '        diagTitle = 'Heating from frictional dissipation'
777        diagCode  = 'UUR     MR      '        diagUnits = 'W/m^2           '
778        diagMate  = diagNum + 2        diagCode  = 'SM      MR      '
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    #endif /* ALLOW_FRICTION_HEATING */
782    
783        diagName  = 'AB_gV   '  #ifdef ALLOW_GENERIC_ADVDIFF
784        diagTitle = 'V momentum tendency from Adams-Bashforth'        diagName  = 'gT_Forc '
785        diagUnits = 'm/s^2           '        diagTitle = 'Potential Temp. forcing tendency'
786        diagCode  = 'VVR     MR      '        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
787        diagMate  = diagNum        diagCode  = 'SMR     MR      '
788        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
789       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
790    
791  #ifdef ALLOW_NONHYDROSTATIC        diagName  = 'gS_Forc '
792        diagName  = 'AB_gW   '        diagTitle = DIAGS_MK_TITLE(
793        diagTitle = 'W momentum tendency from Adams-Bashforth'       &            sTitle//'forcing tendency', myThid )
794        diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
795        diagCode  = 'WM      LR      '        diagCode  = 'SMR     MR      '
796        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
797       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
 #endif /* ALLOW_NONHYDROSTATIC */  
798    
 #ifdef ALLOW_GENERIC_ADVDIFF  
799        diagName  = 'AB_gT   '        diagName  = 'AB_gT   '
800        diagTitle = 'Potential Temp. tendency from Adams-Bashforth'        diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
801        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
# Line 829  c     diagTitle = 'Free-Surface r-Positi Line 827  c     diagTitle = 'Free-Surface r-Positi
827       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
828  #endif /* ALLOW_GENERIC_ADVDIFF */  #endif /* ALLOW_GENERIC_ADVDIFF */
829    
830          diagName  = 'AB_gU   '
831          diagTitle = 'U momentum tendency from Adams-Bashforth'
832          diagUnits = 'm/s^2           '
833          diagCode  = 'UUR     MR      '
834          diagMate  = diagNum + 2
835          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
836         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
837    
838          diagName  = 'AB_gV   '
839          diagTitle = 'V momentum tendency from Adams-Bashforth'
840          diagUnits = 'm/s^2           '
841          diagCode  = 'VVR     MR      '
842          diagMate  = diagNum
843          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
844         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
845    
846    #ifdef ALLOW_NONHYDROSTATIC
847          diagName  = 'AB_gW   '
848          diagTitle = 'W momentum tendency from Adams-Bashforth'
849          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
850          diagCode  = 'WM      LR      '
851          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
852         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
853    #endif /* ALLOW_NONHYDROSTATIC */
854    
855  #ifdef ALLOW_EDDYPSI  #ifdef ALLOW_EDDYPSI
856        diagName  = 'TAUXEDDY'        diagName  = 'TAUXEDDY'
857        diagTitle = 'Zonal Eddy Stress'        diagTitle = 'Zonal Eddy Stress'
858        diagUnits = 'N/m**2          '        diagUnits = 'N/m^2           '
859        diagCode  = 'UU P    MR      '        diagCode  = 'UU      LR      '
860        diagMate  = diagNum        diagMate  = diagNum + 2
861        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
862       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
863            
864        diagName  = 'TAUYEDDY'        diagName  = 'TAUYEDDY'
865        diagTitle = 'Meridional Eddy Stress'        diagTitle = 'Meridional Eddy Stress'
866        diagUnits = 'N/m**2          '        diagUnits = 'N/m^2           '
867        diagCode  = 'VV P    MR      '        diagCode  = 'VV      LR      '
868        diagMate  = diagNum        diagMate  = diagNum
869        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872  #ifdef ALLOW_GMREDI  # ifdef ALLOW_GMREDI
873        diagName  = 'UVELMEAN'        diagName  = 'U_EulerM'
874        diagTitle = 'Zonal Component of Mean Velocity (m/s)'        diagTitle = 'Zonal Eulerian-Mean Velocity (m/s)'
875        diagUnits = 'm/s             '        diagUnits = 'm/s             '
876        diagCode  = 'UUR     MR      '        diagCode  = 'UUR     MR      '
877        diagMate  = diagNum + 2        diagMate  = diagNum + 2
878        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
879       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
880            
881        diagName  = 'VVELMEAN'        diagName  = 'V_EulerM'
882        diagTitle = 'Meridional Component of Mean Velocity (m/s)'        diagTitle = 'Meridional Eulerian-Mean Velocity (m/s)'
883        diagUnits = 'm/s             '        diagUnits = 'm/s             '
884        diagCode  = 'VVR     MR      '        diagCode  = 'VVR     MR      '
885        diagMate  = diagNum        diagMate  = diagNum
886        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
887       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
888    # endif /* ALLOW_GMREDI */
 #endif /* ALLOW_GMREDI */  
   
889  #endif /* ALLOW_EDDYPSI */  #endif /* ALLOW_EDDYPSI */
890    
891        RETURN        RETURN
Legend:
Removed from v.1.49  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.22