pkg/dic/dic_diagnostics_init.F

a1d0e455fd Hann*0001 #include "DIC_OPTIONS.h"
add5eda49e Oliv*0002
                0003 C !INTERFACE: ==========================================================
                0004       SUBROUTINE DIC_DIAGNOSTICS_INIT(myThid )
                0005
                0006 C !DESCRIPTION:
                0007 C define diagnostics for dic package
                0008 C experiment
                0009
                0010 C !USES: ===============================================================
                0011       IMPLICIT NONE
                0012 #include "SIZE.h"
                0013 #include "EEPARAMS.h"
                0014
                0015 C !INPUT PARAMETERS: ===================================================
                0016 C  myThid               :: thread number
                0017       INTEGER myThid
                0018 CEOP
                0019
                0020 #ifdef ALLOW_DIAGNOSTICS
                0021
                0022 C     !LOCAL VARIABLES:
                0023       INTEGER       diagNum
                0024       CHARACTER*8   diagName
                0025       CHARACTER*16  diagCode
                0026       CHARACTER*16  diagUnits
                0027       CHARACTER*(80) diagTitle
                0028
                0029 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0030 C     Define diagnostics Names :
                0031
                0032       diagName  = 'DICBIOA '
                0033       diagTitle = 'Biological Productivity (mol/m3/s)'
                0034       diagUnits = 'mol/m3/sec      '
                0035       diagCode  = 'SMRP    MR      '
                0036       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0037      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0038
                0039       diagName  = 'DICCARB '
                0040       diagTitle = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
                0041       diagUnits = 'mol eq/m3/sec   '
                0042       diagCode  = 'SMRP    MR      '
                0043       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0044      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0045
                0046       diagName  = 'DICTFLX '
                0047       diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
                0048       diagUnits = 'mol/m3/sec      '
                0049       diagCode  = 'SM P    L1      '
                0050       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0051      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0052
                0053       diagName  = 'DICOFLX '
                0054       diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
                0055       diagUnits = 'mol/m3/sec      '
                0056       diagCode  = 'SM P    L1      '
                0057       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0058      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0059
                0060       diagName  = 'DICCFLX '
                0061       diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)'
                0062       diagUnits = 'mol/m2/sec      '
                0063       diagCode  = 'SM P    L1      '
                0064       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0065      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0066
                0067       diagName  = 'DICPCO2 '
                0068       diagTitle = 'Partial Pressure of CO2 (atm)'
                0069       diagUnits = 'atm             '
                0070       diagCode  = 'SM P    M1      '
                0071       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0072      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0073
                0074       diagName  = 'DICPHAV '
                0075       diagTitle = 'pH (dimensionless)'
                0076       diagUnits = 'dimensionless   '
                0077       diagCode  = 'SM P    M1      '
                0078       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0079      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0080
2e3e8c330d Jona*0081 #ifdef DIC_CALCITE_SAT
                0082        diagName  = 'DIC3DSIT'
                0083        diagTitle = 'Three dimensional silicate concentration (mol/m3)'
                0084        diagUnits = 'mol/m3          '
                0085        diagCode  = 'SMRP    MR      '
                0086        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
a1d0e455fd Hann*0087      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0088
2e3e8c330d Jona*0089        diagName  = 'OMEGAC  '
                0090        diagTitle = 'Carbonate saturation'
                0091        diagUnits = 'mol eq/m3/s     '
                0092        diagCode  = 'SMRP    MR      '
                0093        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
a1d0e455fd Hann*0094      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0095
2e3e8c330d Jona*0096        diagName  = 'DIC3DPH '
                0097        diagTitle = 'Three dimensional pH (dimensionless)'
                0098        diagUnits = 'dimensionless   '
                0099        diagCode  = 'SM P    MR      '
                0100        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
a1d0e455fd Hann*0101      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0102
2e3e8c330d Jona*0103        diagName  = 'DIC3DPCO'
                0104        diagTitle = 'Three dimensional CO2 partial pressure (atm)'
                0105        diagUnits = 'atm             '
                0106        diagCode  = 'SMRP    MR      '
                0107        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
a1d0e455fd Hann*0108      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0109
2e3e8c330d Jona*0110        diagName  = 'DIC3DCO3'
                0111        diagTitle = 'Three dimensional carbonate concentration (mol/m3)'
                0112        diagUnits = 'mol/m3          '
                0113        diagCode  = 'SMRP    MR      '
                0114        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
a1d0e455fd Hann*0115      I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
2e3e8c330d Jona*0116 #endif /* DIC_CALCITE_SAT */
a1d0e455fd Hann*0117
add5eda49e Oliv*0118 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0119
                0120 #endif /* ALLOW_DIAGNOSTICS */
                0121
                0122       RETURN
                0123       END