pkg/salt_plume/salt_plume_diagnostics_init.F

e4775240e5 Dimi*0001 #include "SALT_PLUME_OPTIONS.h"
                0002
                0003 CBOP 0
762bde3a3e Dimi*0004 C     !ROUTINE: SALT_PLUME_DIAGNOSTICS_INIT
                0005 C     !INTERFACE:
e4775240e5 Dimi*0006       SUBROUTINE SALT_PLUME_DIAGNOSTICS_INIT( myThid )
                0007
                0008 C     !DESCRIPTION:
                0009 C     Initialize list of all available diagnostics
                0010
                0011 C     !USES:
                0012       IMPLICIT NONE
                0013 #include "EEPARAMS.h"
                0014
                0015 C     !INPUT/OUTPUT PARAMETERS:
                0016 C     myThid ::  my Thread Id number
                0017       INTEGER myThid
                0018 CEOP
                0019
762bde3a3e Dimi*0020 #ifdef ALLOW_SALT_PLUME
e4775240e5 Dimi*0021 #ifdef ALLOW_DIAGNOSTICS
762bde3a3e Dimi*0022
e4775240e5 Dimi*0023 C     !LOCAL VARIABLES:
                0024 C     === Local variables ===
                0025 C     diagNum   :: diagnostics number in the (long) list of available diag.
                0026 C     diagName  :: local short name (8c) of a diagnostics
                0027 C     diagCode  :: local parser field with characteristics of the diagnostics
                0028 C              cf head of S/R DIAGNOSTICS_INIT_EARLY or DIAGNOSTICS_MAIN_INIT
                0029 C     diagUnits :: local string (16c): physical units of a diagnostic field
                0030 C     diagTitle :: local string (80c): description of field in diagnostic
                0031       INTEGER       diagNum
                0032       CHARACTER*8   diagName
                0033       CHARACTER*16  diagCode
                0034       CHARACTER*16  diagUnits
                0035       CHARACTER*(80) diagTitle
                0036
                0037       diagName  = 'PLUMEKB '
                0038       diagTitle = 'fractional plume: [0-1] (unitless)'
                0039       diagUnits = '                '
                0040       diagCode  = 'SM      MR      '
03a64b15f0 Jean*0041       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0042      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e4775240e5 Dimi*0043
2c01608d7b Gael*0044       diagName  = 'oceSPtnd'
                0045       diagTitle = 'salt tendency due to salt plume flux'
                0046      &          //' >0 increases salinity'
                0047       diagUnits = 'g/m^2/s         '
                0048       diagCode  = 'SM      MR      '
                0049       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0050      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0051
1f89baba18 Patr*0052
e4775240e5 Dimi*0053       diagName  = 'oceSPflx'
                0054       diagTitle = 'net surface Salt flux rejected into the ocean during'
                0055      &          //' freezing, (+=down),'
1d220297dd Jean*0056 c    &          //' >0 increases salinity'
e4775240e5 Dimi*0057       diagUnits = 'g/m^2/s         '
                0058       diagCode  = 'SM      U1      '
03a64b15f0 Jean*0059       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0060      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e4775240e5 Dimi*0061
                0062       diagName  = 'oceSPDep'
                0063       diagTitle = 'Salt plume depth based on density criterion (>0)'
                0064       diagUnits = 'm               '
                0065       diagCode  = 'SM      U1      '
03a64b15f0 Jean*0066       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0067      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
e4775240e5 Dimi*0068
1f89baba18 Patr*0069 #ifdef SALT_PLUME_VOLUME
                0070       diagName  = 'oceSPvol'
                0071       diagTitle = ''
                0072      &          //'surf brineVol associated w salt plume flux'
                0073       diagUnits = 'm/s             '
                0074       diagCode  = 'SM      U1      '
                0075       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0076      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0077 c
                0078       diagName  = 'PLUMEKB1'
                0079       diagTitle = 'fractional plume: [0-1] (unitless)'
                0080       diagUnits = '                '
                0081       diagCode  = 'SM      UR      '
                0082       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0083      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0084
                0085       diagName  = 'oceEPtnd'
                0086       diagTitle = 'temperature tendency due to salt plume flux'
                0087      &          //' >0 increases temperature'
                0088       diagUnits = 'W/m^2           '
                0089       diagCode  = 'SM      MR      '
                0090       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0091      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0092 c
                0093       diagName  = 'SPtd_TH '
                0094       diagTitle = 'Energy tendency associated with salt plume flux'
                0095      &          //' >0 increases temperature'
                0096       diagUnits = 'W/m^2           '
                0097       diagCode  = 'SM      MR      '
                0098       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0099      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0100 cSPtd_SLT is a duplicate of oceSPtnd
                0101       diagName  = 'SPtd_SLT'
                0102       diagTitle = 'salt tendency associated with salt plume flux'
                0103      &          //' >0 increases salinity'
                0104       diagUnits = 'g/m^2/s         '
                0105       diagCode  = 'SM      MR      '
                0106       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0107      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0108
                0109       diagName  = 'SPforcT1'
                0110       diagTitle = ''
                0111      &          //' range [1:Nr] '
                0112       diagUnits = 'unitless        '
                0113       diagCode  = 'SM      U1      '
                0114       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0115      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0116
                0117       diagName  = 'SPforcS1'
                0118       diagTitle = ''
                0119      &          //' range [1:Nr] '
                0120       diagUnits = 'unitless        '
                0121       diagCode  = 'SM      U1      '
                0122       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0123      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0124
                0125       diagName  = 'SPVsuf2k'
                0126       diagTitle = 'salt plume volume flux from surface INTO lev k'
                0127      &          //' >0 increases volume(k)'
                0128       diagUnits = 'm/s             '
                0129       diagCode  = 'SM      MR      '
                0130       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0131      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0132
                0133       diagName  = 'SPVk2m  '
                0134       diagTitle = 'vol flux LEAVING k to k-1 due to SPVsuf2k&SPVp2k'
                0135      &          //' <0 decreases volume(k)'
                0136       diagUnits = 'm/s             '
                0137       diagCode  = 'SM      MR      '
                0138       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0139      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0140 #endif /* SALT_PLUME_VOLUME */
                0141
e4775240e5 Dimi*0142 #endif /* ALLOW_DIAGNOSTICS */
762bde3a3e Dimi*0143 #endif /* ALLOW_SALT_PLUME */
e4775240e5 Dimi*0144
                0145       RETURN
                0146       END