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