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 | CDIR$ LIST |
---|
8 | SUBROUTINE h3crem |
---|
9 | #if defined key_passivetrc && defined key_trc_hamocc3 |
---|
10 | CCC--------------------------------------------------------------------- |
---|
11 | CCC |
---|
12 | CCC ROUTINE h3crem |
---|
13 | CCC ****************** |
---|
14 | CCC |
---|
15 | CCC |
---|
16 | CCC |
---|
17 | CCC |
---|
18 | CCC PURPOSE. |
---|
19 | CCC -------- |
---|
20 | CCC *H3CREM* CALCULATES OXYGEN CONSUMPTION DURING REMINERALIZATION |
---|
21 | CCC AND LOSS OF POC TO THE ORGANIC SEDIMENT POOL |
---|
22 | CCC |
---|
23 | CC |
---|
24 | CC METHOD. |
---|
25 | CC ------- |
---|
26 | CC THE OXYGEN AMOUNT AVAILABLE FOR CONSUMPTION DUE TO |
---|
27 | CC REMINERALIZATION OF POC IS DETERMINED. |
---|
28 | CC THE CORRESPONDING DECREASE OF POC AND THE INCREASES OF |
---|
29 | CC INORGANIC CARBON AND PHOSPHATE ARE CALCULATED. ALKALINITY |
---|
30 | CC IS REDUCED BY THE NITRATE SET FREE DURING REMINERALIZATION. |
---|
31 | CC POC SETTLES DOWN TO THE SEDIMENT FROM THE WATER LAYER |
---|
32 | CC DIRECTLY ABOVE THE OCEAN FLOOR (DEPENDS ON BOTTOM DEPTH) |
---|
33 | CC AT A RATE TO BE DEFINED (''OCSED''). |
---|
34 | CC |
---|
35 | CC EXTERNALS. |
---|
36 | CC ---------- |
---|
37 | CC NONE. |
---|
38 | CC |
---|
39 | CC REFERENCE. |
---|
40 | CC ---------- |
---|
41 | CC -. |
---|
42 | CC |
---|
43 | CC* VARIABLE TYPE PURPOSE. |
---|
44 | CC -------- ---- -------- |
---|
45 | CC |
---|
46 | CC *ZREM* REAL MOLES O2 AVAILABLE FOR REMINERALIZATION |
---|
47 | CC (LIMITING INGREDIENT FOR REMINERALIZATION) |
---|
48 | CC *ZREMIK* REAL REMINERALIZATION RATE (FRACTION OF POC, |
---|
49 | CC THAT IS REMINERALISED DURING ONE TIMESTEP) |
---|
50 | CC DUMMY VARIABLE |
---|
51 | CC |
---|
52 | CC MODIFICATIONS: |
---|
53 | CC -------------- |
---|
54 | CC original : 1988-07 E. MAIER-REIMER MPI HAMBURG |
---|
55 | CC additions : 1998 O. Aumont |
---|
56 | CC modifications : 1999 C. Le Quere |
---|
57 | CC --------------------------------------------------------------------- |
---|
58 | CC parameters and commons |
---|
59 | CC ====================== |
---|
60 | CDIR$ NOLIST |
---|
61 | USE oce_trc |
---|
62 | USE trp_trc |
---|
63 | USE sms |
---|
64 | IMPLICIT NONE |
---|
65 | CDIR$ LIST |
---|
66 | CC---------------------------------------------------------------------- |
---|
67 | CC local declarations |
---|
68 | CC ================== |
---|
69 | C |
---|
70 | INTEGER ji, jj, jk |
---|
71 | REAL zremik, zrem |
---|
72 | C |
---|
73 | C* 3. REMINERALIZATION OF POC AND O2 UTILIZATION |
---|
74 | C ---------------------------------------------- |
---|
75 | C |
---|
76 | C* 3.1 LAYER LOOP |
---|
77 | C ---------------- |
---|
78 | C |
---|
79 | DO jk = 1,jpkm1 |
---|
80 | DO jj = 1,jpj |
---|
81 | C |
---|
82 | C* 3.2 SET FRACTION POC TO BE REMINERALISED IN LAYER K |
---|
83 | C ----------------------------------------------------- |
---|
84 | C |
---|
85 | C* 3.3 MERIDIONAL/ZONAL LOOPS |
---|
86 | C ---------------------------- |
---|
87 | C |
---|
88 | DO ji = 1,jpi |
---|
89 | C |
---|
90 | zrem = 0. |
---|
91 | C |
---|
92 | zremik = 50.*spocri*(1.-tmask(ji,jj,jk+1)) |
---|
93 | $ +0.002*(tn(ji,jj,jk)+2.) |
---|
94 | zremik = 0.002*(tn(ji,jj,jk)+2.) |
---|
95 | zremik = zremik/(30.5*3600.*24.)*rfact |
---|
96 | # if defined key_off_degrad |
---|
97 | & *facvol(ji,jj,jk) |
---|
98 | # endif |
---|
99 | C |
---|
100 | C* 3.5 AMOUNT OF O2 AVAILABLE FOR REMINERALIZATION OF POC |
---|
101 | C -------------------------------------------------------- |
---|
102 | C |
---|
103 | zrem = zremik* min( |
---|
104 | & (trn(ji,jj,jk,jpoxy)-oxymin)/o2ut,trn(ji,jj,jk,jppoc)) |
---|
105 | C |
---|
106 | C |
---|
107 | C* 3.6 INCREASE OF PO4 AND SUM((12C)O2) DUE TO RELEASE TO WATER |
---|
108 | C COLUMN DURING REMINERALIZATION |
---|
109 | C -------------------------------------------------------------- |
---|
110 | C |
---|
111 | tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4)+ zrem*rfactr |
---|
112 | tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic)+ zrem*rfactr |
---|
113 | C |
---|
114 | C* 3.7 DECREASE OF ALKALINITY AND O2 DUE TO REMINERALIZATION |
---|
115 | C (ALKALINITY DECREASES DUE TO NITRATE RELEASE) |
---|
116 | C ----------------------------------------------------------- |
---|
117 | C |
---|
118 | tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal)- |
---|
119 | & zrem*rno3*rfactr |
---|
120 | tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy)- |
---|
121 | & rfactr*o2ut*zrem |
---|
122 | C |
---|
123 | C* 3.8 CHANGE IN RARE ISOTOPES OF INORGANIC CARBON AND POC |
---|
124 | C AFTER REMINERALIZATION AND DECREASE OF POC 12C |
---|
125 | C IN THE WATER COLUMN DUE TO REMINERALIZATION |
---|
126 | C ---------------------------------------------------------------------- |
---|
127 | C |
---|
128 | # if defined key_trc_biohamocc13 |
---|
129 | tra(ji,jj,jk,jp13c) = tra(ji,jj,jk,jp13c)+ |
---|
130 | & zrem*pdb*rfactr*plafr13 |
---|
131 | # endif |
---|
132 | tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc)- |
---|
133 | & zrem*rfactr |
---|
134 | C |
---|
135 | ENDDO |
---|
136 | ENDDO |
---|
137 | C |
---|
138 | C END of slab |
---|
139 | C =========== |
---|
140 | C |
---|
141 | ENDDO |
---|
142 | # endif |
---|
143 | RETURN |
---|
144 | END |
---|
145 | |
---|