pkg/ocn_compon_interf/cpl_diagnostics_init.F

9a4f5fc78a Jean*0001 #include "OCN_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 Oceanic 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       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_SLPr'
                0050         diagTitle = 'Atmos Sea-Level pressure anomaly (recv fr Coupler)'
                0051         diagUnits = 'Pa              '
                0052         diagCode  = 'SM      M1      '
                0053         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0054      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0055
                0056         diagName  = 'CPL_Qnet'
                0057         diagTitle = 'Net surface heat flux (+=up) (recv from Coupler)'
                0058         diagUnits = 'W/m^2           '
                0059         diagCode  = 'SM      M1      '
                0060         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0061      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0062
                0063         diagName  = 'CPL_Qsw '
                0064         diagTitle = 'Net shortwave heat flux (+=up)(recv from Coupler)'
                0065         diagUnits = 'W/m^2           '
                0066         diagCode  = 'SM      M1      '
                0067         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0068      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0069
                0070         diagName  = 'CPL_tauX'
                0071         diagTitle = 'Zonal  surface wind-stress (recv from Coupler)'
                0072         diagUnits = 'N/m^2           '
                0073         diagCode  = 'UM      M1      '
                0074         diagMate  = diagNum + 2
                0075         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0076      I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
                0077
                0078         diagName  = 'CPL_tauY'
                0079         diagTitle = 'Merid. surface wind-stress (recv from Coupler)'
                0080         diagUnits = 'N/m^2           '
                0081         diagMate  = diagNum
                0082         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0083      I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
                0084
                0085         diagName  = 'CPL_FWat'
                0086         diagTitle = 'Surface fresh-water flux (+=up)(recv from Coupler)'
                0087         diagUnits = 'kg/m^2/s        '
                0088         diagCode  = 'SM      M1      '
                0089         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0090      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0091
                0092         diagName  = 'CPL_sFlx'
                0093         diagTitle = 'Salt flux from seaice comp (+=up)(recv fr Coupler)'
                0094         diagUnits = 'g/m^2/s         '
                0095         diagCode  = 'SM      M1      '
                0096         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0097      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0098
                0099         diagName  = 'CPL_iceM'
                0100         diagTitle = 'Seaice mass       (recv from Coupler)'
                0101         diagUnits = 'kg/m^2          '
                0102         diagCode  = 'SM      M1      '
                0103         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0104      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0105
                0106 #ifdef ALLOW_SALT_PLUME
                0107         diagName  = 'CPL_sPlm'
                0108         diagTitle = 'Salt-plume flux (+=down) (recv from Coupler)'
                0109         diagUnits = 'g/m^2/s         '
                0110         diagCode  = 'SM      M1      '
                0111         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0112      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0113 #endif /* ALLOW_SALT_PLUME */
                0114
                0115 #ifdef ALLOW_DIC
                0116         diagName  = 'CPL_RnOf'
                0117         diagTitle = 'river Run-Off (for DIC, +=down)(recv from Coupler)'
                0118         diagUnits = 'kg/m^2/s        '
                0119         diagCode  = 'SM      M1      '
                0120         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0121      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0122
                0123         diagName  = 'CPL_aCO2'
                0124         diagTitle = 'Atmospheric CO2 (parts by vol)(recv from Coupler)'
                0125         diagUnits = 'ppm             '
                0126         diagCode  = 'SM      M1      '
                0127         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0128      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0129
                0130         diagName  = 'CPL_wSpd'
                0131         diagTitle = 'Atmospheric surf. wind speed (recv from Coupler)'
                0132         diagUnits = 'm/s             '
                0133         diagCode  = 'SM      M1      '
                0134         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0135      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0136 #endif /* ALLOW_DIC */
                0137
                0138 #if (defined ALLOW_DIC) || (defined ALLOW_THSICE)
                0139         diagName  = 'CPL_icFr'
                0140         diagTitle = 'Seaice fraction (Atmos, recv from Coupler)'
                0141         diagUnits = '0-1             '
                0142         diagCode  = 'SM P    M1      '
                0143         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0144      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0145 #endif /* ALLOW_DIC or ALLOW_THSICE */
                0146 #ifdef ALLOW_THSICE
                0147         numFract  = diagNum
                0148
                0149         diagName  = 'CPL_iceH'
                0150         diagTitle = 'Seaice thickness (Atmos, recv from Coupler)'
                0151         diagUnits = 'm               '
                0152         diagCode  = 'SM PC   M1      '
                0153         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0154      I    diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
                0155         numThick  = diagNum
                0156
                0157         diagName  = 'CPL_snwH'
                0158         diagTitle = 'Snow thickness over seaice (ATM, recv fr Coupler)'
                0159         diagUnits = 'm               '
                0160         diagCode  = 'SM PC   M1      '
                0161         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0162      I    diagName, diagCode, diagUnits, diagTitle, numFract, myThid )
                0163
                0164         diagName  = 'CPL_Qic1'
                0165         diagTitle ='Enthalpy of seaice layer 1 (Atmos, recv fr Coupler)'
                0166         diagUnits = 'J/kg            '
                0167         diagCode  = 'SM  C   M1      '
                0168         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0169      I    diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
                0170
                0171         diagName  = 'CPL_Qic2'
                0172         diagTitle ='Enthalpy of seaice layer 2 (Atmos, recv fr Coupler)'
                0173         diagUnits = 'J/kg            '
                0174         diagCode  = 'SM  C   M1      '
                0175         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0176      I    diagName, diagCode, diagUnits, diagTitle, numThick, myThid )
                0177 #endif /* ALLOW_THSICE */
                0178
                0179 c     ENDIF
                0180
                0181 #endif /* ALLOW_DIAGNOSTICS */
                0182
                0183       RETURN
                0184       END