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 | SUBROUTINE p3zrem |
---|
8 | #if defined key_passivetrc && key_trc_p3zd |
---|
9 | CCC--------------------------------------------------------------------- |
---|
10 | CCC |
---|
11 | CCC ROUTINE p3zrem |
---|
12 | CCC ****************** |
---|
13 | CCC |
---|
14 | CCC |
---|
15 | CCC |
---|
16 | CCC |
---|
17 | CCC PURPOSE. |
---|
18 | CCC -------- |
---|
19 | CC |
---|
20 | CC *P3ZREM* CALCULATES OXYGEN CONSUMPTION DURING REMINERALIZATION |
---|
21 | CC AND LOSS OF POC TO THE ORGANIC SEDIMENT POOL |
---|
22 | CC |
---|
23 | CC* INTERFACE. |
---|
24 | CC ---------- |
---|
25 | CC *CALL* *P3ZREM* |
---|
26 | CC |
---|
27 | CC METHOD. |
---|
28 | CC ------- |
---|
29 | CC THE OXYGEN AMOUNT AVAILABLE FOR CONSUMPTION DUE TO |
---|
30 | CC REMINERALIZATION OF POC IS DETERMINED. |
---|
31 | CC THE CORRESPONDING DECREASE OF POC AND THE INCREASES OF |
---|
32 | CC INORGANIC CARBON AND PHOSPHATE ARE CALCULATED. ALKALINITY |
---|
33 | CC IS REDUCED BY THE NITRATE SET FREE DURING REMINERALIZATION. |
---|
34 | CC POC SETTLES DOWN TO THE SEDIMENT FROM THE WATER LAYER |
---|
35 | CC DIRECTLY ABOVE THE OCEAN FLOOR (DEPENDS ON BOTTOM DEPTH) |
---|
36 | CC AT A RATE TO BE DEFINED (''OCSED''). |
---|
37 | CC |
---|
38 | CC EXTERNALS. |
---|
39 | CC ---------- |
---|
40 | CC NONE. |
---|
41 | CC |
---|
42 | CC REFERENCE. |
---|
43 | CC ---------- |
---|
44 | CC -. |
---|
45 | CC |
---|
46 | CC VARIABLE TYPE PURPOSE. |
---|
47 | CC -------- ---- -------- |
---|
48 | CC |
---|
49 | CC *OLIMI(IE,JE)*REAL MOLES O2 AVAILABLE FOR REMINERALIZATION |
---|
50 | CC (LIMITING INGREDIENT FOR REMINERALIZATION) |
---|
51 | CC *KE1* INTEGER MAX. OF INDEX FOR OCEAN LAYERS MINUS 1 |
---|
52 | CC *OCSED* REAL FRACTION OF POC LOST TO ORGANIC CARBON SEDIMENT |
---|
53 | CC POOL FROM WATER LAYER IMMEDIATELY ABOVE OCEAN |
---|
54 | CC FLOOR |
---|
55 | CC *OCSEDI* REAL FRACTION OF POC REMAINING IN THE WATER LAYER |
---|
56 | CC IMMEDIATELY ABOVE THE BOTTOM AND DOES NOT |
---|
57 | CC ENTER THE ORG. CARBON SEDIMENT POOL |
---|
58 | CC *REMIK* REAL REMINERALIZATION RATE (FRACTION OF POC, |
---|
59 | CC THAT IS REMINERALISED DURING ONE TIMESTEP), |
---|
60 | CC DUMMY VARIABLE |
---|
61 | CC |
---|
62 | CC MODIFICATIONS: |
---|
63 | CC -------------- |
---|
64 | CC original : 1988-07 E. MAIER-REIMER MPI HAMBURG |
---|
65 | CC additions : 1998 O. Aumont |
---|
66 | CC modifications : 1999 C. Le Quere |
---|
67 | CC ---------------------------------------------------------------------- |
---|
68 | CC parameters and commons |
---|
69 | CC ====================== |
---|
70 | CDIR$ NOLIST |
---|
71 | USE oce_trc |
---|
72 | USE trp_trc |
---|
73 | USE sms |
---|
74 | IMPLICIT NONE |
---|
75 | CDIR$ LIST |
---|
76 | CC---------------------------------------------------------------------- |
---|
77 | CC local declarations |
---|
78 | CC ================== |
---|
79 | C |
---|
80 | INTEGER ji, jj, jk, jn, ke1 |
---|
81 | REAL zrem, zorem, zremip, zremik, zphy, zzoo |
---|
82 | C |
---|
83 | C |
---|
84 | C REMINERALIZATION OF POC AND O2 UTILIZATION |
---|
85 | C ------------------------------------------ |
---|
86 | C |
---|
87 | DO jk = 11,jpk-1 |
---|
88 | C |
---|
89 | C SET FRACTION POC TO BE REMINERALISED IN LAYER K |
---|
90 | C ----------------------------------------------- |
---|
91 | C |
---|
92 | DO jj = 1,jpj |
---|
93 | DO ji = 1,jpi |
---|
94 | C |
---|
95 | zremik=50.*spocri*(1.-tmask(ji,jj,jk+1))+xremik |
---|
96 | zremik = zremik/rjjss*rfact*tmask(ji,jj,jk) |
---|
97 | & *trn(ji,jj,jk,jppo4)/(trn(ji,jj,jk,jppo4)+xkdoc1) |
---|
98 | & *trn(ji,jj,jk,jpdoc)/(trn(ji,jj,jk,jpdoc)+xkdoc2) |
---|
99 | # if defined key_off_degrad |
---|
100 | & *facvol(ji,jj,jk) |
---|
101 | # endif |
---|
102 | zremip = xremip/rjjss*rfact*tmask(ji,jj,jk) |
---|
103 | # if defined key_off_degrad |
---|
104 | & *facvol(ji,jj,jk) |
---|
105 | # endif |
---|
106 | zzoo = (resrat+mzrat)/rjjss*rfact*tmask(ji,jj,jk)* |
---|
107 | & trn(ji,jj,jk,jpzoo)*trn(ji,jj,jk,jpzoo)/ |
---|
108 | & (1.E-6+trn(ji,jj,jk,jpzoo)) |
---|
109 | # if defined key_off_degrad |
---|
110 | & *facvol(ji,jj,jk) |
---|
111 | # endif |
---|
112 | zphy = mprat/rjjss*rfact*tmask(ji,jj,jk)* |
---|
113 | & trn(ji,jj,jk,jpphy) |
---|
114 | & *trn(ji,jj,jk,jpphy)/(1.E-6+trn(ji,jj,jk,jpphy)) |
---|
115 | # if defined key_off_degrad |
---|
116 | & *facvol(ji,jj,jk) |
---|
117 | # endif |
---|
118 | C |
---|
119 | C AMOUNT OF O2 AVAILABLE FOR REMINERALIZATION OF POC |
---|
120 | C |
---|
121 | zrem = zremik* trn(ji,jj,jk,jpdoc) |
---|
122 | zorem = zremip*trn(ji,jj,jk,jppoc) |
---|
123 | C |
---|
124 | C |
---|
125 | C INCREASE OF PO4 AND SUM((12C)O2) DUE TO RELEASE TO WATER |
---|
126 | C COLUMN DURING REMINERALIZATION |
---|
127 | C -------------------------------------------------------- |
---|
128 | C |
---|
129 | trn(ji,jj,jk,jpzoo) = trn(ji,jj,jk,jpzoo)-zzoo |
---|
130 | trn(ji,jj,jk,jpphy) = trn(ji,jj,jk,jpphy)-zphy |
---|
131 | trn(ji,jj,jk,jppo4) = trn(ji,jj,jk,jppo4)+zrem |
---|
132 | trn(ji,jj,jk,jpdic) = trn(ji,jj,jk,jpdic)+zrem |
---|
133 | C |
---|
134 | C DECREASE OF ALKALINITY AND O2 DUE TO REMINERALIZATION |
---|
135 | C (ALKALINITY DECREASES DUE TO NITRATE RELEASE) |
---|
136 | C ----------------------------------------------------- |
---|
137 | C |
---|
138 | trn(ji,jj,jk,jptal) = trn(ji,jj,jk,jptal)-zrem*rno3 |
---|
139 | trn(ji,jj,jk,jpoxy) = trn(ji,jj,jk,jpoxy)-o2ut*zrem |
---|
140 | trn(ji,jj,jk,jpoxy) = MAX(4.E-6,trn(ji,jj,jk,jpoxy)) |
---|
141 | C |
---|
142 | # if defined key_trc_biohamocc13 |
---|
143 | C CHANGE IN RARE ISOTOPES OF INORGANIC CARBON AND POC |
---|
144 | C AFTER REMINERALIZATION AND DECREASE OF POC 12C |
---|
145 | C IN THE WATER COLUMN DUE TO REMINERALIZATION |
---|
146 | C |
---|
147 | trn(ji,jj,jk,jp13c) = trn(ji,jj,jk,jp13c)+zrem*pdb* |
---|
148 | & plafr13 |
---|
149 | # endif |
---|
150 | C |
---|
151 | trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc)-zorem |
---|
152 | trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc)+zorem+zphy |
---|
153 | & +zzoo-zrem |
---|
154 | C |
---|
155 | ENDDO |
---|
156 | ENDDO |
---|
157 | ENDDO |
---|
158 | C |
---|
159 | #endif |
---|
160 | C |
---|
161 | RETURN |
---|
162 | END |
---|