/[MITgcm]/MITgcm/pkg/icefront/icefront_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm/pkg/icefront/icefront_diagnostics_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by dimitri, Wed Jan 20 23:33:45 2010 UTC revision 1.2 by dimitri, Fri Jan 22 00:51:54 2010 UTC
# Line 28  CEOP Line 28  CEOP
28  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
29  C     === Local variables ===  C     === Local variables ===
30  C     diagNum   :: diagnostics number in the (long) list of available diag.  C     diagNum   :: diagnostics number in the (long) list of available diag.
 C     diagMate  :: diag. mate number in the (long) list of available diag.  
31  C     diagName  :: local short name (8c) of a diagnostics  C     diagName  :: local short name (8c) of a diagnostics
32  C     diagCode  :: local parser field with characteristics of the diagnostics  C     diagCode  :: local parser field with characteristics of the diagnostics
33  C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT  C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
34  C     diagUnits :: local string (16c): physical units of a diagnostic field  C     diagUnits :: local string (16c): physical units of a diagnostic field
35  C     diagTitle :: local string (80c): description of field in diagnostic  C     diagTitle :: local string (80c): description of field in diagnostic
36        INTEGER       diagNum        INTEGER       diagNum
       INTEGER       diagMate  
37        CHARACTER*8   diagName        CHARACTER*8   diagName
38        CHARACTER*16  diagCode        CHARACTER*16  diagCode
39        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
# Line 43  C     diagTitle :: local string (80c): d Line 41  C     diagTitle :: local string (80c): d
41    
42  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
43    
44  c     IF ( useDiagnotics ) THEN        diagName  = 'ICFfwFlx'
45          diagTitle = 'Ice shelf fresh water flux (positive upward)'
46          diagName  = 'MYPaSta1'        diagUnits = 'kg/m^2/s        '
47          diagTitle = 'Scalar state variable 1 for package ICEFRONT'        diagCode  = 'SM      L1      '
48          diagUnits = 'unit            '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
49          diagCode  = 'SM      MR      '       I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )
50          CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
51       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )        diagName  = 'ICFhtFlx'
52          diagTitle = 'Ice shelf heat flux  (positive upward)'
53          diagName  = 'MYPaSta2'        diagUnits = 'W/m^2           '
54          diagTitle = 'Scalar state variable 2 for package ICEFRONT'        diagCode  = 'SM      L1      '
55          diagUnits = 'unit            '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
56          diagCode  = 'SM      MR      '       I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )
57          CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
58       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )        diagName  = 'ICFForcT'
59          diagTitle = 'Ice shelf forcing for theta, >0 increases theta'
60          diagName  = 'MYPaStaU'        diagUnits = 'W/m^2           '
61          diagTitle = 'uVel-type state variable for package ICEFRONT'        diagCode  = 'SM      L1      '
62          diagUnits = 'unit            '        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
63          diagCode  = 'UM      MR      '       I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )
64          diagMate  = diagNum + 2  
65          IF ( myPa_StaV_Cgrid ) diagCode(2:2) = 'U'        diagName  = 'ICFForcS'
66          CALL DIAGNOSTICS_ADDTOLIST( diagNum,        diagTitle = 'Ice shelf forcing for salt, >0 increases salt'
67       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )        diagUnits = 'g/m^2/s         '
68          diagCode  = 'SM      L1      '
69          diagName  = 'MYPaStaV'        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
70          diagTitle = 'vVel-type state variable for package ICEFRONT'       I     diagName, diagCode, diagUnits, diagTitle, 0, myThid )
         diagUnits = 'unit            '  
         diagCode  = 'VM      MR      '  
         diagMate  = diagNum  
         IF ( myPa_StaV_Cgrid ) diagCode(2:2) = 'V'  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
   
         diagName  = 'MYPaSur1'  
         diagTitle = 'Surface field 1 for package ICEFRONT'  
         diagUnits = 'unit            '  
         diagCode  = 'SM      L1      '  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  
   
         diagName  = 'MYPaSur2'  
         diagTitle = 'Surface field 2 for package ICEFRONT'  
         diagUnits = 'unit            '  
         diagCode  = 'SM      L1      '  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  
   
         diagName  = 'MYPadTdt'  
         diagTitle = ' Pot.Temp. tendency from package ICEFRONT'  
         diagUnits = 'K/s             '  
         diagCode  = 'SM      MR      '  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  
   
         diagName  = 'MYPadSdt'  
         diagTitle = ' Salinity  tendency from package ICEFRONT'  
         diagUnits = 'psu/s           '  
         diagCode  = 'SM      MR      '  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  
   
         diagName  = 'MYPadUdt'  
         diagTitle = 'u-Velocity tendency from package ICEFRONT'  
         diagUnits = 'm/s^2           '  
         diagCode  = 'UM      MR      '  
         diagMate  = diagNum + 2  
         IF ( myPa_Tend_Cgrid ) diagCode(2:2) = 'U'  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
   
         diagName  = 'MYPadVdt'  
         diagTitle = 'v-Velocity tendency from package ICEFRONT'  
         diagUnits = 'm/s^2           '  
         diagCode  = 'VM      MR      '  
         diagMate  = diagNum  
         IF ( myPa_Tend_Cgrid ) diagCode(2:2) = 'V'  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )  
   
 C--     an example of how to set "free levels option" diagnostics  
         diagName  = 'MYPaLev3'  
         diagTitle = '3-levels field for package ICEFRONT'  
         diagUnits = 'unit            '  
         diagCode  = 'SM      MX      '  
         CALL DIAGNOSTICS_ADDTOLIST( diagNum,  
      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  
         CALL DIAGNOSTICS_SETKLEV( diagName, 3, myThid )  
   
 c     ENDIF  
71    
72  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
73    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22