Back to home page

MITgcm

 
 

    


File indexing completed on 2018-03-02 18:37:39 UTC

view on githubraw file Latest commit add29e06 on 2018-01-31 20:35:05 UTC
b2ea1d2979 Jean*0001 #include "ATM_PHYS_OPTIONS.h"
                0002 
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C !ROUTINE: ATM_PHYS_DIAGNOSTICS_INIT
                0006 
                0007 C !INTERFACE:
                0008       SUBROUTINE ATM_PHYS_DIAGNOSTICS_INIT( myThid )
                0009 
                0010 C     !DESCRIPTION:
                0011 C     Initialize list of all available diagnostics for pkg AtmPhys
                0012 
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "EEPARAMS.h"
                0016 #include "SIZE.h"
                0017 c#include "ATM_PHYS_PARAMS.h"
                0018 
                0019 C     !INPUT/OUTPUT PARAMETERS:
                0020 C     myThid ::  my Thread Id number
                0021       INTEGER myThid
                0022 CEOP
                0023 
                0024 #ifdef ALLOW_DIAGNOSTICS
                0025 C     !LOCAL VARIABLES:
                0026 C     === Local variables ===
                0027 C     diagNum   :: diagnostics number in the (long) list of available diag.
                0028 C     diagMate  :: diag. mate number in the (long) list of available diag.
                0029 C     diagName  :: local short name (8c) of a diagnostics
                0030 C     diagCode  :: local parser field with characteristics of the diagnostics
                0031 C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
                0032 C     diagUnits :: local string (16c): physical units of a diagnostic field
                0033 C     diagTitle :: local string (80c): description of field in diagnostic
                0034       INTEGER       diagNum
                0035       INTEGER       diagMate
                0036       CHARACTER*8   diagName
                0037       CHARACTER*16  diagCode
                0038       CHARACTER*16  diagUnits
                0039       CHARACTER*(80) diagTitle
                0040 
                0041 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0042 
                0043 c     IF ( useDiagnotics ) THEN
                0044 
                0045         diagName  = 'AtPh_SST'
                0046         diagTitle = 'Sea Surface Temperature in Atm_Phys package'
                0047         diagUnits = 'K               '
                0048         diagCode  = 'SM      L1      '
                0049         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0050      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0051 
                0052         diagName  = 'AtPhdTdt'
                0053         diagTitle = ' Pot.Temp. tendency from package ATM_PHYS'
                0054         diagUnits = 'K/s             '
                0055         diagCode  = 'SMR     MR      '
                0056         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0057      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0058 
                0059         diagName  = 'AtPhdQdt'
                0060         diagTitle = ' Spec.Humid tendency from package ATM_PHYS'
                0061         diagUnits = '1/s             '
                0062         diagCode  = 'SMR     MR      '
                0063         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0064      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0065 
                0066         diagName  = 'AtPhdUdt'
                0067         diagTitle = 'u-Velocity tendency from package ATM_PHYS'
                0068         diagUnits = 'm/s^2           '
                0069         diagCode  = 'UMR     MR      '
                0070         diagMate  = diagNum + 2
                0071         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0072      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
                0073 
                0074         diagName  = 'AtPhdVdt'
                0075         diagTitle = 'v-Velocity tendency from package ATM_PHYS'
                0076         diagUnits = 'm/s^2           '
                0077         diagCode  = 'VMR     MR      '
                0078         diagMate  = diagNum
                0079         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0080      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
                0081 
                0082         diagName  = 'AtPhDifT'
                0083         diagTitle = 'ATM_PHYS: Vertical Diffusion for Temp & Q'
                0084         diagUnits = 'unit            '
                0085         diagCode  = 'WM      LR      '
                0086         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0087      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0088         diagName  = 'AtPhDifM'
                0089         diagTitle = 'ATM_PHYS: Vertical Diffusion for Momentum'
                0090         diagUnits = 'unit            '
                0091         diagCode  = 'WM      LR      '
                0092         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0093      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0094 
                0095         diagName  = 'AtPhCnvP'
                0096         diagTitle = 'ATM_PHYS: Convective Precip'
                0097         diagUnits = 'kg/m^2/s        '
                0098         diagCode  = 'SM      L1      '
                0099         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0100      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0101         diagName  = 'AtPhLscP'
                0102         diagTitle = 'ATM_PHYS: Large-Scale Precip'
                0103         diagUnits = 'kg/m^2/s        '
                0104         diagCode  = 'SM      L1      '
                0105         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0106      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a690051f17 Jean*0107         diagName  = 'AtPhdTlc'
                0108         diagTitle = 'ATM_PHYS: lsc cond. tendency of Temp'
                0109         diagUnits = 'K/s             '
                0110         diagCode  = 'SMR     MR      '
                0111         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0112      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
