1 | |
---|
2 | CCC $Header$ |
---|
3 | CCC TOP 1.0 , LOCEAN-IPSL (2005) |
---|
4 | C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
5 | C --------------------------------------------------------------------------- |
---|
6 | C $Id$ |
---|
7 | CCC----------------------------------------------------------------- |
---|
8 | CCC |
---|
9 | CCC ROUTINE trcini.hamocc3.h |
---|
10 | CCC ************************ |
---|
11 | CCC |
---|
12 | CCC PURPOSE : |
---|
13 | CCC --------- |
---|
14 | CCC Initialisation of HAMOCC3 biological and chemical variables |
---|
15 | CCC |
---|
16 | CC METHOD : |
---|
17 | CC ------- |
---|
18 | CC 1) SET CONSTANTS FOR CARBONATE CHEMISTRY AS DESCRIBED IN |
---|
19 | CC IN BROECKER ET AL. (1982, GEOSECS) AND EDMOND A. GIESKES |
---|
20 | CC (1970) |
---|
21 | CC 2) INITIATE [CO3--] AND PH-VALUE BY ITERATION |
---|
22 | CC (NEWTON-RAPHSON METHOD FOR SOLVING NONLINEAR SIMULTANEOUS |
---|
23 | CC EQUATIONS, SEE E.G. SCARBOROUGH, J. (1958)) |
---|
24 | CC |
---|
25 | CC This sub-routine merges previous initia2.F bioini.F bioini2.F |
---|
26 | CC |
---|
27 | CC INPUT : |
---|
28 | CC ----- |
---|
29 | CC common |
---|
30 | CC all the common defined in opa |
---|
31 | CC |
---|
32 | CC |
---|
33 | CC OUTPUT : : no |
---|
34 | CC ------ |
---|
35 | CC |
---|
36 | CC WORKSPACE : |
---|
37 | CC --------- |
---|
38 | CC |
---|
39 | CC EXTERNAL : |
---|
40 | CC -------- |
---|
41 | CC SEAICE |
---|
42 | CC WSJCD |
---|
43 | CC WSEK |
---|
44 | CC RHO |
---|
45 | CC |
---|
46 | CC MODIFICATIONS: |
---|
47 | CC -------------- |
---|
48 | CC original : 1988-07 E. MAIER-REIMER MPI HAMBURG |
---|
49 | CC additions : 1999-10 O. Aumont and C. Le Quere |
---|
50 | CC additions : 2001-03 0. Aumont and E. Kastenare : Changes in computation |
---|
51 | CC of the export profiles for silicate and calcite |
---|
52 | CC |
---|
53 | CC REFERENCE for biology: |
---|
54 | CC ---------------------- |
---|
55 | CC |
---|
56 | CC DEGENS, E.T, S. KEMPE, AND A. SPITZY (1984) |
---|
57 | CC CARBON DIOXIDE: A BIOGEOCHEMICAL PORTRAIT. |
---|
58 | CC IN: THE HANDBOOK OF ENVIRONMENTAL CHEMISTRY, VOLUME 1/ |
---|
59 | CC PART C, O. HUTZINGER, ED., SPRINGER-VERLAG, BERLIN, |
---|
60 | CC HEIDELBERG, PP. 127-215. |
---|
61 | CC |
---|
62 | CC DUGDALE. R.C. (1967) |
---|
63 | CC NUTRIENT LIMITATION IN THE SEA: DYNAMICS, IDENTIFICATION |
---|
64 | CC AND SIGNIFICANCE. |
---|
65 | CC LIMNOLOGY AND OCEANOGRAPHY, VOL.12, 685-695. |
---|
66 | CC |
---|
67 | CC PARSONS, T.R., AND M. TAKAHASHI (1973) |
---|
68 | CC BIOLOGICAL OCEANOGRAPHIC PROCESSES. |
---|
69 | CC PERGAMON PRESS, 186 PP. |
---|
70 | CC |
---|
71 | CC TAKAHASHI, T., W.S. BROECKER, AND S. LANGER (1985) |
---|
72 | CC REDFIELD RATIO BASED ON CHEMICAL DATA FROM ISOPYCNAL |
---|
73 | CC SURFACES. |
---|
74 | CC JOURNAL OF GEOPHYSICAL RESEARCH, 90(C4), 6907-6924. |
---|
75 | CC |
---|
76 | CC REFERENCE for chemistry: |
---|
77 | CC ----------------------- |
---|
78 | CC |
---|
79 | CC BERNER, R. A. (1976) |
---|
80 | CC THE SOLUBILITY OF CALCITE AND ARAGONITE IN SEA WATER |
---|
81 | CC AT ATMOSPHERIC PRESSURE AND 34.5 O/OO SALINITY. |
---|
82 | CC AMERICAN JOURNAL OF SCIENCE, VOL. 276, 713-730. |
---|
83 | CC (K'SP(ARAGONITE)=1.45 K'SP(CALCITE)) |
---|
84 | CC |
---|
85 | CC BROECKER, W.S., D.W. SPENCER, AND H. CRAIG (1982) |
---|
86 | CC GEOSECS PACIFIC EXPEDITION. VOL. 3.. HYDROGRAPHIC DATA |
---|
87 | CC 1973-1974, SUPERINTENDANT OF DOCUMENTS, U.S. GOVERNMENT |
---|
88 | CC PRINTING OFFICE, WASHINGTON, D.C., 137 PP.. |
---|
89 | CC |
---|
90 | CC CULBERSON, C.H., AND R.M. PYTKOWICZ (1968) |
---|
91 | CC EFFECT ON PRESSURE ON CARBONIC ACID, BORIC ACID AND THE PH |
---|
92 | CC IN SEA WATER. |
---|
93 | CC LIMNOLOGY AND OCEANOGRAPHY, VOL. 13, 403-417. |
---|
94 | CC |
---|
95 | CC DICKSON, A.G., AND J.P. RILEY (1979) |
---|
96 | CC THE ESTIMATION OF ACID DISSOCIATION CONSTANTS IN SEAWATER |
---|
97 | CC MEDIA FROM POTENTIOMETRIC TITRATIONS WITH STRONG BASE. |
---|
98 | CC I. THE IONIC PRODUCT OF WATER - KW. |
---|
99 | CC MARINE CHEMISTRY, VOL. 7, 89-99. |
---|
100 | CC |
---|
101 | CC EDMOND, J.M., AND J.M.T.M. GIESKES (1970) |
---|
102 | CC ON THE CALCULATION OF THE DEGREE OF SATURATION OF SEA WATER |
---|
103 | CC WITH RESPECT TO CALCIUM CARBONATE UNDER IN SITU CONDITIONS. |
---|
104 | CC GEOCHIM. ET COSMOCHIM. ACTA, 34, 1261-1291. |
---|
105 | CC |
---|
106 | CC INGLE, S.E. (1800) |
---|
107 | CC SOLUBILITY OF CALCITE IN THE OCEAN. |
---|
108 | CC MARINE CHEMISTRY, VOL. 3, 301-319. |
---|
109 | CC |
---|
110 | CC INGLE, S.E., C.H. CULBERSON, J.E. HAWLEY, AND R.M. PYTKOWICZ |
---|
111 | CC (1973) THE SOLUBILITY OF CALCITE IN SEAWATER AT ATMOSPHERIC |
---|
112 | CC PRESSURE AND 35 O/OO SALINITY. |
---|
113 | CC MARINE CHEMISTRY, VOL. 1, 295-307. |
---|
114 | CC |
---|
115 | CC RILEY, J. P., AND G. SKIRROW, EDS. (1965) |
---|
116 | CC CHEMICAL OCEANOGRAPHY. VOL. 1, 712 PP., ACADEMIC PRESS, |
---|
117 | CC LONDON A. NEW YORK. |
---|
118 | CC |
---|
119 | CC SCARBOROUGH, J. (1958) NUMERICAL MATHEMATICAL ANALYSIS. |
---|
120 | CC OXFORD UNIVERSITY PRESS, LONDON, 4TH ED., 576 PP.. |
---|
121 | CC |
---|
122 | CC WEISS, R. F. (1970) THE SOLUBILITY OF NITROGEN |
---|
123 | CC OXYGEN AND ARGON IN WATER AND SEAWATER. |
---|
124 | CC DEEP-SEA RESEARCH, VOL. 17, 721-735. |
---|
125 | CC |
---|
126 | CC WEISS, R. F. (1974) |
---|
127 | CC CARBON DIOXIDE IN WATER AND SEAWATER: THE SOLUBILITY OF A |
---|
128 | CC NON IDEAL GAS. MARINE CHEMISTRY, VOL. 2, 203-215. |
---|
129 | CC |
---|
130 | CC WOOSTER, W.S., A.J. LEE, AND G. DIETRICH (1969) |
---|
131 | CC REDEFINITION OF SALINITY. Z. GEOPHYS., VOL.35, 611-613. |
---|
132 | CC |
---|
133 | CC BROECKER, W.S., D.W. SPENCER, AND H. CRAIG (1982) |
---|
134 | CC GEOSECS PACIFIC EXPEDITION. VOL. 3.. HYDROGRAPHIC DATA |
---|
135 | CC 1973-1974, SUPERINTENDANT OF DOCUMENTS, U.S. GOVERNMENT |
---|
136 | CC PRINTING OFFICE, WASHINGTON, D.C., 137 PP.. |
---|
137 | CC |
---|
138 | CCC--------------------------------------------------------------------- |
---|
139 | CC local declarations |
---|
140 | CC ================== |
---|
141 | INTEGER mo,ke1 |
---|
142 | |
---|
143 | REAL zdm0(jpi,jpj,jpk),zds0(jpi,jpj,jpk) |
---|
144 | REAL zdc0(jpi,jpj,jpk),zrrc(jpi,jpj) |
---|
145 | REAL zrro(jpi,jpj), zrrs(jpi,jpj), tobbb, tob, tun |
---|
146 | REAL zfluo, zfluu |
---|
147 | REAL ztest |
---|
148 | C |
---|
149 | C |
---|
150 | C |
---|
151 | C |
---|
152 | C |
---|
153 | C 1. initialization |
---|
154 | C ----------------- |
---|
155 | C |
---|
156 | C initialisation of local variables |
---|
157 | C |
---|
158 | rfact = rdttra(1) * float(ndttrc) |
---|
159 | rfactr = 1./rfact |
---|
160 | WRITE(numout,*) ' Tracer time step=',rfact,' rdt=',rdt |
---|
161 | C |
---|
162 | CC |
---|
163 | CC Initialize ice cover |
---|
164 | CC |
---|
165 | CC---------------------------------------------------------------------- |
---|
166 | C |
---|
167 | IF (icice .eq. 0) THEN |
---|
168 | WRITE (numout,*) ' No Sea Ice (default value) ' |
---|
169 | cicemo = 0.0 |
---|
170 | ELSEIF (icice .eq. 1) THEN |
---|
171 | WRITE (numout,*) ' Sea Ice from Reynolds surface data' |
---|
172 | ELSEIF (icice .eq. 2) THEN |
---|
173 | WRITE (numout,*) ' Sea Ice coverage from Walsh and Zwally' |
---|
174 | CALL h3cice |
---|
175 | ELSEIF (icice .eq. 3) THEN |
---|
176 | WRITE (numout,*) ' Sea Ice Coverage from coupled run' |
---|
177 | ENDIF |
---|
178 | c |
---|
179 | CC---------------------------------------------------------------------- |
---|
180 | CC |
---|
181 | CC Initialize gas exchange |
---|
182 | CC |
---|
183 | CC---------------------------------------------------------------------- |
---|
184 | C |
---|
185 | C |
---|
186 | IF (igaswind .eq. 0) THEN |
---|
187 | WRITE (numout,*) ' No winds read' |
---|
188 | C |
---|
189 | ELSEIF (igaswind .eq. 1) THEN |
---|
190 | C |
---|
191 | C Hellerman and Rosenstein according to Jean-Claude Dutay |
---|
192 | C |
---|
193 | WRITE (numout,*)' Interannual wind speed' |
---|
194 | C |
---|
195 | ELSEIF (igaswind .eq. 2) THEN |
---|
196 | C |
---|
197 | C Esbensen and Kushnir according to Bonnie Samuels |
---|
198 | C |
---|
199 | WRITE (numout,*)' Wind speed from Esbensen & Kushnir: ' |
---|
200 | CALL h3cwin |
---|
201 | DO mo=1,12 |
---|
202 | #if defined key_vpp |
---|
203 | CALL READ2S(900,kgwanmo(1,1,mo),jpi,jpj) |
---|
204 | #else |
---|
205 | READ (900) kgwanmo(:,:,mo) |
---|
206 | #endif |
---|
207 | ENDDO |
---|
208 | C |
---|
209 | ELSEIF (igaswind .eq. 3) THEN |
---|
210 | C |
---|
211 | WRITE (numout,*) 'Wind speed from coupled run' |
---|
212 | C |
---|
213 | C In case of offline simulations forced by coupled run output |
---|
214 | C |
---|
215 | ELSE |
---|
216 | C |
---|
217 | WRITE (numout,*) ' Variable igaswind invalid!' |
---|
218 | WRITE (numout,*) ' Rerun model with igaswind = 0, 1, 2, 3' |
---|
219 | STOP |
---|
220 | C |
---|
221 | ENDIF |
---|
222 | C |
---|
223 | C |
---|
224 | C |
---|
225 | CC---------------------------------------------------------------------- |
---|
226 | CC |
---|
227 | CC Initialize biological variables |
---|
228 | CC |
---|
229 | CC---------------------------------------------------------------------- |
---|
230 | C |
---|
231 | spocri = 0.003 |
---|
232 | cristl = 0. |
---|
233 | tobox = 50. |
---|
234 | prorca = 0. |
---|
235 | prcaca = 0. |
---|
236 | silpro = 0. |
---|
237 | #if defined key_trc_p3zd |
---|
238 | sinking = 0. |
---|
239 | nu = 0. |
---|
240 | tortz = 0. |
---|
241 | grazp = 0. |
---|
242 | grazpoc = 0. |
---|
243 | tortp = 0. |
---|
244 | respp = 0. |
---|
245 | respz = 0. |
---|
246 | dipn = 55. |
---|
247 | #endif |
---|
248 | C |
---|
249 | C Set biological ratios |
---|
250 | C --------------------- |
---|
251 | C |
---|
252 | rno3 = (16.+2.)/122. |
---|
253 | po4r = 1./122. |
---|
254 | o2ut = 172./122. |
---|
255 | sio2r = 0.75 |
---|
256 | C |
---|
257 | C Set fractionation factors of 13C, 14C |
---|
258 | C ------------------------------------- |
---|
259 | C |
---|
260 | pdb = 0.011112 |
---|
261 | plafr13 = 0.980 |
---|
262 | C |
---|
263 | C Set e-folding depth for penetration of newly formed biogenic poc and |
---|
264 | c caco3 [cm] |
---|
265 | C ---------------------------------------------------------------------- |
---|
266 | c --------- |
---|
267 | C |
---|
268 | C |
---|
269 | #if defined key_off_degrad |
---|
270 | # if defined key_vpp |
---|
271 | CALL READ3S(903,dmin3,jpi,jpj,jpk) |
---|
272 | CALL READ3S(903,diss3,jpi,jpj,jpk) |
---|
273 | CALL READ3S(903,disc3,jpi,jpj,jpk) |
---|
274 | # else |
---|
275 | READ(903) dmin3 |
---|
276 | READ(903) diss3 |
---|
277 | READ(903) disc3 |
---|
278 | # endif |
---|
279 | WRITE (numout,*) 'OK for read of deg.min' |
---|
280 | #else |
---|
281 | C |
---|
282 | C Calculate vertical distribution of newly formed biogenic poc |
---|
283 | C and caco3 particles in the water column in the case of |
---|
284 | C max. possible bottom depth |
---|
285 | C ------------------------------------------------------------ |
---|
286 | C |
---|
287 | C |
---|
288 | zrro = 1. |
---|
289 | zrrs = 1. |
---|
290 | zrrc = 1. |
---|
291 | |
---|
292 | zdm0=0. |
---|
293 | zdc0=0. |
---|
294 | zds0=0. |
---|
295 | |
---|
296 | DO jk = jpkb+1,jpkm1 |
---|
297 | DO jj = 1,jpj |
---|
298 | DO ji = 1,jpi |
---|
299 | |
---|
300 | tun = fsdepw(ji,jj,jk+1) |
---|
301 | tob = fsdepw(ji,jj,jk) |
---|
302 | tobbb=fsdepw(ji,jj,jpkb+1) |
---|
303 | zfluo = (1.E2/tob)**fluexp |
---|
304 | zfluu = (1.E2/tun)**fluexp |
---|
305 | IF (zfluo.gt.1.) zfluo = 1. |
---|
306 | zdm0(ji,jj,jk) = zfluo-zfluu |
---|
307 | zdc0(ji,jj,jk) = exp(-(tob-tobbb)/pendec)- |
---|
308 | & exp(-(tun-tobbb)/pendec) |
---|
309 | zds0(ji,jj,jk) = exp(-(tob-tobbb)/pendes)- |
---|
310 | & exp(-(tun-tobbb)/pendes) |
---|
311 | zrrs(ji,jj) = zrrs(ji,jj)-zds0(ji,jj,jk) |
---|
312 | zrro(ji,jj) = zrro(ji,jj)-zdm0(ji,jj,jk) |
---|
313 | zrrc(ji,jj) = zrrc(ji,jj)-zdc0(ji,jj,jk) |
---|
314 | END DO |
---|
315 | END DO |
---|
316 | END DO |
---|
317 | |
---|
318 | DO jj = 1,jpj |
---|
319 | DO ji = 1,jpi |
---|
320 | zdc0(ji,jj,jpk) = zrrc(ji,jj) |
---|
321 | zdm0(ji,jj,jpk) = zrro(ji,jj) |
---|
322 | zds0(ji,jj,jpk) = zrrs(ji,jj) |
---|
323 | END DO |
---|
324 | END DO |
---|
325 | C |
---|
326 | C Calculate vertical distribution of newly formed biogenic poc and |
---|
327 | C caco3 particles in the water column with realistic topography (lowest |
---|
328 | c |
---|
329 | C wet layer contains total fraction, which has passed the upper layers) |
---|
330 | C ---------------------------------------------------------------------- |
---|
331 | C |
---|
332 | dminl = 0. |
---|
333 | dissl = 0. |
---|
334 | discl = 0. |
---|
335 | C |
---|
336 | dmin3 = zdm0 |
---|
337 | diss3 = zds0 |
---|
338 | disc3 = zdc0 |
---|
339 | C |
---|
340 | DO jk = 1,jpk |
---|
341 | DO jj = 1,jpj |
---|
342 | DO ji = 1,jpi |
---|
343 | IF(tmask(ji,jj,jk).eq.0.) THEN |
---|
344 | |
---|
345 | dminl(ji,jj) = dminl(ji,jj)+dmin3(ji,jj,jk) |
---|
346 | dissl(ji,jj) = dissl(ji,jj)+diss3(ji,jj,jk) |
---|
347 | discl(ji,jj) = discl(ji,jj)+disc3(ji,jj,jk) |
---|
348 | dmin3(ji,jj,jk) = 0.0 |
---|
349 | diss3(ji,jj,jk) = 0.0 |
---|
350 | disc3(ji,jj,jk) = 0.0 |
---|
351 | C |
---|
352 | ENDIF |
---|
353 | ENDDO |
---|
354 | ENDDO |
---|
355 | ENDDO |
---|
356 | C |
---|
357 | C |
---|
358 | DO jj = 1,jpj |
---|
359 | DO ji = 1,jpi |
---|
360 | IF (tmask(ji,jj,1).eq.0.) THEN |
---|
361 | dmin3(ji,jj,1) = 0. |
---|
362 | diss3(ji,jj,1) = 0. |
---|
363 | disc3(ji,jj,1) = 0. |
---|
364 | ENDIF |
---|
365 | ENDDO |
---|
366 | ENDDO |
---|
367 | #endif |
---|
368 | C |
---|
369 | C |
---|
370 | CC---------------------------------------------------------------------- |
---|
371 | CC |
---|
372 | CC Initialize chemical variables |
---|
373 | CC |
---|
374 | CC---------------------------------------------------------------------- |
---|
375 | C |
---|
376 | C set pre-industrial atmospheric [co2] (ppm) and o2/n2 ratio |
---|
377 | C ---------------------------------------------------------- |
---|
378 | C |
---|
379 | atcco20 = 278. |
---|
380 | atcox = 0.20946 |
---|
381 | C |
---|
382 | C Set half precision constants |
---|
383 | C ---------------------------- |
---|
384 | C |
---|
385 | tenm7 = 10.**(-7.0) |
---|
386 | smicr = 1.E-6 |
---|
387 | thousi = 1./1000. |
---|
388 | perc = 0.01 |
---|
389 | fourth = 0.25 |
---|
390 | third = 1./3. |
---|
391 | half = 0.5 |
---|
392 | one = 1. |
---|
393 | two = 2. |
---|
394 | ten = 10. |
---|
395 | C |
---|
396 | C Set lower/upper limits for temperature and salinity |
---|
397 | C --------------------------------------------------- |
---|
398 | C |
---|
399 | tfree = -1.9 |
---|
400 | tboil = 30. |
---|
401 | sweet = 20. |
---|
402 | brine = 40. |
---|
403 | C |
---|
404 | salchl = 1./1.80655 |
---|
405 | temzer = 273.16 |
---|
406 | calcon = 1.03E-2 |
---|
407 | C |
---|
408 | C Set coefficients for apparent solubility equilibrium |
---|
409 | C of calcite (Ingle, 1800, eq. 6) |
---|
410 | C ---------------------------------------------------- |
---|
411 | C |
---|
412 | akcc1 = -34.452 |
---|
413 | akcc2 = -39.866 |
---|
414 | akcc3 = 110.21 |
---|
415 | akcc4 = -7.5752E-6 |
---|
416 | C |
---|
417 | arafra = 0. |
---|
418 | calfra = 1.-arafra |
---|
419 | sucall = arafra+calfra |
---|
420 | aracal = arafra*1.45+calfra |
---|
421 | C |
---|
422 | C Set coefficients for seawater pressure correction |
---|
423 | C ------------------------------------------------- |
---|
424 | C |
---|
425 | devk1 = 24.2 |
---|
426 | devk2 = 16.4 |
---|
427 | devkb = 27.5 |
---|
428 | devk1t = 0.085 |
---|
429 | devk2t = 0.04 |
---|
430 | devkbt = 0.095 |
---|
431 | C |
---|
432 | devkst = 0.23 |
---|
433 | devks = 32.8*arafra+35.4*calfra |
---|
434 | C |
---|
435 | C Set universal gas constants |
---|
436 | C --------------------------- |
---|
437 | C |
---|
438 | rgas = 83.143 |
---|
439 | oxyco = 1./22.4144 |
---|
440 | C |
---|
441 | C Set boron constants |
---|
442 | C ------------------- |
---|
443 | C |
---|
444 | bor1 = 0.00023 |
---|
445 | bor2 = 1./10.82 |
---|
446 | C |
---|
447 | C Set volumetric solubility constants for co2 in ml/l (Weiss, 1974) |
---|
448 | C ----------------------------------------------------------------- |
---|
449 | C |
---|
450 | c00 = -58.0931 |
---|
451 | c01 = 90.5069 |
---|
452 | c02 = 22.2940 |
---|
453 | c03 = 0.027766 |
---|
454 | c04 = -0.025888 |
---|
455 | c05 = 0.0050578 |
---|
456 | C |
---|
457 | C Set coeff. for 1. dissoc. of carbonic acid (Edmond and Gieskes, 1970) |
---|
458 | C --------------------------------------------------------------------- |
---|
459 | C |
---|
460 | c10 = -2307.1266 |
---|
461 | c11 = 2.83655 |
---|
462 | c12 = -1.5529413 |
---|
463 | c13 = -4.0484 |
---|
464 | c14 = -0.20760841 |
---|
465 | c15 = 0.08468345 |
---|
466 | c16 = -0.00654208 |
---|
467 | c17 = -0.001005 |
---|
468 | C |
---|
469 | C Set coeff. for 2. dissoc. of carbonic acid (Edmond and Gieskes, 1970) |
---|
470 | C --------------------------------------------------------------------- |
---|
471 | C |
---|
472 | c20 = -3351.6106 |
---|
473 | c21 = -9.226508 |
---|
474 | c22 = -0.2005743 |
---|
475 | c23 = -23.9722 |
---|
476 | c24 = -0.106901773 |
---|
477 | c25 = 0.1130822 |
---|
478 | c26 = -0.00846934 |
---|
479 | c27 = -0.001005 |
---|
480 | C |
---|
481 | C Set coeff. for 1. dissoc. of boric acid (Edmond and Gieskes, 1970) |
---|
482 | C ------------------------------------------------------------------ |
---|
483 | C |
---|
484 | cb0 = -8966.90 |
---|
485 | cb1 = -2890.53 |
---|
486 | cb2 = -77.942 |
---|
487 | cb3 = 1.728 |
---|
488 | cb4 = -0.0996 |
---|
489 | cb5 = 148.0248 |
---|
490 | cb6 = 137.1942 |
---|
491 | cb7 = 1.62142 |
---|
492 | cb8 = -24.4344 |
---|
493 | cb9 = -25.085 |
---|
494 | cb10 = -0.2474 |
---|
495 | cb11 = 0.053105 |
---|
496 | C |
---|
497 | C Set coeff. for dissoc. of water (Dickson and Riley, 1979, |
---|
498 | C eq. 7, coefficient cw2 corrected from 0.9415 to 0.09415 |
---|
499 | C after pers. commun. to B. Bacastow, 1988) |
---|
500 | C --------------------------------------------------------- |
---|
501 | C |
---|
502 | cw0 = -13847.26 |
---|
503 | cw1 = 148.9652 |
---|
504 | cw2 = -23.6521 |
---|
505 | cw3 = 118.67 |
---|
506 | cw4 = -5.977 |
---|
507 | cw5 = 1.0495 |
---|
508 | cw6 = -0.01615 |
---|
509 | C |
---|
510 | C Set volumetric solubility constants for o2 in ml/l (Weiss, 1970) |
---|
511 | C ---------------------------------------------------------------- |
---|
512 | C |
---|
513 | ox0 = -58.3877 |
---|
514 | ox1 = 85.8079 |
---|
515 | ox2 = 23.8439 |
---|
516 | ox3 = -0.034892 |
---|
517 | ox4 = 0.015568 |
---|
518 | ox5 = -0.0019387 |
---|
519 | C |
---|
520 | #if defined key_off_degrad |
---|
521 | C |
---|
522 | C Read volume for degraded regions (DEGINIT) |
---|
523 | C ------------------------------------------ |
---|
524 | C |
---|
525 | # if defined key_vpp |
---|
526 | CALL READ3S(902,facvol,jpi,jpj,jpk) |
---|
527 | # else |
---|
528 | READ (902) facvol |
---|
529 | # endif |
---|
530 | #endif |
---|
531 | C |
---|
532 | C CALCUL DU MASK DE COTE |
---|
533 | C |
---|
534 | cmask=0. |
---|
535 | do ji=2,jpi-1 |
---|
536 | do jj=2,jpj-1 |
---|
537 | if (tmask(ji,jj,1).eq.1) then |
---|
538 | ztest=tmask(ji+1,jj,1)*tmask(ji-1,jj,1)*tmask(ji,jj+1,1) |
---|
539 | . *tmask(ji,jj-1,1) |
---|
540 | if (ztest.eq.0) cmask(ji,jj)=1. |
---|
541 | endif |
---|
542 | end do |
---|
543 | end do |
---|
544 | |
---|
545 | cmask(1,:)=cmask(jpi-1,:) |
---|
546 | cmask(jpi,:)=cmask(2,:) |
---|
547 | C |
---|
548 | C CALCUL DE LA SURFACE COTIERE |
---|
549 | C |
---|
550 | do ji=2,jpi-1 |
---|
551 | do jj=2,jpj-1 |
---|
552 | areacot=areacot+e1t(ji,jj)*e2t(ji,jj)*cmask(ji,jj) |
---|
553 | end do |
---|
554 | end do |
---|
555 | |
---|
556 | C |
---|
557 | C |
---|
558 | CALL h3cche |
---|
559 | WRITE(numout,*) ' Initialisation of HAMOCC3 done' |
---|