pkg/gchem/gchem_diagnostics_init.F

C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_diagnostics_init.F,v 1.2 2007/07/31 23:01:58 ce107 Exp $
C $Name:  $

#include "GCHEM_OPTIONS.h"

C !INTERFACE: ==========================================================
      SUBROUTINE GCHEM_DIAGNOSTICS_INIT(myThid )

C !DESCRIPTION:
C define diagnostics for gchem package
C experiment

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"

C !INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid
CEOP

#ifdef ALLOW_GCHEM
#ifdef ALLOW_DIAGNOSTICS

C     !LOCAL VARIABLES:
      INTEGER       diagNum
      CHARACTER*8   diagName
      CHARACTER*16  diagCode
      CHARACTER*16  diagUnits
      CHARACTER*(80) diagTitle

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
C     Define diagnostics Names :

      diagName  = 'DICBIOA '
      diagTitle = 'Biological Productivity (mol/m3/s)'
      diagUnits = 'mol/m3/sec      '
      diagCode  = 'SM P    MR      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICCARB '
      diagTitle = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
      diagUnits = 'mol eq/m3/sec   '
      diagCode  = 'SM P    MR      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICTFLX '
      diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
      diagUnits = 'mol/m3/sec      '
      diagCode  = 'SM P    L1      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICOFLX '
      diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
      diagUnits = 'mol/m3/sec      '
      diagCode  = 'SM P    L1      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICCFLX '
      diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)'
      diagUnits = 'mol/m2/sec      '
      diagCode  = 'SM P    L1      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICPCO2 '
      diagTitle = 'Partial Pressure of CO2 (atm)'
      diagUnits = 'atm             '
      diagCode  = 'SM P    M1      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

      diagName  = 'DICPHAV '
      diagTitle = 'pH (dimensionless)'
      diagUnits = 'dimensionless   '
      diagCode  = 'SM P    M1      '
      CALL DIAGNOSTICS_ADD2LIST( diagNum,
     I          diagName, diagCode, diagUnits, diagTitle, myThid )

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

#endif /* ALLOW_DIAGNOSTICS */
#endif /* ALLOW_GCHEM */

      RETURN
      END
1	jahn	1.1	C $Header: /u/gcmpack/MITgcm/pkg/gchem/gchem_diagnostics_init.F,v 1.2 2007/07/31 23:01:58 ce107 Exp $
2			C $Name: $
3
4			#include "GCHEM_OPTIONS.h"
5
6			C !INTERFACE: ==========================================================
7			SUBROUTINE GCHEM_DIAGNOSTICS_INIT(myThid )
8
9			C !DESCRIPTION:
10			C define diagnostics for gchem package
11			C experiment
12
13			C !USES: ===============================================================
14			IMPLICIT NONE
15			#include "SIZE.h"
16			#include "EEPARAMS.h"
17
18			C !INPUT PARAMETERS: ===================================================
19			C myThid :: thread number
20			INTEGER myThid
21			CEOP
22
23			#ifdef ALLOW_GCHEM
24			#ifdef ALLOW_DIAGNOSTICS
25
26			C !LOCAL VARIABLES:
27			INTEGER diagNum
28			CHARACTER*8 diagName
29			CHARACTER*16 diagCode
30			CHARACTER*16 diagUnits
31			CHARACTER*(80) diagTitle
32
33			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
34			C Define diagnostics Names :
35
36			diagName = 'DICBIOA '
37			diagTitle = 'Biological Productivity (mol/m3/s)'
38			diagUnits = 'mol/m3/sec '
39			diagCode = 'SM P MR '
40			CALL DIAGNOSTICS_ADD2LIST( diagNum,
41			I diagName, diagCode, diagUnits, diagTitle, myThid )
42
43			diagName = 'DICCARB '
44			diagTitle = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
45			diagUnits = 'mol eq/m3/sec '
46			diagCode = 'SM P MR '
47			CALL DIAGNOSTICS_ADD2LIST( diagNum,
48			I diagName, diagCode, diagUnits, diagTitle, myThid )
49
50			diagName = 'DICTFLX '
51			diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
52			diagUnits = 'mol/m3/sec '
53			diagCode = 'SM P L1 '
54			CALL DIAGNOSTICS_ADD2LIST( diagNum,
55			I diagName, diagCode, diagUnits, diagTitle, myThid )
56
57			diagName = 'DICOFLX '
58			diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
59			diagUnits = 'mol/m3/sec '
60			diagCode = 'SM P L1 '
61			CALL DIAGNOSTICS_ADD2LIST( diagNum,
62			I diagName, diagCode, diagUnits, diagTitle, myThid )
63
64			diagName = 'DICCFLX '
65			diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)'
66			diagUnits = 'mol/m2/sec '
67			diagCode = 'SM P L1 '
68			CALL DIAGNOSTICS_ADD2LIST( diagNum,
69			I diagName, diagCode, diagUnits, diagTitle, myThid )
70
71			diagName = 'DICPCO2 '
72			diagTitle = 'Partial Pressure of CO2 (atm)'
73			diagUnits = 'atm '
74			diagCode = 'SM P M1 '
75			CALL DIAGNOSTICS_ADD2LIST( diagNum,
76			I diagName, diagCode, diagUnits, diagTitle, myThid )
77
78			diagName = 'DICPHAV '
79			diagTitle = 'pH (dimensionless)'
80			diagUnits = 'dimensionless '
81			diagCode = 'SM P M1 '
82			CALL DIAGNOSTICS_ADD2LIST( diagNum,
83			I diagName, diagCode, diagUnits, diagTitle, myThid )
84
85			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
86
87			#endif /* ALLOW_DIAGNOSTICS */
88			#endif /* ALLOW_GCHEM */
89
90			RETURN
91			END