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