b2ea1d2979 Jean*0113 
a690051f17 Jean*0114 C--     Radiation & Surface fluxes:
                0115         diagName  = 'AtPhdTrd'
                0116         diagTitle = 'ATM_PHYS: radiation tendency of Temp'
                0117         diagUnits = 'K/s             '
                0118         diagCode  = 'SMR     MR      '
                0119         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0120      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0121         diagName  = 'AtPhNSR '
                0122         diagTitle = 'ATM_PHYS: Net Shortwave Radiation (3D,+=dw)'
                0123         diagUnits = 'W/m^2           '
                0124         diagCode  = 'SM      LX      '
                0125         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0126      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0127         CALL DIAGNOSTICS_SETKLEV( diagName, Nr+1, myThid )
                0128         diagName  = 'AtPhNLR '
                0129         diagTitle = 'ATM_PHYS: Net Longwave  Radiation (3D,+=up)'
                0130         diagUnits = 'W/m^2           '
                0131         diagCode  = 'SM      LX      '
                0132         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0133      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0134         CALL DIAGNOSTICS_SETKLEV( diagName, Nr+1, myThid )
                0135         diagName  = 'AtPhDSR '
                0136         diagTitle = 'ATM_PHYS: Downward Shortwave Radiation (3D,+=dw)'
                0137         diagUnits = 'W/m^2           '
                0138         diagCode  = 'SM P    LX      '
                0139         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0140      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0141         CALL DIAGNOSTICS_SETKLEV( diagName, Nr+1, myThid )
                0142         diagName  = 'AtPhDLR '
                0143         diagTitle = 'ATM_PHYS: Downward Longwave  Radiation (3D,+=dw)'
                0144         diagUnits = 'W/m^2           '
                0145         diagCode  = 'SM P    LX      '
                0146         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0147      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0148         CALL DIAGNOSTICS_SETKLEV( diagName, Nr+1, myThid )
b2ea1d2979 Jean*0149         diagName  = 'AtPhInSR'
                0150         diagTitle = 'ATM_PHYS: Incoming Shortwave Radiation (+=dw)'
                0151         diagUnits = 'W/m^2           '
                0152         diagCode  = 'SM P    U1      '
                0153         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0154      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
31f8711b60 Jean*0155         diagName  = 'AtPhNTSR'
                0156         diagTitle = 'ATM_PHYS: Net Top Shortwave Radiation (+=dw)'
                0157         diagUnits = 'W/m^2           '
                0158         diagCode  = 'SM P    U1      '
                0159         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0160      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
b2ea1d2979 Jean*0161         diagName  = 'AtPhOLR '
                0162         diagTitle = 'ATM_PHYS: Outgoing Longwave  Radiation (+=up)'
                0163         diagUnits = 'W/m^2           '
                0164         diagCode  = 'SM P    U1      '
                0165         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0166      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
31f8711b60 Jean*0167         diagName  = 'AtPhDSSR'
                0168         diagTitle = 'ATM_PHYS: Downward Surface Shortwave Rad. (+=dw)'
                0169         diagUnits = 'W/m^2           '
                0170         diagCode  = 'SM P    U1      '
                0171         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0172      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
b2ea1d2979 Jean*0173         diagName  = 'AtPhNSSR'
                0174         diagTitle = 'ATM_PHYS: Net Surface Shortwave Radiation (+=dw)'
                0175         diagUnits = 'W/m^2           '
                0176         diagCode  = 'SM P    L1      '
                0177         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0178      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0179         diagName  = 'AtPhDSLR'
                0180         diagTitle = 'ATM_PHYS: Downward Surface Longwave Rad. (+=dw)'
                0181         diagUnits = 'W/m^2           '
                0182         diagCode  = 'SM P    L1      '
                0183         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0184      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0185         diagName  = 'AtPhUSLR'
                0186         diagTitle = 'ATM_PHYS: Upward Surface Longwave Rad. (+=up)'
                0187         diagUnits = 'W/m^2           '
                0188         diagCode  = 'SM P    L1      '
                0189         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0190      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0191         diagName  = 'AtPhSens'
                0192         diagTitle = 'ATM_PHYS: surface sensible heat flux (+=up)'
                0193         diagUnits = 'W/m^2           '
                0194         diagCode  = 'SM      L1      '
                0195         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0196      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0197         diagName  = 'AtPhEvap'
                0198         diagTitle = 'ATM_PHYS: surface evaporation'
                0199         diagUnits = 'kg/m^2/s        '
                0200         diagCode  = 'SM      L1      '
                0201         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0202      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0203         diagName  = 'AtPhTauX'
                0204         diagTitle = 'ATM_PHYS: Zonal comp of surface stress'
