bling/pkg/bling_diagnostics_init.F

C $Header:  $
C $Name:  $

#include "GCHEM_OPTIONS.h"

CBOP
      subroutine BLING_DIAGNOSTICS_INIT(myThid )

C     =================================================================
C     | subroutine bling_diagnostics_init
C     | o define diagnostics for bling pkg
C     =================================================================

      implicit none
      
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"

C  myThid               :: thread number
      INTEGER myThid
CEOP

#ifdef ALLOW_DIAGNOSTICS

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

C     Define diagnostics Names :

      diagName  = 'BLGBIOA '
      diagTitle = 'Tendency of carbon due to bio prod (mol C/m3/s)'
      diagUnits = 'mol C/m3/sec    '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGREMI '
      diagTitle = 'Tendency of carbon due to bio remin (mol C/m3/s)'
      diagUnits = 'mol C/m3/sec    '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGCARB '
      diagTitle = 'Carbonate tendency due to biology (mol C/m3/s)'
      diagUnits = 'mol C/m3/sec    '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGPSM  '
      diagTitle = 'Small phytoplankton biomass (mol C/m3/s)'
      diagUnits = 'mol C/m3        '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGPLG  '
      diagTitle = 'Large phytoplankton biomass (mol C/m3/s)'
      diagUnits = 'mol C/m3        '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGCHL  '
      diagTitle = 'Chlorophyll concentration'
      diagUnits = 'mol chl/m3      '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

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

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

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

      diagName  = 'BLGPCO2 '
      diagTitle = 'Surface ocean partial pressure of CO2'
      diagUnits = 'atm             '
      diagCode  = 'SM P    M1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGPH3D '
      diagTitle = 'pH'
      diagUnits = 'dimensionless   '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'BLGOMAR '
      diagTitle = 'Aragonite saturation'
      diagUnits = 'dimensionless   '
      diagCode  = 'SMRP    MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )

#endif /* ALLOW_DIAGNOSTICS */

      RETURN
      END
1	C $Header: $
2	C $Name: $
3
4	#include "GCHEM_OPTIONS.h"
5
6	CBOP
7	subroutine BLING_DIAGNOSTICS_INIT(myThid )
8
9	C =================================================================
10	C \| subroutine bling_diagnostics_init
11	C \| o define diagnostics for bling pkg
12	C =================================================================
13
14	implicit none
15
16	C === Global variables ===
17	#include "SIZE.h"
18	#include "EEPARAMS.h"
19
20	C myThid :: thread number
21	INTEGER myThid
22	CEOP
23
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 Define diagnostics Names :
34
35	diagName = 'BLGBIOA '
36	diagTitle = 'Tendency of carbon due to bio prod (mol C/m3/s)'
37	diagUnits = 'mol C/m3/sec '
38	diagCode = 'SMRP MR '
39	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
40	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
41
42	diagName = 'BLGREMI '
43	diagTitle = 'Tendency of carbon due to bio remin (mol C/m3/s)'
44	diagUnits = 'mol C/m3/sec '
45	diagCode = 'SMRP MR '
46	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
47	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
48
49	diagName = 'BLGCARB '
50	diagTitle = 'Carbonate tendency due to biology (mol C/m3/s)'
51	diagUnits = 'mol C/m3/sec '
52	diagCode = 'SMRP MR '
53	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
54	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
55
56	diagName = 'BLGPSM '
57	diagTitle = 'Small phytoplankton biomass (mol C/m3/s)'
58	diagUnits = 'mol C/m3 '
59	diagCode = 'SMRP MR '
60	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
61	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
62
63	diagName = 'BLGPLG '
64	diagTitle = 'Large phytoplankton biomass (mol C/m3/s)'
65	diagUnits = 'mol C/m3 '
66	diagCode = 'SMRP MR '
67	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
68	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
69
70	diagName = 'BLGCHL '
71	diagTitle = 'Chlorophyll concentration'
72	diagUnits = 'mol chl/m3 '
73	diagCode = 'SMRP MR '
74	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
75	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
76
77	diagName = 'BLGTFLX '
78	diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
79	diagUnits = 'mol C/m3/sec '
80	diagCode = 'SM P L1 '
81	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
82	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
83
84	diagName = 'BLGOFLX '
85	diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
86	diagUnits = 'mol O2/m3/sec '
87	diagCode = 'SM P L1 '
88	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
89	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
90
91	diagName = 'BLGCFLX '
92	diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)'
93	diagUnits = 'mol/m2/sec '
94	diagCode = 'SM P L1 '
95	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
96	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
97
98	diagName = 'BLGPCO2 '
99	diagTitle = 'Surface ocean partial pressure of CO2'
100	diagUnits = 'atm '
101	diagCode = 'SM P M1 '
102	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
103	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
104
105	diagName = 'BLGPH3D '
106	diagTitle = 'pH'
107	diagUnits = 'dimensionless '
108	diagCode = 'SMRP MR '
109	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
110	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
111
112	diagName = 'BLGOMAR '
113	diagTitle = 'Aragonite saturation'
114	diagUnits = 'dimensionless '
115	diagCode = 'SMRP MR '
116	CALL DIAGNOSTICS_ADDTOLIST( diagNum,
117	I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
118
119	#endif /* ALLOW_DIAGNOSTICS */
120
121	RETURN
122	END