/[MITgcm]/MITgcm_contrib/dcarroll/highres_darwin/code/dic_diagnostics_init.F
ViewVC logotype

Contents of /MITgcm_contrib/dcarroll/highres_darwin/code/dic_diagnostics_init.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Sun Sep 22 21:23:46 2019 UTC (5 years, 10 months ago) by dcarroll
Branch: MAIN
CVS Tags: HEAD
Initial check in of high resolution Darwin simulation code

1 C $Header: /u/gcmpack/MITgcm_contrib/ecco_darwin/v4_llc270/code_darwin/dic_diagnostics_init.F,v 1.4 2019/08/30 00:13:34 dcarroll Exp $
2 C $Name: $
3
4 #include "GCHEM_OPTIONS.h"
5
6 C !INTERFACE: ==========================================================
7 SUBROUTINE DIC_DIAGNOSTICS_INIT(myThid )
8
9 C !DESCRIPTION:
10 C define diagnostics for dic 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_DIAGNOSTICS
24
25 C !LOCAL VARIABLES:
26 INTEGER diagNum
27 CHARACTER*8 diagName
28 CHARACTER*16 diagCode
29 CHARACTER*16 diagUnits
30 CHARACTER*(80) diagTitle
31
32 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
33 C Define diagnostics Names :
34
35 diagName = 'DICBIOA '
36 diagTitle = 'Biological Productivity (mol/m3/s)'
37 diagUnits = 'mol/m3/sec '
38 diagCode = 'SMRP MR '
39 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
40 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
41
42 diagName = 'DICCARB '
43 diagTitle = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
44 diagUnits = 'mol eq/m3/sec '
45 diagCode = 'SMRP MR '
46 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
47 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
48
49 diagName = 'DICTFLX '
50 diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
51 diagUnits = 'mol/m3/sec '
52 diagCode = 'SM P L1 '
53 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
54 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
55
56 diagName = 'DICOFLX '
57 diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
58 diagUnits = 'mol/m3/sec '
59 diagCode = 'SM P L1 '
60 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
61 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
62
63 diagName = 'DICCFLX '
64 diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)'
65 diagUnits = 'mol/m2/sec '
66 diagCode = 'SM P L1 '
67 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
68 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
69
70 diagName = 'DICFGCO2 '
71 diagTitle = 'Fugacity of CO2 (atm)'
72 diagUnits = 'atm '
73 diagCode = 'SM P L1 '
74 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
75 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
76
77 diagName = 'DICPCO2 '
78 diagTitle = 'Partial Pressure of CO2 (atm)'
79 diagUnits = 'atm '
80 diagCode = 'SM P M1 '
81 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
82 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
83
84 diagName = 'DICPHAV '
85 diagTitle = 'pH (dimensionless)'
86 diagUnits = 'dimensionless '
87 diagCode = 'SM P M1 '
88 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
89 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
90
91 diagName = 'DICSFLX '
92 diagTitle = 'DIC flux from sediments (mol/m2/s)'
93 diagUnits = 'mol/m2/sec '
94 diagCode = 'SM P M1 '
95 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
96 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
97
98 diagName = 'ALKSFLX '
99 diagTitle = 'ALK flux from sediments (mol/m2/s)'
100 diagUnits = 'mol/m2/sec '
101 diagCode = 'SM P M1 '
102 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
103 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
104
105 diagName = 'CO3SW '
106 diagTitle = 'CO3 at seafloor (mol/m3)'
107 diagUnits = 'mol/m3 '
108 diagCode = 'SM P M1 '
109 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
110 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
111
112 diagName = 'CO3SED '
113 diagTitle = 'CO3 in sediment (mol/m3)'
114 diagUnits = 'mol/m3 '
115 diagCode = 'SM P M1 '
116 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
117 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
118
119 diagName = 'OMEGAC '
120 diagTitle = 'calcite saturation state (dimensionless)'
121 diagUnits = 'dimensionless '
122 diagCode = 'SM P MR '
123 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
124 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
125
126 diagName = 'DISSC '
127 diagTitle = 'calcite dissolution rate (mol/m3/s)'
128 diagUnits = 'mol/m3/s '
129 diagCode = 'SM P MR '
130 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
131 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
132
133 diagName = 'DCO2FLX '
134 diagTitle = 'delta flux of CO2 - air-sea exch (mol/m2/s)'
135 diagUnits = 'mol/m2/s '
136 diagCode = 'SM P M1 '
137 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
138 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
139
140 diagName = 'DCO2FLXT'
141 diagTitle = 'delta CO2 flux due to temperature (mol/m2/s)'
142 diagUnits = 'mol/m2/s '
143 diagCode = 'SM P M1 '
144 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
145 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
146
147 diagName = 'DCO2FLXS'
148 diagTitle = 'delta CO2 flux due to salinity (mol/m2/s)'
149 diagUnits = 'mol/m2/s '
150 diagCode = 'SM P M1 '
151 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
152 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
153
154 diagName = 'DCO2FLXA'
155 diagTitle = 'delta CO2 flux due to alkalinity (mol/m2/s)'
156 diagUnits = 'mol/m2/s '
157 diagCode = 'SM P M1 '
158 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
159 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
160
161 C diagName = 'DCO2FLXD'
162 C diagTitle = 'delta CO2 flux due to DIC (mol/m2/s)'
163 C diagUnits = 'mol/m2/s '
164 C diagCode = 'SM P M1 '
165 C CALL DIAGNOSTICS_ADDTOLIST( diagNum,
166 C I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
167
168 C diagName = 'DCO2FLXP'
169 C diagTitle = 'delta CO2 flux due to piston velocity (mol/m2/s)'
170 C diagUnits = 'mol/m2/s '
171 C diagCode = 'SM P M1 '
172 C CALL DIAGNOSTICS_ADDTOLIST( diagNum,
173 C I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
174
175 diagName = 'DCO2FLXC'
176 diagTitle = 'delta CO2 flux due to apCO2 (mol/m2/s)'
177 diagUnits = 'mol/m2/s '
178 diagCode = 'SM P M1 '
179 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
180 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
181
182 diagName = 'DCO2FLXR'
183 diagTitle = 'residual delta CO2 flux due to
184 & bio and physics (mol/m2/s)'
185 diagUnits = 'mol/m2/s '
186 diagCode = 'SM P M1 '
187 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
188 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
189
190 diagName = 'DCO2FLXB'
191 diagTitle = 'delta CO2 flux due to biology (mol/m2/s)'
192 diagUnits = 'mol/m2/s '
193 diagCode = 'SM P M1 '
194 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
195 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
196
197 diagName = 'DCO2FLXP'
198 diagTitle = 'delta CO2 flux due to physics (mol/m2/s)'
199 diagUnits = 'mol/m2/s '
200 diagCode = 'SM P M1 '
201 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
202 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
203
204 diagName = 'BCO2MIXD'
205 diagTitle = 'CO2 flux budget mixing depth (m)'
206 diagUnits = 'm '
207 diagCode = 'SM P M1 '
208 CALL DIAGNOSTICS_ADDTOLIST( diagNum,
209 I diagName, diagCode, diagUnits, diagTitle, 0, myThid )
210
211 C---1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
212
213 #endif /* ALLOW_DIAGNOSTICS */
214
215 RETURN
216 END

  ViewVC Help
Powered by ViewVC 1.1.22