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