pkg/atm_compon_interf/cpl_diagnostics_init.F

9a4f5fc78a Jean*0001 #include "ATM_CPL_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C !ROUTINE: CPL_DIAGNOSTICS_INIT
                0006
                0007 C !INTERFACE:
                0008       SUBROUTINE CPL_DIAGNOSTICS_INIT( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize list of available diagnostics
                0012 C     for Coupled Atmospheric Component
                0013
                0014 C     !USES:
                0015       IMPLICIT NONE
                0016 #include "EEPARAMS.h"
                0017 #include "SIZE.h"
                0018 #include "CPL_PARAMS.h"
                0019
                0020 C     !INPUT/OUTPUT PARAMETERS:
                0021 C     myThid ::  my Thread Id number
                0022       INTEGER myThid
                0023 CEOP
                0024
                0025 #ifdef ALLOW_DIAGNOSTICS
                0026 C     !LOCAL VARIABLES:
                0027 C     === Local variables ===
                0028 C     diagNum   :: diagnostics number in the (long) list of available diag.
                0029 C     diagMate  :: diag. mate number in the (long) list of available diag.
                0030 C     diagName  :: local short name (8c) of a diagnostics
                0031 C     diagCode  :: local parser field with characteristics of the diagnostics
                0032 C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
                0033 C     diagUnits :: local string (16c): physical units of a diagnostic field
                0034 C     diagTitle :: local string (80c): description of field in diagnostic
                0035       INTEGER       diagNum
                0036 c     INTEGER       diagMate
                0037       CHARACTER*8   diagName
                0038       CHARACTER*16  diagCode
                0039       CHARACTER*16  diagUnits
                0040       CHARACTER*(80) diagTitle
                0041 #ifdef ALLOW_THSICE
                0042       INTEGER       numFract, numThick
                0043 #endif /* ALLOW_THSICE */
                0044
                0045 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0046
                0047 c     IF ( useDiagnotics ) THEN
                0048
                0049         diagName  = 'CPL_MxlD'
                0050         diagTitle = 'Ocean mixed-layer depth [m] (recv from Coupler)'
                0051         diagUnits = 'm               '
                0052         diagCode  = 'SM      M1      '
                0053         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0054      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0055
                0056         diagName  = 'CPL_SST '
                0057         diagTitle = 'Ocean surface temperature (recv from Coupler)'
                0058         diagUnits = 'degC            '
                0059         diagCode  = 'SM      M1      '
                0060         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0061      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0062
                0063         diagName  = 'CPL_SSS '
                0064         diagTitle = 'Ocean surface salinity    (recv from Coupler)'
                0065         diagUnits = 'g/kg            '
                0066         diagCode  = 'SM      M1      '
                0067         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0068      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0069
                0070         diagName  = 'CPL_vSq '
                0071         diagName  = 'CPL_Vel2'
                0072         diagTitle = 'Ocean surface velocity square (recv from Coupler)'
                0073         diagUnits = 'm^2/s^2         '
                0074         diagCode  = 'SM      M1      '
                0075         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0076      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0077
                0078 #ifdef ALLOW_AIM
                0079         diagName  = 'CPL_fCO2'
                0080         diagTitle ='Ocean air-sea flux of CO2 (+=down)(recv fr Coupler)'
                0081         diagUnits = 'mol/m^2/s       '
                0082         diagCode  = 'SM      M1      '
                0083         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0084      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0085 #endif /* ALLOW_AIM */
                0086
                0087 #ifdef ALLOW_THSICE
                0088         diagName  = 'CPL_icFr'
                0089         diagTitle = 'Seaice fraction (Ocean, recv from Coupler)'
                0090         diagUnits = '0-1             '
                0091         diagCode  = 'SM P    M1      '
                0092         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0093      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0094         numFract  = diagNum
                0095
                0096         diagName  = 'CPL_iceH'
                0097         diagTitle = 'Seaice thickness (Ocean, recv from Coupler)'
                0098         diagUnits = 'm               '
                0099         diagCode  = 'SM PC   M1      '
                0100         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0101      I    diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
                0102         numThick  = diagNum
                0103
                0104         diagName  = 'CPL_snwH'
                0105         diagTitle = 'Snow thickness over seaice (OCN, recv fr Coupler)'
                0106         diagUnits = 'm               '
                0107         diagCode  = 'SM PC   M1      '
                0108         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0109      I    diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
                0110
                0111         diagName  = 'CPL_Qic1'
                0112         diagTitle ='Enthalpy of seaice layer 1 (Ocean, recv fr Coupler)'
                0113         diagUnits = 'J/kg            '
                0114         diagCode  = 'SM  C   M1      '
                0115         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0116      I    diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
                0117
                0118         diagName  = 'CPL_Qic2'
                0119         diagTitle ='Enthalpy of seaice layer 2 (Ocean, recv fr Coupler)'
                0120         diagUnits = 'J/kg            '
                0121         diagCode  = 'SM  C   M1      '
                0122         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0123      I    diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
                0124 #endif /* ALLOW_THSICE */
                0125
                0126 c     ENDIF
                0127
                0128 #endif /* ALLOW_DIAGNOSTICS */
                0129
                0130       RETURN
                0131       END