1 | CCC $Header$ |
---|
2 | CCC TOP 1.0 , LOCEAN-IPSL (2005) |
---|
3 | C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
4 | C --------------------------------------------------------------------------- |
---|
5 | CDIR$ LIST |
---|
6 | SUBROUTINE h3cbio |
---|
7 | #if defined key_passivetrc && defined key_trc_hamocc3 |
---|
8 | CCC--------------------------------------------------------------------- |
---|
9 | CCC |
---|
10 | CCC ROUTINE h3cbio |
---|
11 | CCC ****************** |
---|
12 | CCC |
---|
13 | CCC |
---|
14 | CCC PURPOSE. |
---|
15 | CCC -------- |
---|
16 | CCC *H3CBIO* MODELS PRODUCTION OF BIOGENIC MATTER (POC ''SOFT |
---|
17 | CCC TISSUE'' AND CACO3 PARTICLES ''HARD PARTS'') |
---|
18 | CCC AND ITS DISTRIBUTION IN WATER COLUMN |
---|
19 | CCC |
---|
20 | CC |
---|
21 | CC METHOD. |
---|
22 | CC ------- |
---|
23 | CC IN THE SURFACE LAYER POC IS PRODUCED ACCORDING TO |
---|
24 | CC NURTRIENTS AVAILABLE AND GROWTH CONDITIONS. NUTRIENT UPTAKE |
---|
25 | CC KINETICS FOLLOW MICHAELIS-MENTON FORMULATION. PROPORTIONAL |
---|
26 | CC TO THE AMOUNT OF ORGANIC MATTER, CACO3 HARD PARTS ARE PRODUCED. |
---|
27 | CC THE TOTAL PARTICLE AMOUNT PRODUCED, IS DISTRIBUTED IN THE WATER |
---|
28 | CC COLUMN BELOW THE SURFACE LAYER. |
---|
29 | CC |
---|
30 | CC EXTERNALS. |
---|
31 | CC ---------- |
---|
32 | CC NONE. |
---|
33 | CC |
---|
34 | CC REFERENCE. |
---|
35 | CC ---------- |
---|
36 | CC |
---|
37 | CC BACASTOW, R., AND E. MAIER-REIMER (1985) |
---|
38 | CC CIRCULATION MODEL OF THE OCEAN CARBON CYCLE. |
---|
39 | CC 1. DESCRIPTION OF THE MODEL, PP. 224-232. |
---|
40 | CC 2. COMPARISON OF THE MODEL RESULTS WITH OBSERVATIONAL DATA, |
---|
41 | CC PP. 233-240. |
---|
42 | CC IN: "ATMOSPHERIC CARBON DIOXIDE - ITS SOURCES, SINKS, AND |
---|
43 | CC GLOBAL TRANSPORT", KANDERSTEG, 2 TO 6 SEPTEMBER 1985, |
---|
44 | CC COMMISSION ON ATMOSPHERIC CHEMISTRY AND GLOBAL POLLUTION, |
---|
45 | CC INTERNATIONAL ASSOCIATION OF METEOROLOGY AND ATMOSPHERIC PHYSICS. |
---|
46 | CC |
---|
47 | CC DUGDALE. R.C. (1967) |
---|
48 | CC NUTRIENT LIMITATION IN THE SEA: DYNAMICS, IDENTIFICATION |
---|
49 | CC AND SIGNIFICANCE. |
---|
50 | CC LIMNOLOGY AND OCEANOGRAPHY, VOL.12, 685-695. |
---|
51 | CC |
---|
52 | CC PARSONS, T.R., AND M. TAKAHASHI (1973) |
---|
53 | CC BIOLOGICAL OCEANOGRAPHIC PROCESSES. |
---|
54 | CC PERGAMON PRESS, 186 PP. |
---|
55 | CC |
---|
56 | CC MODIFICATIONS: |
---|
57 | CC -------------- |
---|
58 | CC original : 1988-07 E. MAIER-REIMER MPI HAMBURG |
---|
59 | CC additions : 1998 O. Aumont |
---|
60 | CC modifications : 1999 C. Le Quere |
---|
61 | CC modifications : 2001 O. Aumont |
---|
62 | CC --------------------------------------------------------------------------- |
---|
63 | CC parameters and commons |
---|
64 | CC ====================== |
---|
65 | CDIR$ NOLIST |
---|
66 | USE oce_trc |
---|
67 | USE trp_trc |
---|
68 | USE sms |
---|
69 | IMPLICIT NONE |
---|
70 | CDIR$ LIST |
---|
71 | CC---------------------------------------------------------------------- |
---|
72 | CC local declarations |
---|
73 | CC ================== |
---|
74 | C |
---|
75 | INTEGER ji,jj,jk |
---|
76 | REAL propro, brupro, silpot, expofa |
---|
77 | REAL calred, calpot, topot, silfra |
---|
78 | C |
---|
79 | C ============================================================ |
---|
80 | C* 2. CALCULATE PRODUCTION OF SOFT TISSUE AND RELATED CACO3 |
---|
81 | C HARD PARTS USING MICHAELIS MENTON KINETICS (CF. DUGDALE, |
---|
82 | C 1967; PARSONS AND TAKAHASHI, 1973) |
---|
83 | C ============================================================ |
---|
84 | C |
---|
85 | C* 1. SET HALF PRECISION CONSTANTS |
---|
86 | C ----------------------------------------------------------------- |
---|
87 | C |
---|
88 | zero = 0. |
---|
89 | one = 1. |
---|
90 | two = 2. |
---|
91 | C |
---|
92 | CALL h3copt |
---|
93 | C |
---|
94 | C vertical slab |
---|
95 | C =============== |
---|
96 | C |
---|
97 | C |
---|
98 | DO jk = 1,jpkb |
---|
99 | DO jj = 1,jpj |
---|
100 | DO ji = 1,jpi |
---|
101 | C |
---|
102 | C* 2.2 PRIMARY PRODUCTION OF POC (SOFT TISSUE) USING |
---|
103 | C MICHAELIS-MENTON NUTRIENT UPTAKE KINETICS [MOLE/L] |
---|
104 | C -------------------------------------------------------- |
---|
105 | C |
---|
106 | propro = strn(ji,jj)*0.25 |
---|
107 | IF (propro.le.0.) propro = 0. |
---|
108 | C |
---|
109 | C for T lt -2. no production |
---|
110 | C --------------------------- |
---|
111 | C |
---|
112 | IF (tn(ji,jj,jk).lt.-2.) THEN |
---|
113 | prlat(ji,jj) = 0. |
---|
114 | ELSE |
---|
115 | prlat(ji,jj) = propro*((tn(ji,jj,jk)+2.)/(tn(ji,jj,jk)+10.)) |
---|
116 | $ *tmask(ji,jj,jk)*(1.-freeze(ji,jj)) |
---|
117 | C |
---|
118 | # if defined key_off_degrad |
---|
119 | $ *facvol(ji,jj,jk) |
---|
120 | # endif |
---|
121 | ENDIF |
---|
122 | C |
---|
123 | IF (hmld(ji,jj).gt.50.) THEN |
---|
124 | prlat(ji,jj) = prlat(ji,jj)*(50./hmld(ji,jj)) |
---|
125 | ENDIF |
---|
126 | C |
---|
127 | C* 2.3 POC PRODUCTION RELATED BIOGENIC CACO3 PARTICLE |
---|
128 | C PRODUCTION [MOLE/L] |
---|
129 | C ---------------------------------------------------- |
---|
130 | C |
---|
131 | prorca(ji,jj,jk) = prlat(ji,jj)*trn(ji,jj,jk,jppo4)**2 |
---|
132 | & /(conc0+trn(ji,jj,jk,jppo4))*rfact |
---|
133 | |
---|
134 | brupro = grosip*prorca(ji,jj,jk) |
---|
135 | tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4)- |
---|
136 | & prorca(ji,jj,jk)*rfactr |
---|
137 | tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy)+ |
---|
138 | & o2ut*prorca(ji,jj,jk)*rfactr |
---|
139 | silpot = min(brupro,trn(ji,jj,jk,jpsil)) |
---|
140 | expofa = exp(0.1*(tn(ji,jj,jk)-10.)) |
---|
141 | calred = expofa/(1.+expofa) |
---|
142 | calpot = calred*prorca(ji,jj,jk) |
---|
143 | topot = silpot+calpot+1.E-20 |
---|
144 | silfra = silpot/topot |
---|
145 | C ENDE ERSATZ |
---|
146 | prcaca(ji,jj,jk) = caco3r*calpot*silic0/(silic0+silpot) |
---|
147 | silpro(ji,jj,jk) = silpot*silfra |
---|
148 | $ *trn(ji,jj,jk,jpsil)/(trn(ji,jj,jk,jpsil)+silic0) |
---|
149 | C |
---|
150 | C* 5.1 SURFACE LAYER C13/C12, C14/C12 TOTAL CO2 RATIOS |
---|
151 | C ----------------------------------------------------- |
---|
152 | C |
---|
153 | tra(ji,jj,jk,jpsil) = tra(ji,jj,jk,jpsil)- |
---|
154 | & silpro(ji,jj,jk)*rfactr |
---|
155 | C |
---|
156 | C* 5.2 CONSUMPTION OF TOTAL (12C)O2 |
---|
157 | C ---------------------------------- |
---|
158 | C |
---|
159 | tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - |
---|
160 | & (prorca(ji,jj,jk)+prcaca(ji,jj,jk)) * rfactr |
---|
161 | C |
---|
162 | C* 5.3 CONSUMPTION OF ALKALINITY DUE TO CA++ UPTAKE |
---|
163 | C AND INCREASE OF ALKALINITY DUE TO NITRATE |
---|
164 | C CONSUMPTION DURING ORGANIC SOFT TISSUE PRODUCTION |
---|
165 | C ------------------------------------------------------- |
---|
166 | C |
---|
167 | tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal)+ |
---|
168 | & (RNO3*PRORCA(ji,jj,jk)-TWO*PRCACA(ji,jj,jk))*rfactr |
---|
169 | C |
---|
170 | C* 5.4 CONSUMPTION OF PHOSPHATE |
---|
171 | C PHOSPHATE REDUCTION HAPPENS IN SUBROUTINE ADVECT |
---|
172 | c trn(ji,jj,jk,jppo4) = trn(ji,jj,jk,jppo4)-PRORCA(ji,jj,jk) |
---|
173 | C ------------------------------------------------------ |
---|
174 | C |
---|
175 | # if defined key_trc_biohamocc13 |
---|
176 | C |
---|
177 | C* 5.6 CONSUMPTION OF TOTAL (13C)O2 |
---|
178 | C ---------------------------------- |
---|
179 | C |
---|
180 | tra(ji,jj,jk,jp13c) = tra(ji,jj,jk,jp13c) - |
---|
181 | & (pdb*plafr13*prorca(ji,jj,jk)+pdb*prcaca(ji,jj,jk))*rfactr |
---|
182 | C |
---|
183 | # endif |
---|
184 | C |
---|
185 | ENDDO |
---|
186 | END DO |
---|
187 | END DO |
---|
188 | C |
---|
189 | # if defined key_trc_diaadd |
---|
190 | do jj=1,jpj |
---|
191 | do ji=1,jpi |
---|
192 | trc2d(ji,jj,8)=tn(ji,jj,1) |
---|
193 | trc2d(ji,jj,9)=hmld(ji,jj) |
---|
194 | end do |
---|
195 | end do |
---|
196 | |
---|
197 | # endif |
---|
198 | #endif |
---|
199 | RETURN |
---|
200 | END |
---|
201 | |
---|