1 | !!--------------------------------------------------------------------- |
---|
2 | !! |
---|
3 | !! COMMON passivetrc.hamocc3.h |
---|
4 | !! ****************************** |
---|
5 | !! |
---|
6 | !! purpose : |
---|
7 | !! --------- |
---|
8 | !! INCLUDE COMMON FILE for HAMO!!3 biological model |
---|
9 | !! |
---|
10 | !! modifications : |
---|
11 | !! ------------- |
---|
12 | !! original : 00-02 (O. Aumont) |
---|
13 | !! |
---|
14 | !!--------------------------------------------------------------------- |
---|
15 | !! opa8, ipsl (11/96) |
---|
16 | !!--------------------------------------------------------------------- |
---|
17 | !! |
---|
18 | #if defined key_trc_hamocc3 |
---|
19 | !! |
---|
20 | !!---------------------------------------------------------------------- |
---|
21 | !! |
---|
22 | !! COMMON/cchem1/ : Variable for chemistry of the CO2 cycle |
---|
23 | !! |
---|
24 | !! --------------------------------------------------------------------- |
---|
25 | !! |
---|
26 | REAL akb3(jpi,jpj,jpk), ak13(jpi,jpj,jpk), ak23(jpi,jpj,jpk) |
---|
27 | REAL aksp(jpi,jpj,jpk), co3(jpi,jpj,jpk), hi(jpi,jpj,jpk) |
---|
28 | REAL borat(jpi,jpj,jpk) |
---|
29 | |
---|
30 | COMMON/cchem1/ akb3,ak13,ak23,aksp,co3,hi,borat |
---|
31 | !! |
---|
32 | !!---------------------------------------------------------------------- |
---|
33 | !! |
---|
34 | !! COMMON/cchem2/ : Variable for chemistry of the CO2 cycle |
---|
35 | !! |
---|
36 | !! --------------------------------------------------------------------- |
---|
37 | !! |
---|
38 | !! dmin3 : vertical distribution of newly formed biogenic poc |
---|
39 | !! disc3 : vertical distribution of newly formed biogenic caco3 |
---|
40 | !! diss3 : vertical distribution of newly formed biogenic sio4 |
---|
41 | |
---|
42 | REAL atcco2, atcco20, atcox,h2co3(jpi,jpj),bobsx(jpi,jpj) |
---|
43 | REAL discl(jpi,jpj),akw3(jpi,jpj,jpk) |
---|
44 | REAL dissl(jpi,jpj),dmin(jpi,jpj),dminl(jpi,jpj),fugaci(jpi,jpj) |
---|
45 | REAL diss3(jpi,jpj,jpk),dmin3(jpi,jpj,jpk),disc3(jpi,jpj,jpk) |
---|
46 | REAL str,strn(jpi,jpj),fbod(jpi,jpj),sedlam,sedpoc(jpi,jpj) |
---|
47 | REAL sedcal(jpi,jpj),sedsil(jpi,jpj) |
---|
48 | REAL sedlostpoc, sedlostsil,sedlostcal |
---|
49 | |
---|
50 | COMMON/cchem2/ atcco2, atcco20, atcox, h2co3, bobsx, & |
---|
51 | &, akw3, dmin, fugaci, str, strn |
---|
52 | COMMON/cotsedim/sedlam,dmin3,dminl,discl,dissl,diss3,disc3,& |
---|
53 | & sedpoc,sedcal,sedsil,fbod,& |
---|
54 | & sedlostpoc, sedlostsil,sedlostcal |
---|
55 | |
---|
56 | !! |
---|
57 | !!---------------------------------------------------------------------- |
---|
58 | !! |
---|
59 | !! COMMON/cchem3/ : Variable for chemistry of the CO2 cycle |
---|
60 | !! |
---|
61 | !! --------------------------------------------------------------------- |
---|
62 | !! |
---|
63 | !! VARIABLE TYPE PURPOSE. |
---|
64 | !! -------- ---- -------- |
---|
65 | !! |
---|
66 | !! *THIRD* REAL 1./3. HALF PRECISION |
---|
67 | !! *TFREE* REAL APPROXIMATE VALUE OF FREEZING POINT OF |
---|
68 | !! SEAWATER |
---|
69 | !! *TBOIL* REAL SETS UPPER LIMIT FOR SEAWATER TEMPERATURE |
---|
70 | !! *SWEET* REAL SETS LOWER LIMIT FOR SALINITY |
---|
71 | !! *BRINE* REAL SETS UPPER LIMIT FOR SALINITY |
---|
72 | !! *TEN* REAL 10. HALF PRECISION |
---|
73 | !! *THOUSI* REAL 0.001 HALF PRECICION |
---|
74 | !! *SMICR* REAL 1E-6 HALF PRECISION |
---|
75 | !! *SALCHL* REAL CONVERSION FACTOR TO CALCULATE CHLORID |
---|
76 | !! CONCENTRATION |
---|
77 | !! S(O/OO)=1.80655*CL(O/OO) |
---|
78 | !! KALLE/DIETRICH , P. 60, OR |
---|
79 | !! WOOSTER ET AL., 1969. |
---|
80 | !! *TEMZER* REAL ZERO DEG C EXPRESSED AS ABS. TEMPERATURE |
---|
81 | !! (DEG KELVIN) |
---|
82 | !! *PERC* REAL 0.01 HALF PRECISION |
---|
83 | !! *AKCC1* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED |
---|
84 | !! APPARENT SOLUBILITY PRODUCT K'SP OF |
---|
85 | !! CALCITE (INGLE, 1800, EQ. 6) |
---|
86 | !! (CF. BROECKER ET AL., 1982) |
---|
87 | !! *AKCC2* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED |
---|
88 | !! APPARENT SOLUBILITY PRODUCT K'SP OF |
---|
89 | !! CALCITE (INGLE, 1800, EQ. 6) |
---|
90 | !! (CF. BROECKER ET AL., 1982) |
---|
91 | !! *AKCC3* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED |
---|
92 | !! APPARENT SOLUBILITY PRODUCT K'SP OF |
---|
93 | !! CALCITE (INGLE, 1800, EQ. 6) |
---|
94 | !! (CF. BROECKER ET AL., 1982) |
---|
95 | !! *AKCC4* REAL COEFFICIENT FOR EMPIRICALLY DETERMINED |
---|
96 | !! APPARENT SOLUBILITY PRODUCT K'SP OF |
---|
97 | !! CALCITE (INGLE, 1800, EQ. 6) |
---|
98 | !! (CF. BROECKER ET AL., 1982) |
---|
99 | !! *ARAFRA* REAL FRACTION OF ARAGONITE IN BIOGENIC CACO3 |
---|
100 | !! PARTICLES (E.G. 0.3 FOR 30 PERCENT) |
---|
101 | !! *CALFRA* REAL FRACTION OF CALCITE IN BIOGENIC CACO3 |
---|
102 | !! PARTICLES (E.G. 0.3 FOR 30 PERCENT) |
---|
103 | !! *ARACAL* REAL FACTOR TO CONVERT APP. CALCITE SOLUBILITY |
---|
104 | !! PRODUCT (0 DBAR) TO THE APP. SOLUBILITY |
---|
105 | !! PRODUCT OF ARAGONITE (BERNER, 1976; |
---|
106 | !! CF. BROECKER ET AL., 1982) |
---|
107 | !! *DEVK1* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
108 | !! OF 1. DISSOCIATION CONSTANT OF CARBONIC |
---|
109 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
110 | !! (CF. BROECKER ET AL., 1982) |
---|
111 | !! *DEVK2* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
112 | !! OF 2. DISSOCIATION CONSTANT OF CARBONIC |
---|
113 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
114 | !! (CF. BROECKER ET AL., 1982) |
---|
115 | !! *DEVKB* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
116 | !! OF 1. DISSOCIATION CONSTANT OF BORIC |
---|
117 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
118 | !! (CF. BROECKER ET AL., 1982) |
---|
119 | !! *DEVK1T* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
120 | !! OF FIRST DISSOCIATION CONSTANT OF CARBONIC |
---|
121 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
122 | !! (CF. BROECKER ET AL., 1982) |
---|
123 | !! *DEVK2T* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
124 | !! OF SECOND DISSOCIATION CONSTANT OF CARBONIC |
---|
125 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
126 | !! (CF. BROECKER ET AL., 1982) |
---|
127 | !! *DEVKBT* REAL COEFFICIENT FOR SEAWATER PRESSURE CORRECTION |
---|
128 | !! OF DISSOCIATION CONSTANT OF BORIC |
---|
129 | !! ACID AFTER CULBERSON AND PYTKOWICZ, 1968 |
---|
130 | !! (CF. BROECKER ET AL., 1982) |
---|
131 | !! *DEVKS* REAL COEFFICIENT FOR PRESSURE CORRECTION OF |
---|
132 | !! SOLUBILITY PRODUCT OF CALCITE OR ARAGONITE |
---|
133 | !! AFTER EDMOND AND GIESKES (1970), P. 1285 |
---|
134 | !! (REFERENCE TO CULBERSON AND PYTKOWICZ, 1968, |
---|
135 | !! AS DONE IN BROECKER ET AL., 1982, IS |
---|
136 | !! NOT CORRECT) |
---|
137 | !! *DEVKST* REAL COEFFICIENT FOR PRESSURE CORRECTION OF |
---|
138 | !! SOLUBILITY PRODUCT OF CALCITE OR ARAGONITE |
---|
139 | !! AFTER EDMOND AND GIESKES (1970), P. 1285, |
---|
140 | !! IN TERM WITH TEMPERATURE |
---|
141 | !! (REFERENCE TO CULBERSON AND PYTKOWICZ, 1968, |
---|
142 | !! AS DONE IN BROECKER ET AL., 1982, IS |
---|
143 | !! NOT CORRECT) |
---|
144 | !! *RGAS* REAL UNIVERSAL GAS CONSTANT (BOLTZMANN'S CONSTANT |
---|
145 | !! TIMES AVOGADRO'S CONSTANT = |
---|
146 | !! 1.3804E-16*6.023*10E+23=8.3143E+7 ERG/GRD*MOL= |
---|
147 | !! 83.143E+6 ERG/GRD*MOL=8.3143 J/K*MOL) |
---|
148 | !! MULTIPLIED WITH 10 (TO A!!OUNT FOR |
---|
149 | !! CHANGE FROM BAR TO DBAR) |
---|
150 | !! (CF. EDMOND AND GIESKES, 1970. P. 1285, |
---|
151 | !! BROECKER ET AL., 1982, P. 79) |
---|
152 | !! *BOR1* REAL TOTAL BORON CONTENT IN G/KG AT CL=19 O/OO |
---|
153 | !! (S=35) |
---|
154 | !! (CF. RILEY AND SKIRROW, VOL. 1, P. 648) |
---|
155 | !! *BOR2* REAL INVERSE OF ATOMIC WEIGHT OF BORON FOR |
---|
156 | !! CONVERTING SPECIFIC TOTAL BORATE IN |
---|
157 | !! CONCENTRATIONS |
---|
158 | !! *OXYCO* REAL INVERS OF NORMAL MOLAL VOLUME OF AN |
---|
159 | !! IDEAL GAS [CM**-3] |
---|
160 | !! *C00* REAL VOLUMETRIC SOLUBILITY CONSTANT A1 FOR |
---|
161 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
162 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
163 | !! *C01* REAL VOLUMETRIC SOLUBILITY CONSTANT A2 FOR |
---|
164 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
165 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
166 | !! *C02* REAL VOLUMETRIC SOLUBILITY CONSTANT A3 FOR |
---|
167 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
168 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
169 | !! *C03* REAL VOLUMETRIC SOLUBILITY CONSTANT B1 FOR |
---|
170 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
171 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
172 | !! *C04* REAL VOLUMETRIC SOLUBILITY CONSTANT B2 FOR |
---|
173 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
174 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
175 | !! *C05* REAL VOLUMETRIC SOLUBILITY CONSTANT B3 FOR |
---|
176 | !! THE SOLUBILITY OF CO2 IN ML/L FROM AIR |
---|
177 | !! AT ONE ATMOSPHERE (WEISS, 1974) |
---|
178 | !! *C10* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER |
---|
179 | !! EDMOND AND GIESKES (1970) |
---|
180 | !! *C11* REAL COEFF. FOR F1. H2CO3 DISSOC. CONST. AFTER |
---|
181 | !! EDMOND AND GIESKES (1970) |
---|
182 | !! *C12* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER |
---|
183 | !! EDMOND AND GIESKES (1970) |
---|
184 | !! *C13* REAL COEFF. FOR 1. H2CO3 DISSOC. CONST. AFTER |
---|
185 | !! EDMOND AND GIESKES (1970) |
---|
186 | !! *C20* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER |
---|
187 | !! EDMOND AND GIESKES (1970) |
---|
188 | !! *C21* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER |
---|
189 | !! EDMOND AND GIESKES (1970) |
---|
190 | !! *C22* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER |
---|
191 | !! EDMOND AND GIESKES (1970) |
---|
192 | !! *C23* REAL COEFF. FOR 2. H2CO3 DISSOC. CONST. AFTER |
---|
193 | !! EDMOND AND GIESKES (1970) |
---|
194 | !! *CB0* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER |
---|
195 | !! EDMOND AND GIESKES (1970) |
---|
196 | !! *CB1* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER |
---|
197 | !! EDMOND AND GIESKES (1970) |
---|
198 | !! *CB2* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER |
---|
199 | !! EDMOND AND GIESKES (1970) |
---|
200 | !! *CB3* REAL COEFF. FOR 1. H3BO3 DISSOC. CONST. AFTER |
---|
201 | !! EDMOND AND GIESKES (1970) |
---|
202 | !! *CW0* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) |
---|
203 | !! *CW1* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) |
---|
204 | !! *CW2* REAL COEFF. FOR KW (DICKSON AND RILEY, 1979) |
---|
205 | !! (CORRECTED A!!ORDING TO B. BACASTOW, |
---|
206 | !! PERS. COMMUN., 1988) |
---|
207 | !! *OX0* REAL VOLUMETRIC SOLUBILITY CONSTANT A1 FOR |
---|
208 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
209 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
210 | !! *OX1* REAL VOLUMETRIC SOLUBILITY CONSTANT A2 FOR |
---|
211 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
212 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
213 | !! *OX2* REAL VOLUMETRIC SOLUBILITY CONSTANT A3 FOR |
---|
214 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
215 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
216 | !! *OX3* REAL VOLUMETRIC SOLUBILITY CONSTANT A4 FOR |
---|
217 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
218 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
219 | !! *OX4* REAL VOLUMETRIC SOLUBILITY CONSTANT B1 FOR |
---|
220 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
221 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
222 | !! *OX5* REAL VOLUMETRIC SOLUBILITY CONSTANT B2 FOR |
---|
223 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
224 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
225 | !! *OX6* REAL VOLUMETRIC SOLUBILITY CONSTANT B3 FOR |
---|
226 | !! THE SOLUBILITY OF OXYGEN IN ML/L FROM |
---|
227 | !! MOIST AIR AT ONE ATMOSPHERE (WEISS, 1970) |
---|
228 | !! *T* REAL DUMMY VARIABLE, ABSOLUTE SEAWATER TEMP. |
---|
229 | !! *qtt* REAL DUMMY VARIABLE, ABSOLUTE SEAWATER TEMP., |
---|
230 | !! DIVIDED BY 100. |
---|
231 | !! *S* REAL DUMMY VARIABLE, SALINITY |
---|
232 | !! *CL* REAL CHLORINITY (CL(O/OO)=S(O/OO)/1.80655) |
---|
233 | !! AFTER WOOSTER ET AL., 1969 |
---|
234 | !! (C.F. KALLE/DIETRICH , P. 60) |
---|
235 | !! *CEK0* REAL LN(K0), LOGARITHM OF CO2 SOLUBILITY IN |
---|
236 | !! SEAWATER IN VOLUMETRIC UNITS (EQ. 12 IN |
---|
237 | !! WEISS, 1974) |
---|
238 | !! *CK1* REAL PK1-VALUE (K1= 1. H2CO3 DISSOC. CONST..), |
---|
239 | !! AFTER EDMOND AND GIESKES (1970) |
---|
240 | !! *CK2* REAL PK2-VALUE (K2= 2. H2CO3 DISSOC. CONST..), |
---|
241 | !! AFTER EDMOND AND GIESKES (1970) |
---|
242 | !! *CKB* REAL PKB-VALUE (KB= 1. H3BO3 DISSOC. CONST..) |
---|
243 | !! AFTER EDMOND AND GIESKES (1970) |
---|
244 | !! *CKW* REAL PKW-VALUE (KW=H2O DISSOC. CONST.) AFTER |
---|
245 | !! DICKSON AND RILEY (1979) |
---|
246 | !! *OXY* REAL LN(C*), LOGARITHM OF O2 SOLUBILITY IN |
---|
247 | !! SEAWATER IN VOLUMETRIC UNITS (EQ. 4 IN |
---|
248 | !! WEISS, 1970) |
---|
249 | !! *AK1* REAL K1, 1. H2CO3 DISSOC. CONSTANT |
---|
250 | !! (EDMOND AND GIESKES, 1970) |
---|
251 | !! *AK2* REAL K2, 2. H2CO3 DISSOC. CONSTANT |
---|
252 | !! (EDMOND AND GIESKES, 1970) |
---|
253 | !! *AKB* REAL KB, 1. H3BO3 DISSOC. CONSTANT |
---|
254 | !! (EDMOND AND GIESKES, 1970) |
---|
255 | !! *AKW* REAL KW, H2O DISSOC. CONSTANT, LIT ? |
---|
256 | !! *AK0* REAL EXP(LN(K0))=K0 CO2 SOLUBILITY IN SEAWATER |
---|
257 | !! IN VOLUMETRIC UNITS (ML/L)(WEISS, 1974, |
---|
258 | !! CF. EQ. 12) |
---|
259 | !! *RRR* REAL SIGMA-T IN OCEAN MODEL, DUMMY VARIABLE |
---|
260 | !! (USED FOR CALCULATION OF TOTAL BORAT |
---|
261 | !! CONCENTRATION) (SIGMA-T=RHO(S,T,0)/1000.) |
---|
262 | !! *BOR* REAL TOTAL BORAT CONCENTRATION , DUMMY VAR. |
---|
263 | !! *TENM7* REAL 10.**-7 HALF PRECISION |
---|
264 | !! *P* REAL APPROXIMATE PRESSURE AT DEPTH OF U-POINTS |
---|
265 | !! IN BAR, DUMMY VARIABLE |
---|
266 | !! *AKSP0* REAL CACO3 SOLUBILITY PRODUCT AT P=0 DBAR |
---|
267 | !! A!!ORDING TO INGLE (1800), EQ. 6; THE |
---|
268 | !! CITATION OF CULBERSON AND PYTKOWICZ, 1968, |
---|
269 | !! IN BROECKER ET AL., 1982, IS PRESUMABLY |
---|
270 | !! NOT CORRECT) |
---|
271 | !! *AKSP(jpi,jpj,jpk)* REAL CACO3 SOLUBILITY PRODUCT AT IN SITU PRESSURE |
---|
272 | !! FOLLOWING THE PROCEDURE DESCRIBED IN EDMOND |
---|
273 | !! AND GIESKES (1970), P. 1285 |
---|
274 | !! *CP* REAL TERM IN EXPONENT OF EQUATIONS FOR PRESSURE |
---|
275 | !! CORRECTION OF DISSOC. CONSTANTS (CARB., |
---|
276 | !! BOR. ACID) AND CALCITE SOLUB. PRODUCT |
---|
277 | !! (CF. BROECKER ET AL., 1982, EDMOND AND |
---|
278 | !! GIESKES, 1970) |
---|
279 | !! *TC* REAL TEMPERATURE AT OCEAN GRID POINTS (DEG C), |
---|
280 | !! DUMMY VARIABLE |
---|
281 | !! *KI* INTEGER COUNTS ITERATIONS FOR NEWTON-RAPHSON METHOD |
---|
282 | !! FOR INITIATION OF [CO3--] AND [H+] |
---|
283 | !! *H* REAL [H+], DUMMY VARIABLE |
---|
284 | !! *R* REAL [CO3--] [MOLE/L], DUMMY VARIABLE |
---|
285 | !! *ALKA* REAL GIVEN ALKALINITY [EQV/L], DUMMY VARIABLE |
---|
286 | !! *C* REAL GIVEN [SUM(12C)O2] [MOLE/L], DUMMY VARIABLE |
---|
287 | !! *A* REAL ALKALINITY [EQV/L] AS FUNCTION OF [CO3--] |
---|
288 | !! AND [H+], DUMMY VARIABLE |
---|
289 | !! *DCDS* REAL LOCAL DERIVATIVE |
---|
290 | !! [SUM(CO2)]([CO3--],H+]) -> [CO3--] |
---|
291 | !! *DADS* REAL LOCAL DERIVATIVE |
---|
292 | !! [ALK]([CO3--],H+]) -> [CO3--] |
---|
293 | !! *DCDH* REAL LOCAL DERIVATIVE |
---|
294 | !! [SUM(CO2)]([CO3--],H+]) -> [H+] |
---|
295 | !! *DADH* REAL LOCAL DERIVATIVE |
---|
296 | !! [ALK]([CO3--],H+]) -> [H+] |
---|
297 | !! *F* REAL FUNCTION [ALK] GIVEN MINUS [ALK] IN TERMS |
---|
298 | !! OF [CO3--] AND [H+] |
---|
299 | !! *!!* REAL [SUM(CO2)] [MOLE/L] AS FUNCTION OF |
---|
300 | !! [CO3--] AND [H+], DUMMY VARIABLE |
---|
301 | !! *GG* REAL FUNCTION [SUM(CO2)] GIVEN MINUS [SUM(CO2)] |
---|
302 | !! IN TERMS OF [CO3--] AND [H+] |
---|
303 | !! *DETI* REAL DETERMINANT WITH LOCAL DERIVATIVES FOR |
---|
304 | !! NEWTON-RAPHSON ITERATION |
---|
305 | !! |
---|
306 | REAL brine, akcc1, akcc2, akcc3, akcc4 |
---|
307 | REAL arafra, calfra, aracal, devk1, devk2, devkb |
---|
308 | REAL devk1t, devk2t, devkbt, devkst, devks |
---|
309 | REAL bor1, bor2, c00, c01, c02, c03, c04, c05, c10, c11 |
---|
310 | REAL c12, c13, c20, c21, c22, c23, cb0, cb1, cb2, cb3 |
---|
311 | REAL c14, c15, c16, c17, c24, c25, c26, c27 |
---|
312 | REAL cb4, cb5, cb6, cb7, cb8, cb9, cb10, cb11 |
---|
313 | REAL cw3, cw4, cw5, cw6 |
---|
314 | REAL cw0, cw1, cw2, ox0, ox1, ox2, ox3, ox4 |
---|
315 | REAL cek0, ckb, ck1, ck2, ckw, ak1, ak2, ak0 |
---|
316 | REAL bor, aksp0, eak0, co2, ah, ba |
---|
317 | REAL tenm7, smicr, thousi, perc, third, ten, tfree, tboil |
---|
318 | REAL sweet, salchl, temzer, sucall, rgas, oxyco, ox5 |
---|
319 | REAL vi0, vi1, vi2, vi3, vi4, vi5, oxy, vit, rrr |
---|
320 | REAL soloxy, tco2, qtt |
---|
321 | REAL chemc(jpi,jpj,3), cpexp |
---|
322 | |
---|
323 | COMMON/cchem3/ brine, akcc1, akcc2, akcc3, akcc4 & |
---|
324 | &, arafra, calfra, aracal, devk1, devk2, devkb & |
---|
325 | &, devk1t, devk2t, devkbt, devkst, devks & |
---|
326 | &, bor1, bor2, c00, c01, c02, c03, c04, c05, c10, c11 & |
---|
327 | &, c12, c13, c20, c21, c22, c23, cb0, cb1, cb2, cb3 & |
---|
328 | &, c14, c15, c16, c17, c24, c25, c26, c27 & |
---|
329 | &, cb4, cb5, cb6, cb7, cb8, cb9, cb10, cb11 & |
---|
330 | &, cw3, cw4, cw5, cw6 & |
---|
331 | &, cw0, cw1, cw2, ox0, ox1, ox2, ox3, ox4 & |
---|
332 | &, cek0, ckb, ck1, ck2, ckw, ak1, ak2, ak0 & |
---|
333 | &, bor, aksp0, eak0, co2, ah, ba & |
---|
334 | &, tenm7, smicr, thousi, perc, third, ten, tfree, tboil & |
---|
335 | &, sweet, salchl, temzer, sucall, rgas, oxyco, ox5 & |
---|
336 | &, vi0, vi1, vi2, vi3, vi4, vi5, oxy, vit, rrr & |
---|
337 | &, soloxy, tco2, qtt, chemc, cpexp |
---|
338 | !! |
---|
339 | !!---------------------------------------------------------------------- |
---|
340 | !! |
---|
341 | !! COMMON/cotcon/ : Time variables |
---|
342 | !! |
---|
343 | !! --------------------------------------------------------------------- |
---|
344 | !! |
---|
345 | INTEGER nyrbeg, nyrend, nyr, nyrfend, nyrf |
---|
346 | PARAMETER (nyrbeg=1950,nyrend=1990) |
---|
347 | PARAMETER (nyr=nyrend+1-nyrbeg) |
---|
348 | PARAMETER (nyrfend=2300,nyrf=nyrfend+1-nyrbeg) |
---|
349 | |
---|
350 | INTEGER ntspmn, iabsyr |
---|
351 | REAL ttsec, ttyear, tspyr, absyr |
---|
352 | REAL rfact, rfactr |
---|
353 | |
---|
354 | INTEGER ilecvit |
---|
355 | REAL xtvit |
---|
356 | |
---|
357 | COMMON/cottim/ ttsec, ttyear, tspyr, ntspmn, absyr, iabsyr & |
---|
358 | &, rfact, rfactr, xtvit, ilecvit |
---|
359 | |
---|
360 | !! |
---|
361 | !!---------------------------------------------------------------------- |
---|
362 | !! |
---|
363 | !! COMMON/cotgas/ : Gas exchange |
---|
364 | !! |
---|
365 | !! --------------------------------------------------------------------- |
---|
366 | !! |
---|
367 | INTEGER igaswind,icice, nvit1p,nutrstm |
---|
368 | REAL gasfac,qcumul(jptra),tdepart(jptra),tcumul(jptra) |
---|
369 | REAL kgwanin(jpi,jpj), wsmo(jpi,jpj,12), kgwanmo(jpi,jpj,12) |
---|
370 | REAL cicemo(jpi,jpj,12) |
---|
371 | REAL patm(jpi,jpj),ppres(jpi,jpj,12) |
---|
372 | |
---|
373 | COMMON/cotgas/ gasfac, igaswind, nvit1p, icice & |
---|
374 | &,kgwanin, wsmo, kgwanmo, cicemo & |
---|
375 | &,patm,ppres,qcumul,tdepart,tcumul |
---|
376 | |
---|
377 | !! |
---|
378 | !!---------------------------------------------------------------------- |
---|
379 | !! |
---|
380 | !! COMMON/cotcon/ : Variable for testing convergence |
---|
381 | !! |
---|
382 | !! --------------------------------------------------------------------- |
---|
383 | !! |
---|
384 | REAL rconv(jpj), rconvs |
---|
385 | REAL rsumds, rsumdx |
---|
386 | |
---|
387 | COMMON/cotcon/ rconv, rconvs, rsumds, rsumdx |
---|
388 | !! |
---|
389 | !!---------------------------------------------------------------------- |
---|
390 | !! |
---|
391 | !! COMMON/cotham/ : biological parameters specific to HAMOCC3 |
---|
392 | !! |
---|
393 | !! --------------------------------------------------------------------- |
---|
394 | !! |
---|
395 | !! calcon = MEAN TOTAL [CA++] IN SEAWATER (MOLES/KG) |
---|
396 | !! (SEE BROECKER A. PENG, 1982, P. 26) |
---|
397 | !! ([CA++](MOLES/KG)=1.026E-2*(S/35.) AFTER |
---|
398 | !! CULKIN(1965), CF. BROECKER ET AL. 1982) |
---|
399 | !! spocri = variable for sediments |
---|
400 | !! conc0 = MICHAELIS MENTON CONSTANT OF PRODUCTIVITY (HALF SATU- |
---|
401 | !! RATION CONSTANT; NUTRIENT CONCENTRATION, WHERE NUTRIENT |
---|
402 | !! UPTAKE VELOCITY IS HALF ITS POTENTIAL MAXIMUM VALUE; |
---|
403 | !! E.G. DUGDALE (1967), PARSONS AND TAKAHASHI (1973)) |
---|
404 | !! rno3 = REDFIELD RATIO MOLES N (+ HPO4-) / MOLES C |
---|
405 | !! (FOR CHANGE IN ALKALINITY DUE TO PRODUCTION/REMINE- |
---|
406 | !! RALIZATION OF ORGANIC MATTER) |
---|
407 | !! N:C=16:122, SEE TAKAHASHI ET AL. (1985), P. 6923 |
---|
408 | !! |
---|
409 | !! COMPOSITION OF PLANKTONIC MATERIAL AND RESPIRATION: |
---|
410 | !! C106H263O110N16P + 138 O2 -> |
---|
411 | !! -> 106 CO2 + 16 NO3- + HPO4(2-) + 122 H2O + 18 H+ , |
---|
412 | !! SEE DEGENS ET AL. (1984), P. 152 |
---|
413 | !! *PO4R* REAL RATIO (MOLES P)/(MOLES C) (REDFIELD RATIO P:C) |
---|
414 | !! caco3r = ratio c-c(organic)/c-cacao3 in biogenic particles |
---|
415 | !! ----------------------------------------------------- |
---|
416 | !! achtung caco3r = 0.3 hat gut funktioniert bis auf zu wenig eqpco2 |
---|
417 | !! in expofa 0.5 auf 0.2 gesetzt |
---|
418 | !! sio2r = ratio moles c / moles sio4 in biogenic particles |
---|
419 | !! cristl = CRISTALLIZATION RATE IN CACO3 SUPERSATURATED WATER |
---|
420 | !! fecpef = FRACTION OF NEWLY FORMED BIOGENIC POC THAT FALLS |
---|
421 | !! IMMEDIATELY TO THE BOTTOM LAYER |
---|
422 | !! tobox = FRACTION OF NEWLY FORMED BIOGENIC CACO3 PARTICLES |
---|
423 | !! THAT FALLS IMMEDIATELY TO THE BOTTOM LAYER |
---|
424 | !! dispo0 = inertia concentration of caco3 for spontaneous dissolution [moles/l] |
---|
425 | !! oxymin = minimum [o2] [moles/l] for bacterial decomposition organic matter |
---|
426 | !! |
---|
427 | REAL caco3r, rno3, o2ut, sio2r, tobox ,po4r |
---|
428 | REAL prlat(jpi,jpj) |
---|
429 | REAL prcaca(jpi,jpj,jpk), silpro(jpi,jpj,jpk), prorca(jpi,jpj,jpk) |
---|
430 | REAL sco2, dispo0, conc0, silic0 |
---|
431 | REAL cristl, calcon, plafr13, pdb, c13at |
---|
432 | REAL grosip, oxymin, fecpef, bigshl, spocri |
---|
433 | REAL one, half, zero, two, four, fourth |
---|
434 | REAL pendec,pendes,fluexp |
---|
435 | REAL cmask(jpi,jpj),areacot |
---|
436 | INTEGER jpkb |
---|
437 | |
---|
438 | COMMON/cotham/ caco3r, rno3, o2ut, sio2r & |
---|
439 | &, tobox, po4r, prlat & |
---|
440 | &, prcaca, silpro, prorca & |
---|
441 | &, sco2, dispo0, conc0, silic0 & |
---|
442 | &, cristl, calcon, plafr13, pdb, c13at & |
---|
443 | &, grosip, oxymin, fecpef, bigshl, spocri & |
---|
444 | &, one, half, zero, two, four, fourth & |
---|
445 | &, pendec,pendes,fluexp,jpkb & |
---|
446 | & ,cmask, areacot |
---|
447 | |
---|
448 | #endif |
---|
449 | |
---|