pkg/shap_filt/shap_filt_diagnostics_init.F

64e64319fa Jean*0001 #include "SHAP_FILT_OPTIONS.h"
                0002
                0003 CBOP
                0004 C     !ROUTINE: SHAP_FILT_DIAGNOSTICS_INIT
                0005 C     !INTERFACE:
                0006       SUBROUTINE SHAP_FILT_DIAGNOSTICS_INIT( myThid )
                0007
                0008 C     !DESCRIPTION: \bv
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE SHAP_FILT_DIAGNOSTICS_INIT
                0011 C     | o Routine to initialize list of all available diagnostics
                0012 C     |   for SHAP_FILT package
                0013 C     *==========================================================*
                0014 C     \ev
                0015 C     !USES:
                0016       IMPLICIT NONE
                0017
                0018 C     === Global variables ===
41bf0a383d Jean*0019 #include "SIZE.h"
64e64319fa Jean*0020 #include "EEPARAMS.h"
41bf0a383d Jean*0021 #include "PARAMS.h"
64e64319fa Jean*0022
                0023 C     !INPUT/OUTPUT PARAMETERS:
                0024 C     === Routine arguments ===
                0025 C     myThid ::  my thread Id number
                0026       INTEGER myThid
                0027 CEOP
                0028
                0029 #ifdef ALLOW_DIAGNOSTICS
41bf0a383d Jean*0030 C     ! FUNCTIONS:
                0031       CHARACTER*16 DIAGS_MK_UNITS
                0032       EXTERNAL     DIAGS_MK_UNITS
                0033
64e64319fa Jean*0034 C     !LOCAL VARIABLES:
                0035 C     === Local variables ===
                0036 C     diagNum   :: diagnostics number in the (long) list of available diag.
a51fbd3807 Jean*0037 C     diagMate  :: diag. mate number in the (long) list of available diag.
64e64319fa Jean*0038 C     diagName  :: local short name (8c) of a diagnostics
                0039 C     diagCode  :: local parser field with characteristics of the diagnostics
                0040 C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
                0041 C     diagUnits :: local string (16c): physical units of a diagnostic field
                0042 C     diagTitle :: local string (80c): description of field in diagnostic
                0043       INTEGER       diagNum
a51fbd3807 Jean*0044       INTEGER       diagMate
64e64319fa Jean*0045       CHARACTER*8   diagName
                0046       CHARACTER*16  diagCode
                0047       CHARACTER*16  diagUnits
                0048       CHARACTER*(80) diagTitle
41bf0a383d Jean*0049       CHARACTER*10  trUnits
64e64319fa Jean*0050
                0051 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0052
                0053 c     IF ( useDiagnostics ) THEN
                0054
                0055 C-    Shapiro Filter tendency:
                0056
                0057         diagName  = 'SHAP_dT '
                0058         diagTitle = 'Temperature Tendency due to Shapiro Filter'
                0059         diagUnits = 'K/s             '
a51fbd3807 Jean*0060         diagCode  = 'SMR     MR      '
                0061         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0062      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
64e64319fa Jean*0063
41bf0a383d Jean*0064         trUnits = 'kg/kg'
                0065         IF ( fluidIsWater .OR. useAIM ) trUnits = 'g/kg'
64e64319fa Jean*0066         diagName  = 'SHAP_dS '
                0067         diagTitle = 'Specific Humidity Tendency due to Shapiro Filter'
41bf0a383d Jean*0068         diagUnits =  DIAGS_MK_UNITS( trUnits//'/s', myThid )
a51fbd3807 Jean*0069         diagCode  = 'SMR     MR      '
                0070         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0071      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
64e64319fa Jean*0072
                0073         diagName  = 'SHAP_dU '
                0074         diagTitle = 'Zonal Wind Tendency due to Shapiro Filter'
                0075         diagUnits = 'm/s^2           '
a51fbd3807 Jean*0076         diagCode  = 'UUR     MR      '
                0077         diagMate  = diagNum + 2
                0078         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0079      I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
64e64319fa Jean*0080
                0081         diagName  = 'SHAP_dV '
                0082         diagTitle = 'Meridional Wind Tendency due to Shapiro Filter'
                0083         diagUnits = 'm/s^2           '
a51fbd3807 Jean*0084         diagCode  = 'VVR     MR      '
                0085         diagMate  = diagNum
                0086         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0087      I    diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
64e64319fa Jean*0088
73e3454c50 Davi*0089         diagName  = 'SHAP_dKE'
a51fbd3807 Jean*0090         diagTitle =
73e3454c50 Davi*0091      &  'Kinetic Energy Tendency due to Shapiro Filter (Mass Weighted)'
                0092         diagUnits = 'm^2/s^3         '
a51fbd3807 Jean*0093         diagCode  = 'SMr     MR      '
                0094         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0095      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
73e3454c50 Davi*0096
64e64319fa Jean*0097 c     ENDIF
                0098
                0099 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0100
                0101 #endif /* ALLOW_DIAGNOSTICS */
                0102
                0103       RETURN
                0104       END