[341] | 1 | |
---|
[725] | 2 | CCC $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/SMS/p4zrem.F,v 1.8 2007/10/12 09:28:41 opalod Exp $ |
---|
[341] | 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 --------------------------------------------------------------------------- |
---|
[186] | 6 | CDIR$ LIST |
---|
| 7 | SUBROUTINE p4zrem |
---|
[772] | 8 | #if defined key_top && defined key_pisces |
---|
[186] | 9 | CCC--------------------------------------------------------------------- |
---|
| 10 | CCC |
---|
| 11 | CCC ROUTINE p4zrem : PISCES MODEL |
---|
| 12 | CCC ***************************** |
---|
| 13 | CCC |
---|
| 14 | CCC PURPOSE : |
---|
| 15 | CCC --------- |
---|
| 16 | CCC Compute remineralization/scavenging of organic compounds |
---|
| 17 | CCC |
---|
| 18 | CC INPUT : |
---|
| 19 | CC ----- |
---|
| 20 | CC common |
---|
| 21 | CC all the common defined in opa |
---|
| 22 | CC |
---|
| 23 | CC |
---|
| 24 | CC OUTPUT : : no |
---|
| 25 | CC ------ |
---|
| 26 | CC |
---|
| 27 | CC EXTERNAL : |
---|
| 28 | CC -------- |
---|
| 29 | CC None |
---|
| 30 | CC |
---|
| 31 | CC MODIFICATIONS: |
---|
| 32 | CC -------------- |
---|
| 33 | CC original : 2004 - O. Aumont |
---|
| 34 | CC---------------------------------------------------------------------- |
---|
| 35 | CC parameters and commons |
---|
| 36 | CC ====================== |
---|
| 37 | CDIR$ NOLIST |
---|
| 38 | USE oce_trc |
---|
| 39 | USE trp_trc |
---|
| 40 | USE sms |
---|
| 41 | IMPLICIT NONE |
---|
[339] | 42 | #include "domzgr_substitute.h90" |
---|
[186] | 43 | CDIR$ LIST |
---|
| 44 | CC---------------------------------------------------------------------- |
---|
| 45 | CC local declarations |
---|
| 46 | CC ================== |
---|
| 47 | INTEGER ji, jj, jk |
---|
| 48 | REAL remip,remik,xlam1b |
---|
| 49 | REAL xkeq,xfeequi,siremin |
---|
[339] | 50 | REAL zsatur,zsatur2,znusil,zdepbac(jpi,jpj,jpk) |
---|
| 51 | REAL zlamfac,zstep,fesatur(jpi,jpj,jpk) |
---|
[725] | 52 | REAL ztempbac(jpi,jpj) |
---|
[186] | 53 | C |
---|
[339] | 54 | C Time step duration for the biology |
---|
| 55 | C |
---|
| 56 | zstep=rfact2/rjjss |
---|
| 57 | C |
---|
[186] | 58 | C Computation of the mean phytoplankton concentration as |
---|
| 59 | C a crude estimate of the bacterial biomass |
---|
| 60 | C -------------------------------------------------- |
---|
| 61 | C |
---|
[725] | 62 | DO jk=1,jpkm1 |
---|
| 63 | DO jj = 1, jpj |
---|
| 64 | DO ji = 1, jpi |
---|
| 65 | IF (fsdept(ji,jj,jk).lt.120.) THEN |
---|
| 66 | zdepbac(ji,jj,jk)=min(0.7*(trn(ji,jj,jk,jpzoo) |
---|
| 67 | & +2*trn(ji,jj,jk,jpmes)),4E-6) |
---|
| 68 | ztempbac(ji,jj)=zdepbac(ji,jj,jk) |
---|
| 69 | ELSE |
---|
| 70 | zdepbac(ji,jj,jk)=min(1.,120./fsdept(ji,jj,jk)) |
---|
| 71 | & *ztempbac(ji,jj) |
---|
| 72 | ENDIF |
---|
| 73 | END DO |
---|
| 74 | END DO |
---|
[339] | 75 | END DO |
---|
[186] | 76 | |
---|
[339] | 77 | DO jk = 1,jpkm1 |
---|
[186] | 78 | DO jj = 1,jpj |
---|
| 79 | DO ji = 1,jpi |
---|
| 80 | C |
---|
| 81 | C DENITRIFICATION FACTOR COMPUTED FROM O2 LEVELS |
---|
| 82 | C ---------------------------------------------- |
---|
| 83 | C |
---|
| 84 | nitrfac(ji,jj,jk)= |
---|
| 85 | & max(0.,0.4*(6.E-6-trn(ji,jj,jk,jpoxy))/(oxymin+ |
---|
| 86 | & trn(ji,jj,jk,jpoxy))) |
---|
| 87 | END DO |
---|
| 88 | END DO |
---|
| 89 | END DO |
---|
| 90 | |
---|
[339] | 91 | nitrfac(:,:,:)=min(1.,nitrfac(:,:,:)) |
---|
| 92 | |
---|
[186] | 93 | DO jk = 1,jpkm1 |
---|
| 94 | DO jj = 1,jpj |
---|
| 95 | DO ji = 1,jpi |
---|
| 96 | C |
---|
| 97 | C DOC ammonification. Depends on depth, phytoplankton biomass |
---|
| 98 | C and a limitation term which is supposed to be a parameterization |
---|
| 99 | C of the bacterial activity. |
---|
| 100 | C ---------------------------------------------------------------- |
---|
[339] | 101 | C |
---|
| 102 | remik = xremik*zstep/1E-6*xlimbac(ji,jj,jk) |
---|
| 103 | & *zdepbac(ji,jj,jk) |
---|
[186] | 104 | # if defined key_off_degrad |
---|
| 105 | & *facvol(ji,jj,jk) |
---|
| 106 | # endif |
---|
[339] | 107 | remik=max(remik,5.5E-4*zstep) |
---|
[186] | 108 | C |
---|
| 109 | C Ammonification in oxic waters with oxygen consumption |
---|
| 110 | C ----------------------------------------------------- |
---|
| 111 | C |
---|
| 112 | olimi(ji,jj,jk)=min((trn(ji,jj,jk,jpoxy)-rtrn)/o2ut, |
---|
| 113 | & remik*(1.-nitrfac(ji,jj,jk))*trn(ji,jj,jk,jpdoc)) |
---|
| 114 | C |
---|
| 115 | C Ammonification in suboxic waters with denitrification |
---|
| 116 | C ------------------------------------------------------- |
---|
| 117 | C |
---|
[339] | 118 | denitr(ji,jj,jk)=min((trn(ji,jj,jk,jpno3)-rtrn)/rdenit, |
---|
[186] | 119 | & remik*nitrfac(ji,jj,jk)*trn(ji,jj,jk,jpdoc)) |
---|
| 120 | END DO |
---|
| 121 | END DO |
---|
| 122 | END DO |
---|
[339] | 123 | C |
---|
| 124 | olimi(:,:,:)=max(0.,olimi(:,:,:)) |
---|
| 125 | denitr(:,:,:)=max(0.,denitr(:,:,:)) |
---|
| 126 | C |
---|
[186] | 127 | DO jk = 1,jpkm1 |
---|
| 128 | DO jj = 1,jpj |
---|
| 129 | DO ji = 1,jpi |
---|
| 130 | C |
---|
| 131 | C NH4 nitrification to NO3. Ceased for oxygen concentrations |
---|
| 132 | C below 2 umol/L. Inhibited at strong light |
---|
| 133 | C ---------------------------------------------------------- |
---|
| 134 | C |
---|
[339] | 135 | onitr(ji,jj,jk)=nitrif*zstep*trn(ji,jj,jk,jpnh4)/(1. |
---|
| 136 | & +emoy(ji,jj,jk))*(1.-nitrfac(ji,jj,jk)) |
---|
[186] | 137 | # if defined key_off_degrad |
---|
| 138 | & *facvol(ji,jj,jk) |
---|
| 139 | # endif |
---|
| 140 | END DO |
---|
| 141 | END DO |
---|
| 142 | END DO |
---|
| 143 | |
---|
| 144 | DO jk = 1,jpkm1 |
---|
| 145 | DO jj = 1,jpj |
---|
| 146 | DO ji = 1,jpi |
---|
| 147 | C |
---|
| 148 | C Bacterial uptake of iron. No iron is available in DOC. So |
---|
| 149 | C Bacteries are obliged to take up iron from the water. Some |
---|
| 150 | C studies (especially at Papa) have shown this uptake to be |
---|
| 151 | C significant |
---|
| 152 | C ---------------------------------------------------------- |
---|
| 153 | C |
---|
[339] | 154 | xbactfer(ji,jj,jk)=15E-6*rfact2*4.*0.4*prmax(ji,jj,jk) |
---|
| 155 | & *(xlimphy(ji,jj,jk)*zdepbac(ji,jj,jk))**2 |
---|
| 156 | & /(xkgraz2+zdepbac(ji,jj,jk)) |
---|
| 157 | & *(0.5+sign(0.5,trn(ji,jj,jk,jpfer)-2E-11)) |
---|
| 158 | C |
---|
[186] | 159 | END DO |
---|
| 160 | END DO |
---|
| 161 | END DO |
---|
| 162 | C |
---|
| 163 | DO jk = 1,jpkm1 |
---|
| 164 | DO jj = 1,jpj |
---|
| 165 | DO ji = 1,jpi |
---|
| 166 | C |
---|
| 167 | C POC disaggregation by turbulence and bacterial activity. |
---|
| 168 | C ------------------------------------------------------------- |
---|
| 169 | C |
---|
[339] | 170 | remip=xremip*zstep*tgfunc(ji,jj,jk)*(1.-0.5*nitrfac(ji,jj,jk)) |
---|
[186] | 171 | # if defined key_off_degrad |
---|
| 172 | & *facvol(ji,jj,jk) |
---|
| 173 | # endif |
---|
| 174 | C |
---|
| 175 | C POC disaggregation rate is reduced in anoxic zone as shown by |
---|
| 176 | C sediment traps data. In oxic area, the exponent of the martin's |
---|
| 177 | C law is around -0.87. In anoxic zone, it is around -0.35. This |
---|
| 178 | C means a disaggregation constant about 0.5 the value in oxic zones |
---|
| 179 | C ----------------------------------------------------------------- |
---|
| 180 | C |
---|
| 181 | orem(ji,jj,jk)=remip*trn(ji,jj,jk,jppoc) |
---|
[617] | 182 | ofer(ji,jj,jk)=remip*trn(ji,jj,jk,jpsfe) |
---|
[772] | 183 | #if ! defined key_kriest |
---|
[186] | 184 | orem2(ji,jj,jk)=remip*trn(ji,jj,jk,jpgoc) |
---|
| 185 | ofer2(ji,jj,jk)=remip*trn(ji,jj,jk,jpbfe) |
---|
[617] | 186 | #else |
---|
| 187 | orem2(ji,jj,jk)=remip*trn(ji,jj,jk,jpnum) |
---|
| 188 | #endif |
---|
[186] | 189 | C |
---|
| 190 | END DO |
---|
| 191 | END DO |
---|
| 192 | END DO |
---|
| 193 | |
---|
| 194 | DO jk = 1,jpkm1 |
---|
| 195 | DO jj = 1,jpj |
---|
| 196 | DO ji = 1,jpi |
---|
| 197 | C |
---|
| 198 | C Remineralization rate of BSi depedant on T and saturation |
---|
| 199 | C --------------------------------------------------------- |
---|
| 200 | C |
---|
| 201 | zsatur=(sio3eq(ji,jj,jk)-trn(ji,jj,jk,jpsil))/ |
---|
| 202 | & (sio3eq(ji,jj,jk)+rtrn) |
---|
[339] | 203 | zsatur=max(rtrn,zsatur) |
---|
| 204 | zsatur2=zsatur*(1.+tn(ji,jj,jk)/400.)**4 |
---|
| 205 | znusil=0.225*(1.+tn(ji,jj,jk)/15.)*zsatur+0.775*zsatur2**9 |
---|
[260] | 206 | |
---|
[339] | 207 | siremin=xsirem*zstep*znusil |
---|
[186] | 208 | # if defined key_off_degrad |
---|
| 209 | & *facvol(ji,jj,jk) |
---|
| 210 | # endif |
---|
| 211 | C |
---|
| 212 | osil(ji,jj,jk)=siremin*trn(ji,jj,jk,jpdsi) |
---|
| 213 | END DO |
---|
| 214 | END DO |
---|
| 215 | END DO |
---|
[339] | 216 | C |
---|
| 217 | fesatur(:,:,:)=0.6E-9 |
---|
| 218 | C |
---|
[186] | 219 | DO jk = 1,jpkm1 |
---|
| 220 | DO jj = 1,jpj |
---|
| 221 | DO ji = 1,jpi |
---|
| 222 | C |
---|
| 223 | C scavenging rate of iron. this scavenging rate depends on the |
---|
| 224 | C load in particles on which they are adsorbed. The |
---|
| 225 | C parameterization has been taken from studies on Th |
---|
| 226 | C ------------------------------------------------------------ |
---|
| 227 | C |
---|
| 228 | xkeq=fekeq(ji,jj,jk) |
---|
| 229 | xfeequi=(-(1.+fesatur(ji,jj,jk)*xkeq-xkeq*trn(ji,jj,jk,jpfer))+ |
---|
| 230 | & sqrt((1.+fesatur(ji,jj,jk)*xkeq-xkeq*trn(ji,jj,jk,jpfer))**2 |
---|
| 231 | & +4.*trn(ji,jj,jk,jpfer)*xkeq))/(2.*xkeq) |
---|
| 232 | |
---|
[772] | 233 | #if ! defined key_kriest |
---|
[186] | 234 | xlam1b=3E-5+xlam1*(trn(ji,jj,jk,jppoc) |
---|
| 235 | & +trn(ji,jj,jk,jpgoc)+trn(ji,jj,jk,jpcal)+ |
---|
| 236 | & trn(ji,jj,jk,jpdsi))*1E6 |
---|
[617] | 237 | #else |
---|
| 238 | xlam1b=3E-5+xlam1*(trn(ji,jj,jk,jppoc) |
---|
| 239 | & +trn(ji,jj,jk,jpcal)+trn(ji,jj,jk,jpdsi))*1E6 |
---|
| 240 | #endif |
---|
[339] | 241 | xscave(ji,jj,jk)=xfeequi*xlam1b*zstep |
---|
[186] | 242 | # if defined key_off_degrad |
---|
| 243 | & *facvol(ji,jj,jk) |
---|
| 244 | # endif |
---|
| 245 | C |
---|
| 246 | C Increased scavenging for very high iron concentrations |
---|
| 247 | C found near the coasts due to increased lithogenic particles |
---|
| 248 | C and let's say it unknown processes (precipitation, ...) |
---|
| 249 | C ----------------------------------------------------------- |
---|
| 250 | C |
---|
[339] | 251 | zlamfac=max(0.,(gphit(ji,jj)+55.)/30.) |
---|
| 252 | zlamfac=min(1.,zlamfac) |
---|
[772] | 253 | #if ! defined key_kriest |
---|
[725] | 254 | xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698. |
---|
[339] | 255 | & *trn(ji,jj,jk,jppoc)+1.05E4*trn(ji,jj,jk,jpgoc)) |
---|
[725] | 256 | & *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0., |
---|
[339] | 257 | & (trn(ji,jj,jk,jpfer)*1E9-1.)) |
---|
[617] | 258 | #else |
---|
[725] | 259 | xlam1b=(80.*(trn(ji,jj,jk,jpdoc)+35E-6)+698. |
---|
[617] | 260 | & *trn(ji,jj,jk,jppoc)) |
---|
[725] | 261 | & *zdiss(ji,jj,jk)+1E-4*(1.-zlamfac)+xlam1*max(0., |
---|
[617] | 262 | & (trn(ji,jj,jk,jpfer)*1E9-1.)) |
---|
| 263 | #endif |
---|
[339] | 264 | |
---|
[617] | 265 | |
---|
[725] | 266 | xaggdfe(ji,jj,jk)=xlam1b*zstep*0.5*(trn(ji,jj,jk,jpfer) |
---|
[339] | 267 | & -xfeequi) |
---|
[186] | 268 | # if defined key_off_degrad |
---|
| 269 | & *facvol(ji,jj,jk) |
---|
| 270 | # endif |
---|
[339] | 271 | |
---|
[186] | 272 | C |
---|
| 273 | END DO |
---|
| 274 | END DO |
---|
| 275 | END DO |
---|
| 276 | C |
---|
| 277 | #endif |
---|
| 278 | RETURN |
---|
| 279 | END |
---|