/[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.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 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_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  #ifdef ALLOW_DIAGNOSTICS        diagName  = 'AB_gU   '
831          diagName  = 'TAUXEDDY'        diagTitle = 'U momentum tendency from Adams-Bashforth'
832          diagTitle = 'Zonal Eddy Stress'        diagUnits = 'm/s^2           '
833          diagUnits = 'N/m**2          '        diagCode  = 'UUR     MR      '
834          diagCode  = 'UU P    MR      '        diagMate  = diagNum + 2
835          diagMate  = diagNum        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
836          CALL DIAGNOSTICS_ADDTOLIST( diagNum,       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
837       I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
838                  diagName  = 'AB_gV   '
839          diagName  = 'TAUYEDDY'        diagTitle = 'V momentum tendency from Adams-Bashforth'
840          diagTitle = 'Meridional Eddy Stress'        diagUnits = 'm/s^2           '
841          diagUnits = 'N/m**2          '        diagCode  = 'VVR     MR      '
842          diagCode  = 'VV P    MR      '        diagMate  = diagNum
843          diagMate  = diagNum        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
844          CALL DIAGNOSTICS_ADDTOLIST( diagNum,       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
845       I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
846  #endif /* ALLOW_DIAGNOSTICS */  #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
856          diagName  = 'TAUXEDDY'
857          diagTitle = 'Zonal Eddy Stress'
858          diagUnits = 'N/m^2           '
859          diagCode  = 'UU      LR      '
860          diagMate  = diagNum + 2
861          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
862         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
863    
864          diagName  = 'TAUYEDDY'
865          diagTitle = 'Meridional Eddy Stress'
866          diagUnits = 'N/m^2           '
867          diagCode  = 'VV      LR      '
868          diagMate  = diagNum
869          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872    # ifdef ALLOW_GMREDI
873          diagName  = 'U_EulerM'
874          diagTitle = 'Zonal Eulerian-Mean Velocity (m/s)'
875          diagUnits = 'm/s             '
876          diagCode  = 'UUR     MR      '
877          diagMate  = diagNum + 2
878          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
879         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
880    
881          diagName  = 'V_EulerM'
882          diagTitle = 'Meridional Eulerian-Mean Velocity (m/s)'
883          diagUnits = 'm/s             '
884          diagCode  = 'VVR     MR      '
885          diagMate  = diagNum
886          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
887         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
888    # endif /* ALLOW_GMREDI */
889    #endif /* ALLOW_EDDYPSI */
890    
891        RETURN        RETURN
892        END        END
Legend:
Removed from v.1.46  
changed lines
  Added in v.1.51
  ViewVC Help
Powered by ViewVC 1.1.22