8fcaa83b07 Jean*0205         diagUnits = 'N/m^2           '
ca5d514420 Jean*0206         diagCode  = 'UM      L1      '
                0207         diagMate  = diagNum + 2
b2ea1d2979 Jean*0208         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
ca5d514420 Jean*0209      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
b2ea1d2979 Jean*0210         diagName  = 'AtPhTauY'
                0211         diagTitle = 'ATM_PHYS: Meridional comp of surface stress'
8fcaa83b07 Jean*0212         diagUnits = 'N/m^2           '
ca5d514420 Jean*0213         diagCode  = 'VM      L1      '
                0214         diagMate  = diagNum
b2ea1d2979 Jean*0215         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
ca5d514420 Jean*0216      I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
31f8711b60 Jean*0217 
                0218 C--     Dargan BettsMiller Convection:
                0219         diagName  = 'AtPhCAPE'
                0220         diagTitle = 'ATM_PHYS: Convective Avail. Pot. Energy: CAPE'
                0221         diagUnits = 'J/kg            '
                0222         diagCode  = 'SM      L1      '
                0223         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0224      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0225         diagName  = 'AtPhCnIn'
                0226         diagTitle = 'ATM_PHYS: Convection inhibition'
                0227         diagUnits = '1               '
                0228         diagCode  = 'SM      L1      '
                0229         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0230      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0231         diagName  = 'AtPhKlzb'
                0232         diagTitle = 'ATM_PHYS: Level of zero buoyancy (=klzb, Convect)'
                0233         diagUnits = '1               '
                0234         diagCode  = 'SM      L1      '
                0235         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0236      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0237         diagName  = 'AtPhConv'
                0238         diagTitle = 'ATM_PHYS: B.M. Convection Flag '
                0239      &            //'(0:no, 1:shallow, 2:Deep)'
                0240         diagUnits = '1               '
                0241         diagCode  = 'SM      L1      '
                0242         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0243      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0244         diagName  = 'AtPhRlxT'
                0245         diagTitle = 'ATM_PHYS: Temp. Relax. inv time-scale (Convection)'
                0246         diagUnits = '1/s             '
                0247         diagCode  = 'SM      L1      '
                0248         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0249      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0250         diagName  = 'AtPhRlxQ'
                0251         diagTitle = 'ATM_PHYS: Humid Relax. inv time-scale (Convection)'
                0252         diagUnits = '1/s             '
                0253         diagCode  = 'SM      L1      '
                0254         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0255      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0256         diagName  = 'AtPh_Trf'
                0257         diagTitle = 'ATM_PHYS: Ref Temp for convection (Moist Adiabat)'
                0258         diagUnits = 'K               '
                0259         diagCode  = 'SM      MR      '
                0260         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0261      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0262         diagName  = 'AtPh_Qrf'
                0263         diagTitle = 'ATM_PHYS: Ref Humid for convection (Saturated)'
                0264         diagUnits = 'kg/kg           '
                0265         diagCode  = 'SM      MR      '
                0266         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
b2ea1d2979 Jean*0267      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
a690051f17 Jean*0268         diagName  = 'AtPhdTcv'
                0269         diagTitle = 'ATM_PHYS: convection tendency of Temp'
                0270         diagUnits = 'K/s             '
                0271         diagCode  = 'SMR     MR      '
                0272         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0273      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
b2ea1d2979 Jean*0274 
                0275         diagName  = 'AtPhdtTg'
a690051f17 Jean*0276         diagTitle = 'ATM_PHYS: vertical diff. tendency of Temp.'
b2ea1d2979 Jean*0277         diagUnits = 'K/s             '
                0278         diagCode  = 'SMR     MR      '
                0279         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0280      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0281         diagName  = 'AtPhdtQg'
                0282         diagTitle = 'ATM_PHYS: vertical diff. tendency of spec.humid'
                0283         diagUnits = '1/s             '
                0284         diagCode  = 'SMR     MR      '
                0285         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0286      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0287         diagName  = 'AtPhDisH'
                0288         diagTitle = 'ATM_PHYS: heating from dissipation'
                0289         diagUnits = 'K/s             '
                0290         diagCode  = 'SMR     MR      '
                0291         CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0292      I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0293 
                0294 c     ENDIF
                0295 
                0296 #endif /* ALLOW_DIAGNOSTICS */
                0297 
                0298       RETURN
                0299       END