pkg/land/land_diagnostics_init.F

4c4ae8300c Jean*0001 #include "LAND_OPTIONS.h"
                0002
                0003 CBOP
                0004 C     !ROUTINE: LAND_DIAGNOSTICS_INIT
                0005 C     !INTERFACE:
                0006       SUBROUTINE LAND_DIAGNOSTICS_INIT( myThid )
                0007
                0008 C     !DESCRIPTION: \bv
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE LAND_DIAGNOSTICS_INIT
                0011 C     | o Routine to initialize list of all available diagnostics
                0012 C     |   for LAND package
                0013 C     *==========================================================*
                0014 C     \ev
                0015 C     !USES:
                0016       IMPLICIT NONE
                0017
                0018 C     === Global variables ===
                0019 #include "EEPARAMS.h"
                0020 c #include "SIZE.h"
                0021 c #include "PARAMS.h"
                0022 c #include "LAND_PARAMS.h"
                0023
                0024 C     !INPUT/OUTPUT PARAMETERS:
                0025 C     === Routine arguments ===
                0026 C     myThid ::  my thread Id number
                0027       INTEGER myThid
                0028 CEOP
                0029
                0030 #ifdef ALLOW_DIAGNOSTICS
                0031 C     !LOCAL VARIABLES:
                0032 C     === Local variables ===
                0033 C     diagNum   :: diagnostics number in the (long) list of available diag.
                0034 C     diagName  :: local short name (8c) of a diagnostics
                0035 C     diagCode  :: local parser field with characteristics of the diagnostics
                0036 C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
                0037 C     diagUnits :: local string (16c): physical units of a diagnostic field
                0038 C     diagTitle :: local string (80c): description of field in diagnostic
                0039       INTEGER       diagNum
                0040       CHARACTER*8   diagName
                0041       CHARACTER*16  diagCode
                0042       CHARACTER*16  diagUnits
                0043       CHARACTER*(80) diagTitle
                0044
                0045 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0046
                0047 c     IF ( useDiagnostics ) THEN
                0048
                0049 C--    Land state variables:
                0050         diagName  = 'GrdSurfT'
                0051         diagTitle = 'Surface Temperature over land'
                0052         diagUnits = 'degC            '
1d6f9f4568 Jean*0053         diagCode  = 'SM      Lg      '
03a64b15f0 Jean*0054         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0055      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0056
                0057         diagName  = 'GrdTemp '
                0058         diagTitle = 'Ground Temperature at each level'
                0059         diagUnits = 'degC            '
                0060         diagCode  = 'SM      MG      '
03a64b15f0 Jean*0061         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0062      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0063
                0064         diagName  = 'GrdEnth '
                0065         diagTitle = 'Ground Enthalpy at each level'
                0066         diagUnits = 'J/m3            '
                0067         diagCode  = 'SM      MG      '
03a64b15f0 Jean*0068         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0069      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0070
                0071         diagName  = 'GrdWater'
1d6f9f4568 Jean*0072         diagTitle =
4c4ae8300c Jean*0073      &       'Ground Water (vs Field Capacity) Fraction at each level'
                0074         diagUnits = '0-1             '
                0075         diagCode  = 'SM P    MG      '
03a64b15f0 Jean*0076         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0077      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0078
                0079 c       cdiag(160) = 'SNOW    '
                0080 c       tdiag(160) = 'Snow Depth (mm water equivalent)'
                0081 c       udiag(160) = 'mm water equiv  '
                0082 c       gdiag(160) = 'SM      L1      '
                0083 c---
                0084         diagName  = 'LdSnowH '
                0085         diagTitle = 'Snow Thickness over land'
                0086         diagUnits = 'm               '
1d6f9f4568 Jean*0087         diagCode  = 'SM P    Lg      '
03a64b15f0 Jean*0088         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0089      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0090
                0091         diagName  = 'LdSnwAge'
                0092         diagTitle = 'Snow Age over land'
                0093         diagUnits = 's               '
1d6f9f4568 Jean*0094         diagCode  = 'SM P    Lg      '
03a64b15f0 Jean*0095         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0096      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0097
                0098 C--    Fluxes:
                0099
                0100         diagName  = 'RUNOFF  '
                0101         diagTitle = 'Run-Off per surface unit'
fe3429c1e0 Davi*0102         diagUnits = 'kg/m2/s         '
4c4ae8300c Jean*0103         diagCode  = 'SM      L1      '
03a64b15f0 Jean*0104         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0105      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0106
                0107         diagName  = 'EnRunOff'
                0108         diagTitle = 'Energy flux associated with run-Off'
                0109         diagUnits = 'W/m^2           '
                0110         diagCode  = 'SM      L1      '
03a64b15f0 Jean*0111         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0112      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0113
                0114         diagName  = 'landHFlx'
                0115         diagTitle = 'net surface downward Heat flux over land'
                0116         diagUnits = 'W/m^2           '
1d6f9f4568 Jean*0117         diagCode  = 'SM      Lg      '
03a64b15f0 Jean*0118         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0119      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0120
                0121         diagName  = 'landPmE '
                0122         diagTitle = 'Precipitation minus Evaporation over land'
                0123         diagUnits = 'kg/m^2/s        '
1d6f9f4568 Jean*0124         diagCode  = 'SM      Lg      '
03a64b15f0 Jean*0125         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0126      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0127
                0128         diagName  = 'ldEnFxPr'
1d6f9f4568 Jean*0129         diagTitle =
4c4ae8300c Jean*0130      &    'Energy flux (over land) associated with Precip (snow,rain)'
                0131         diagUnits = 'W/m^2           '
1d6f9f4568 Jean*0132         diagCode  = 'SM      Lg      '
03a64b15f0 Jean*0133         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0134      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
4c4ae8300c Jean*0135
                0136 c     ENDIF
                0137
                0138 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0139
                0140 #endif /* ALLOW_DIAGNOSTICS */
                0141
                0142       RETURN
                0143       END