| 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 | SUBROUTINE trcbio(kt) |
|---|
| 7 | #if defined key_passivetrc && defined key_trc_lobster1 |
|---|
| 8 | CCC--------------------------------------------------------------------- |
|---|
| 9 | CCC |
|---|
| 10 | CCC ROUTINE trcbio |
|---|
| 11 | CCC ******************* |
|---|
| 12 | CCC |
|---|
| 13 | CCC PURPOSE : |
|---|
| 14 | CCC --------- |
|---|
| 15 | CCC compute the now trend due to biogeochemical processes |
|---|
| 16 | CCC and add it to the general trend of passive tracers equations. |
|---|
| 17 | CCC |
|---|
| 18 | CCC Three options: |
|---|
| 19 | CCC Default option : no biological trend |
|---|
| 20 | CCC IF 'key_trc_lobster1' : LOBSTER1 model |
|---|
| 21 | CCC |
|---|
| 22 | CC METHOD : |
|---|
| 23 | CC ------- |
|---|
| 24 | CC each now biological flux is calculated in FUNCTION of now |
|---|
| 25 | CC concentrations of tracers. |
|---|
| 26 | CC depending on the tracer, these fluxes are sources or sinks. |
|---|
| 27 | CC the total of the sources and sinks for each tracer |
|---|
| 28 | CC is added to the general trend. |
|---|
| 29 | CC |
|---|
| 30 | CC tra = tra + zf...tra - zftra... |
|---|
| 31 | CC | | |
|---|
| 32 | CC | | |
|---|
| 33 | CC source sink |
|---|
| 34 | CC |
|---|
| 35 | CC |
|---|
| 36 | CC IF 'key_trc_diabio' key is activated, the biogeochemical |
|---|
| 37 | CC trends for passive tracers are saved for futher diagnostics. |
|---|
| 38 | CC |
|---|
| 39 | CC multitasked on vertical slab (jj-loop) |
|---|
| 40 | CC |
|---|
| 41 | CC ----- |
|---|
| 42 | CC argument |
|---|
| 43 | CC ktask : task identificator |
|---|
| 44 | CC kt : time step |
|---|
| 45 | CC COMMON |
|---|
| 46 | CC /comcoo/ : orthogonal curvilinear coordinates |
|---|
| 47 | CC and scale factors |
|---|
| 48 | CC depths |
|---|
| 49 | CC /cottrp/ : present and next fields for passive |
|---|
| 50 | CC : tracers |
|---|
| 51 | CC /comtsk/ : multitasking |
|---|
| 52 | CC /comtke/ : emin, en() |
|---|
| 53 | CC /cotbio/ : biological parameters |
|---|
| 54 | CC |
|---|
| 55 | CC OUTPUT : |
|---|
| 56 | CC ------ |
|---|
| 57 | CC COMMON |
|---|
| 58 | CC /cottrp/ tra : general tracer trend increased by the |
|---|
| 59 | CC now horizontal tracer advection trend |
|---|
| 60 | CC /cottbd/ trbio : now horizontal tracer advection trend |
|---|
| 61 | CC (IF 'key_trc_diabio' is activated) |
|---|
| 62 | CC |
|---|
| 63 | CC WORKSPACE : |
|---|
| 64 | CC --------- |
|---|
| 65 | CC local |
|---|
| 66 | CC zdet,zzoo,zphy,znh4,zno3,zdom : now concentrations |
|---|
| 67 | CC zlt,zlno3,zlnh4,zle : limitation terms for phyto |
|---|
| 68 | CC zfno3phy and so on.. : fluxes between bio boxes |
|---|
| 69 | CC zphya,zzooa,zdeta, ... : after bio trends |
|---|
| 70 | CC zppz, zpdz, zpppz, zppdz, zfood : preferences terms |
|---|
| 71 | CC zfilpz, zfilpd : filtration terms |
|---|
| 72 | CC COMMON |
|---|
| 73 | CC |
|---|
| 74 | CC EXTERNAL : no |
|---|
| 75 | CC -------- |
|---|
| 76 | CC |
|---|
| 77 | CC REFERENCES : no |
|---|
| 78 | CC ---------- |
|---|
| 79 | CC |
|---|
| 80 | CC MODIFICATIONS: |
|---|
| 81 | CC -------------- |
|---|
| 82 | CC original : 99-07 (M. Levy) |
|---|
| 83 | CC 00-12 (E. Kestenare): assign a parameter |
|---|
| 84 | CC to name individual tracers |
|---|
| 85 | CC 01-03 (M. Levy) LNO3 + dia2d |
|---|
| 86 | CC---------------------------------------------------------------------- |
|---|
| 87 | CC---------------------------------------------------------------------- |
|---|
| 88 | USE oce_trc |
|---|
| 89 | USE trp_trc |
|---|
| 90 | USE sms |
|---|
| 91 | USE lbclnk |
|---|
| 92 | IMPLICIT NONE |
|---|
| 93 | CC local declarations |
|---|
| 94 | CC ================== |
|---|
| 95 | INTEGER kt |
|---|
| 96 | INTEGER ji,jj,jk,jn |
|---|
| 97 | REAL ztot(jpi), ze3t(jpk) |
|---|
| 98 | REAL zdet,zzoo,zphy,zno3,znh4,zdom,zlno3,zlnh4,zle,zlt |
|---|
| 99 | REAL zno3phy, znh4phy, zphynh4, zphydom, zphydet, zphyzoo, zdetzoo |
|---|
| 100 | $ ,zzoonh4, zzoodom, zzoodet, zdetnh4, zdetdom, znh4no3, zdomnh4 |
|---|
| 101 | $ ,zppz,zpdz,zpppz,zppdz,zfood,zfilpz,zfildz,zphya,zzooa,zno3a |
|---|
| 102 | $ ,znh4a,zdeta,zdoma, ztra, zzoobod, zboddet, zdomaju |
|---|
| 103 | |
|---|
| 104 | CC---------------------------------------------------------------------- |
|---|
| 105 | CC statement functions |
|---|
| 106 | CC =================== |
|---|
| 107 | CDIR$ NOLIST |
|---|
| 108 | #include "domzgr_substitute.h90" |
|---|
| 109 | CDIR$ LIST |
|---|
| 110 | CCC--------------------------------------------------------------------- |
|---|
| 111 | CCC OPA8, LODYC (07/99) |
|---|
| 112 | CCC--------------------------------------------------------------------- |
|---|
| 113 | C | --------------| |
|---|
| 114 | C | LOBSTER1 MODEL| |
|---|
| 115 | C | --------------| |
|---|
| 116 | |
|---|
| 117 | #if defined key_trc_diaadd |
|---|
| 118 | C convert fluxes in per day |
|---|
| 119 | DO jk=1,jpkbm1 |
|---|
| 120 | ze3t(jk)=e3t(jk)*86400. |
|---|
| 121 | END DO |
|---|
| 122 | DO jk=jpkb,jpk |
|---|
| 123 | ze3t(jk)=0. |
|---|
| 124 | END DO |
|---|
| 125 | #endif |
|---|
| 126 | C |
|---|
| 127 | C vertical slab |
|---|
| 128 | C ============= |
|---|
| 129 | C |
|---|
| 130 | DO 1000 jj = 2,jpjm1 |
|---|
| 131 | C |
|---|
| 132 | C 1. biological level |
|---|
| 133 | C =================== |
|---|
| 134 | C |
|---|
| 135 | DO ji = 2,jpim1 |
|---|
| 136 | fbod(ji,jj)=0. |
|---|
| 137 | #if defined key_trc_diaadd |
|---|
| 138 | DO jn=1,jpdia2d |
|---|
| 139 | trc2d(ji,jj,jn)=0. |
|---|
| 140 | END DO |
|---|
| 141 | #endif |
|---|
| 142 | END DO |
|---|
| 143 | |
|---|
| 144 | DO jk=1,jpkbm1 |
|---|
| 145 | DO ji = 2,jpim1 |
|---|
| 146 | C |
|---|
| 147 | C |
|---|
| 148 | C 1.1 trophic variables( det, zoo, phy, no3, nh4, dom) |
|---|
| 149 | C --------------------------------------------------- |
|---|
| 150 | C |
|---|
| 151 | C negative trophic variables DO not contribute to the fluxes |
|---|
| 152 | C |
|---|
| 153 | zdet = max(0.,trn(ji,jj,jk,jpdet)) |
|---|
| 154 | zzoo = max(0.,trn(ji,jj,jk,jpzoo)) |
|---|
| 155 | zphy = max(0.,trn(ji,jj,jk,jpphy)) |
|---|
| 156 | zno3 = max(0.,trn(ji,jj,jk,jpno3)) |
|---|
| 157 | znh4 = max(0.,trn(ji,jj,jk,jpnh4)) |
|---|
| 158 | zdom = max(0.,trn(ji,jj,jk,jpdom)) |
|---|
| 159 | C |
|---|
| 160 | C |
|---|
| 161 | C 1.2 Limitations |
|---|
| 162 | C ---------------- |
|---|
| 163 | C |
|---|
| 164 | zlt = 1. |
|---|
| 165 | zle = 1. - exp( -xpar(ji,jj,jk)/aki/zlt) |
|---|
| 166 | C psinut,akno3,aknh4 added by asklod AS Kremeur 2005-03 |
|---|
| 167 | zlno3 = zno3* exp(-psinut*znh4) / (akno3+zno3) |
|---|
| 168 | zlnh4 = znh4 / (znh4+aknh4) |
|---|
| 169 | |
|---|
| 170 | C |
|---|
| 171 | C |
|---|
| 172 | C 1.3 sinks and sources |
|---|
| 173 | C --------------------- |
|---|
| 174 | C |
|---|
| 175 | C |
|---|
| 176 | C 1. phytoplankton production and exsudation |
|---|
| 177 | C |
|---|
| 178 | zno3phy = tmumax * zle * zlt * zlno3 * zphy |
|---|
| 179 | znh4phy = tmumax * zle * zlt * zlnh4 * zphy |
|---|
| 180 | |
|---|
| 181 | C fphylab added by asklod AS Kremeur 2005-03 |
|---|
| 182 | zphydom = rgamma * (1 - fphylab) * (zno3phy + znh4phy) |
|---|
| 183 | zphynh4 = rgamma * fphylab * (zno3phy + znh4phy) |
|---|
| 184 | |
|---|
| 185 | C |
|---|
| 186 | C 2. zooplankton production |
|---|
| 187 | C |
|---|
| 188 | C preferences |
|---|
| 189 | C |
|---|
| 190 | zppz = rppz |
|---|
| 191 | zpdz = 1. - rppz |
|---|
| 192 | zpppz = ( zppz * zphy ) / |
|---|
| 193 | $ ( ( zppz * zphy + zpdz * zdet ) + 1.e-13 ) |
|---|
| 194 | zppdz = ( zpdz * zdet ) / |
|---|
| 195 | $ ( ( zppz * zphy + zpdz * zdet ) + 1.e-13 ) |
|---|
| 196 | zfood = zpppz * zphy + zppdz * zdet |
|---|
| 197 | C |
|---|
| 198 | C filtration |
|---|
| 199 | C |
|---|
| 200 | zfilpz = taus * zpppz / (aks + zfood) |
|---|
| 201 | zfildz = taus * zppdz / (aks + zfood) |
|---|
| 202 | C |
|---|
| 203 | C grazing |
|---|
| 204 | C |
|---|
| 205 | zphyzoo = zfilpz * zphy * zzoo |
|---|
| 206 | zdetzoo = zfildz * zdet * zzoo |
|---|
| 207 | C |
|---|
| 208 | C 3. fecal pellets production |
|---|
| 209 | C |
|---|
| 210 | zzoodet = rpnaz * zphyzoo + rdnaz * zdetzoo |
|---|
| 211 | C |
|---|
| 212 | C 4. zooplankton liquide excretion |
|---|
| 213 | C |
|---|
| 214 | zzoonh4 = tauzn * fzoolab * zzoo |
|---|
| 215 | zzoodom = tauzn * (1 - fzoolab) * zzoo |
|---|
| 216 | C |
|---|
| 217 | C 5. mortality |
|---|
| 218 | C |
|---|
| 219 | C phytoplankton mortality |
|---|
| 220 | C |
|---|
| 221 | zphydet = tmminp * zphy |
|---|
| 222 | C |
|---|
| 223 | C |
|---|
| 224 | C zooplankton mortality |
|---|
| 225 | c closure : flux fbod is redistributed below level jpkbio |
|---|
| 226 | C |
|---|
| 227 | zzoobod = tmminz * zzoo * zzoo |
|---|
| 228 | fbod(ji,jj) = fbod(ji,jj) |
|---|
| 229 | $ + (1-fdbod) * zzoobod * fse3t(ji,jj,jk) |
|---|
| 230 | zboddet = fdbod * zzoobod |
|---|
| 231 | C |
|---|
| 232 | C |
|---|
| 233 | C 6. detritus and dom breakdown |
|---|
| 234 | C |
|---|
| 235 | C |
|---|
| 236 | zdetnh4 = taudn * fdetlab * zdet |
|---|
| 237 | zdetdom = taudn * (1 - fdetlab) * zdet |
|---|
| 238 | |
|---|
| 239 | zdomnh4 = taudomn * zdom |
|---|
| 240 | C |
|---|
| 241 | C flux added to express how the excess of nitrogen from |
|---|
| 242 | C PHY, ZOO and DET to DOM goes directly to NH4 (flux of ajustment) |
|---|
| 243 | zdomaju = (1 - redf/reddom) * (zphydom + zzoodom + zdetdom) |
|---|
| 244 | C |
|---|
| 245 | C 7. Nitrification |
|---|
| 246 | C |
|---|
| 247 | znh4no3 = taunn * znh4 |
|---|
| 248 | C |
|---|
| 249 | C |
|---|
| 250 | C |
|---|
| 251 | C 1.4 determination of trends |
|---|
| 252 | C --------------------------- |
|---|
| 253 | C |
|---|
| 254 | C total trend for each biological tracer |
|---|
| 255 | C |
|---|
| 256 | zphya = zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo |
|---|
| 257 | $ - zphydet |
|---|
| 258 | zzooa = zphyzoo + zdetzoo - zzoodet - zzoodom - zzoonh4 |
|---|
| 259 | $ - zzoobod |
|---|
| 260 | zno3a = - zno3phy + znh4no3 |
|---|
| 261 | znh4a = - znh4phy - znh4no3 + zphynh4 + zzoonh4 + zdomnh4 |
|---|
| 262 | $ + zdetnh4 + zdomaju |
|---|
| 263 | zdeta = zphydet + zzoodet - zdetzoo - zdetnh4 - zdetdom + |
|---|
| 264 | $ zboddet |
|---|
| 265 | zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju |
|---|
| 266 | C |
|---|
| 267 | #if defined key_trc_diabio |
|---|
| 268 | trbio(ji,jj,jk,1) = zno3phy |
|---|
| 269 | trbio(ji,jj,jk,2) = znh4phy |
|---|
| 270 | trbio(ji,jj,jk,3) = zphynh4 |
|---|
| 271 | trbio(ji,jj,jk,4) = zphydom |
|---|
| 272 | trbio(ji,jj,jk,5) = zphyzoo |
|---|
| 273 | trbio(ji,jj,jk,6) = zphydet |
|---|
| 274 | trbio(ji,jj,jk,7) = zdetzoo |
|---|
| 275 | trbio(ji,jj,jk,9) = zzoodet |
|---|
| 276 | trbio(ji,jj,jk,10) = zzoobod |
|---|
| 277 | trbio(ji,jj,jk,11) = zzoonh4 |
|---|
| 278 | trbio(ji,jj,jk,12) = zzoodom |
|---|
| 279 | trbio(ji,jj,jk,13) = znh4no3 |
|---|
| 280 | trbio(ji,jj,jk,14) = zdomnh4 |
|---|
| 281 | trbio(ji,jj,jk,15) = zdetnh4 |
|---|
| 282 | #endif |
|---|
| 283 | #if defined key_trc_diaadd |
|---|
| 284 | trc2d(ji,jj,1)=trc2d(ji,jj,1)+zno3phy*ze3t(jk) |
|---|
| 285 | trc2d(ji,jj,2)=trc2d(ji,jj,2)+znh4phy*ze3t(jk) |
|---|
| 286 | trc2d(ji,jj,3)=trc2d(ji,jj,3)+zphydom*ze3t(jk) |
|---|
| 287 | trc2d(ji,jj,4)=trc2d(ji,jj,4)+zphynh4*ze3t(jk) |
|---|
| 288 | trc2d(ji,jj,5)=trc2d(ji,jj,5)+zphyzoo*ze3t(jk) |
|---|
| 289 | trc2d(ji,jj,6)=trc2d(ji,jj,6)+zphydet*ze3t(jk) |
|---|
| 290 | trc2d(ji,jj,7)=trc2d(ji,jj,7)+zdetzoo*ze3t(jk) |
|---|
| 291 | c trend number 8 is in trcsed.F |
|---|
| 292 | trc2d(ji,jj,9)=trc2d(ji,jj,9)+zzoodet*ze3t(jk) |
|---|
| 293 | trc2d(ji,jj,10)=trc2d(ji,jj,10)+zzoobod*ze3t(jk) |
|---|
| 294 | trc2d(ji,jj,11)=trc2d(ji,jj,11)+zzoonh4*ze3t(jk) |
|---|
| 295 | trc2d(ji,jj,12)=trc2d(ji,jj,12)+zzoodom*ze3t(jk) |
|---|
| 296 | trc2d(ji,jj,13)=trc2d(ji,jj,13)+znh4no3*ze3t(jk) |
|---|
| 297 | trc2d(ji,jj,14)=trc2d(ji,jj,14)+zdomnh4*ze3t(jk) |
|---|
| 298 | trc2d(ji,jj,15)=trc2d(ji,jj,15)+zdetnh4*ze3t(jk) |
|---|
| 299 | |
|---|
| 300 | trc2d(ji,jj,16)=trc2d(ji,jj,16)+(zno3phy+znh4phy-zphynh4 |
|---|
| 301 | $ -zphydom-zphyzoo-zphydet)*ze3t(jk) |
|---|
| 302 | trc2d(ji,jj,17)=trc2d(ji,jj,17)+(zphyzoo+zdetzoo-zzoodet |
|---|
| 303 | $ -zzoobod-zzoonh4-zzoodom) *ze3t(jk) |
|---|
| 304 | trc2d(ji,jj,18)=trc2d(ji,jj,18)+zdetdom*ze3t(jk) |
|---|
| 305 | c trend number 19 is in trcexp.F |
|---|
| 306 | trc3d(ji,jj,jk,1)= zno3phy *86400 |
|---|
| 307 | trc3d(ji,jj,jk,2)= znh4phy *86400 |
|---|
| 308 | trc3d(ji,jj,jk,3)= znh4no3 *86400 |
|---|
| 309 | #endif |
|---|
| 310 | C |
|---|
| 311 | C tracer flux at totox-point added to the general trend |
|---|
| 312 | C |
|---|
| 313 | tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta |
|---|
| 314 | tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa |
|---|
| 315 | tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya |
|---|
| 316 | tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a |
|---|
| 317 | tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a |
|---|
| 318 | tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma |
|---|
| 319 | C |
|---|
| 320 | END DO |
|---|
| 321 | END DO |
|---|
| 322 | C |
|---|
| 323 | C 2. under biological level |
|---|
| 324 | C ========================= |
|---|
| 325 | C |
|---|
| 326 | DO jk = jpkb,jpk |
|---|
| 327 | C |
|---|
| 328 | C 2.1 compute the remineralisation of all quantities towards nitrate |
|---|
| 329 | C ------------------------------------------------------------------ |
|---|
| 330 | C |
|---|
| 331 | DO ji = 2,jpim1 |
|---|
| 332 | C |
|---|
| 333 | C 2.1.1 trophic variables( det, zoo, phy, no3, nh4, dom) |
|---|
| 334 | C ----------------------------------------------------- |
|---|
| 335 | C |
|---|
| 336 | C negative trophic variables DO not contribute to the fluxes |
|---|
| 337 | C |
|---|
| 338 | zdet = max(0.,trn(ji,jj,jk,jpdet)) |
|---|
| 339 | zzoo = max(0.,trn(ji,jj,jk,jpzoo)) |
|---|
| 340 | zphy = max(0.,trn(ji,jj,jk,jpphy)) |
|---|
| 341 | zno3 = max(0.,trn(ji,jj,jk,jpno3)) |
|---|
| 342 | znh4 = max(0.,trn(ji,jj,jk,jpnh4)) |
|---|
| 343 | zdom = max(0.,trn(ji,jj,jk,jpdom)) |
|---|
| 344 | CC |
|---|
| 345 | CC 2.1.2 Limitations |
|---|
| 346 | CC ---------------- |
|---|
| 347 | CC |
|---|
| 348 | zlt = 0. |
|---|
| 349 | zle = 0. |
|---|
| 350 | zlno3 = 0. |
|---|
| 351 | zlnh4 = 0. |
|---|
| 352 | CC |
|---|
| 353 | CC |
|---|
| 354 | CC 2.1.3 sinks and sources |
|---|
| 355 | CC --------------------- |
|---|
| 356 | CC |
|---|
| 357 | CC |
|---|
| 358 | CC 1. phytoplankton production and exsudation |
|---|
| 359 | CC |
|---|
| 360 | zno3phy = 0. |
|---|
| 361 | znh4phy = 0. |
|---|
| 362 | C |
|---|
| 363 | zphydom = 0. |
|---|
| 364 | zphynh4 = 0. |
|---|
| 365 | CC |
|---|
| 366 | CC 2. zooplankton production |
|---|
| 367 | CC |
|---|
| 368 | CC grazing |
|---|
| 369 | CC |
|---|
| 370 | zphyzoo = 0. |
|---|
| 371 | zdetzoo = 0. |
|---|
| 372 | CC |
|---|
| 373 | CC 3. fecal pellets production |
|---|
| 374 | CC |
|---|
| 375 | zzoodet = 0. |
|---|
| 376 | CC |
|---|
| 377 | CC 4. zooplankton liquide excretion |
|---|
| 378 | CC |
|---|
| 379 | zzoonh4 = tauzn * fzoolab * zzoo |
|---|
| 380 | zzoodom = tauzn * (1 - fzoolab) * zzoo |
|---|
| 381 | CC |
|---|
| 382 | CC 5. mortality |
|---|
| 383 | CC |
|---|
| 384 | CC phytoplankton mortality |
|---|
| 385 | CC |
|---|
| 386 | zphydet = tmminp * zphy |
|---|
| 387 | CC |
|---|
| 388 | CC |
|---|
| 389 | CC zooplankton mortality |
|---|
| 390 | Cc closure : flux fbod is redistributed below level jpkbio |
|---|
| 391 | CC |
|---|
| 392 | zzoobod = 0. |
|---|
| 393 | zboddet = 0. |
|---|
| 394 | CC |
|---|
| 395 | CC |
|---|
| 396 | CC 6. detritus and dom breakdown |
|---|
| 397 | CC |
|---|
| 398 | zdetnh4 = taudn * fdetlab * zdet |
|---|
| 399 | zdetdom = taudn * (1 - fdetlab) * zdet |
|---|
| 400 | C |
|---|
| 401 | zdomnh4 = taudomn * zdom |
|---|
| 402 | zdomaju = (1 - redf/reddom) * (zphydom + zzoodom + zdetdom) |
|---|
| 403 | CC |
|---|
| 404 | CC 7. Nitrification |
|---|
| 405 | CC |
|---|
| 406 | znh4no3 = taunn * znh4 |
|---|
| 407 | CC |
|---|
| 408 | CC |
|---|
| 409 | CC 2.1.4 determination of trends |
|---|
| 410 | CC --------------------------- |
|---|
| 411 | CC |
|---|
| 412 | CC total trend for each biological tracer |
|---|
| 413 | CC |
|---|
| 414 | zphya = zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo |
|---|
| 415 | $ - zphydet |
|---|
| 416 | zzooa = zphyzoo + zdetzoo - zzoodet - zzoodom - zzoonh4 |
|---|
| 417 | $ - zzoobod |
|---|
| 418 | zno3a = - zno3phy + znh4no3 |
|---|
| 419 | znh4a = - znh4phy - znh4no3 + zphynh4 + zzoonh4 + zdomnh4 |
|---|
| 420 | $ + zdetnh4 + zdomaju |
|---|
| 421 | zdeta = zphydet + zzoodet - zdetzoo - zdetnh4 - zdetdom + |
|---|
| 422 | $ zboddet |
|---|
| 423 | zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju |
|---|
| 424 | CC |
|---|
| 425 | #if defined key_trc_diabio |
|---|
| 426 | trbio(ji,jj,jk,1) = zno3phy |
|---|
| 427 | trbio(ji,jj,jk,2) = znh4phy |
|---|
| 428 | trbio(ji,jj,jk,3) = zphynh4 |
|---|
| 429 | trbio(ji,jj,jk,4) = zphydom |
|---|
| 430 | trbio(ji,jj,jk,5) = zphyzoo |
|---|
| 431 | trbio(ji,jj,jk,6) = zphydet |
|---|
| 432 | trbio(ji,jj,jk,7) = zdetzoo |
|---|
| 433 | trbio(ji,jj,jk,9) = zzoodet |
|---|
| 434 | trbio(ji,jj,jk,10) = zzoobod |
|---|
| 435 | trbio(ji,jj,jk,11) = zzoonh4 |
|---|
| 436 | trbio(ji,jj,jk,12) = zzoodom |
|---|
| 437 | trbio(ji,jj,jk,13) = znh4no3 |
|---|
| 438 | trbio(ji,jj,jk,14) = zdomnh4 |
|---|
| 439 | trbio(ji,jj,jk,15) = zdetnh4 |
|---|
| 440 | #endif |
|---|
| 441 | #if defined key_trc_diaadd |
|---|
| 442 | trc2d(ji,jj,1)=trc2d(ji,jj,1)+zno3phy*ze3t(jk) |
|---|
| 443 | trc2d(ji,jj,2)=trc2d(ji,jj,2)+znh4phy*ze3t(jk) |
|---|
| 444 | trc2d(ji,jj,3)=trc2d(ji,jj,3)+zphydom*ze3t(jk) |
|---|
| 445 | trc2d(ji,jj,4)=trc2d(ji,jj,4)+zphynh4*ze3t(jk) |
|---|
| 446 | trc2d(ji,jj,5)=trc2d(ji,jj,5)+zphyzoo*ze3t(jk) |
|---|
| 447 | trc2d(ji,jj,6)=trc2d(ji,jj,6)+zphydet*ze3t(jk) |
|---|
| 448 | trc2d(ji,jj,7)=trc2d(ji,jj,7)+zdetzoo*ze3t(jk) |
|---|
| 449 | Cc trend number 8 is in trcsed.F |
|---|
| 450 | trc2d(ji,jj,9)=trc2d(ji,jj,9)+zzoodet*ze3t(jk) |
|---|
| 451 | trc2d(ji,jj,10)=trc2d(ji,jj,10)+zzoobod*ze3t(jk) |
|---|
| 452 | trc2d(ji,jj,11)=trc2d(ji,jj,11)+zzoonh4*ze3t(jk) |
|---|
| 453 | trc2d(ji,jj,12)=trc2d(ji,jj,12)+zzoodom*ze3t(jk) |
|---|
| 454 | trc2d(ji,jj,13)=trc2d(ji,jj,13)+znh4no3*ze3t(jk) |
|---|
| 455 | trc2d(ji,jj,14)=trc2d(ji,jj,14)+zdomnh4*ze3t(jk) |
|---|
| 456 | trc2d(ji,jj,15)=trc2d(ji,jj,15)+zdetnh4*ze3t(jk) |
|---|
| 457 | |
|---|
| 458 | trc2d(ji,jj,16)=trc2d(ji,jj,16)+(zno3phy+znh4phy-zphynh4 |
|---|
| 459 | $ -zphydom-zphyzoo-zphydet)*ze3t(jk) |
|---|
| 460 | trc2d(ji,jj,17)=trc2d(ji,jj,17)+(zphyzoo+zdetzoo-zzoodet |
|---|
| 461 | $ -zzoobod-zzoonh4-zzoodom) *ze3t(jk) |
|---|
| 462 | trc2d(ji,jj,18)=trc2d(ji,jj,18)+zdetdom*ze3t(jk) |
|---|
| 463 | |
|---|
| 464 | trc3d(ji,jj,jk,1)= zno3phy *86400 |
|---|
| 465 | trc3d(ji,jj,jk,2)= znh4phy *86400 |
|---|
| 466 | trc3d(ji,jj,jk,3)= znh4no3 *86400 |
|---|
| 467 | #endif |
|---|
| 468 | CC |
|---|
| 469 | CC tracer flux at totox-point added to the general trend |
|---|
| 470 | CC |
|---|
| 471 | tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta |
|---|
| 472 | tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa |
|---|
| 473 | tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya |
|---|
| 474 | tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a |
|---|
| 475 | tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a |
|---|
| 476 | tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma |
|---|
| 477 | CC |
|---|
| 478 | END DO |
|---|
| 479 | END DO |
|---|
| 480 | |
|---|
| 481 | |
|---|
| 482 | |
|---|
| 483 | |
|---|
| 484 | c$$$ DO jk = jpkb,jpk |
|---|
| 485 | c$$$C |
|---|
| 486 | c$$$C 2.1 Old way to compute the remineralisation : asklod AS Kremeur (before 2005-03) |
|---|
| 487 | c$$$C ------------------------------------------------------------------ |
|---|
| 488 | c$$$C |
|---|
| 489 | c$$$ DO ji=2,jpim1 |
|---|
| 490 | c$$$ ztot(ji) = 0. |
|---|
| 491 | c$$$ END DO |
|---|
| 492 | c$$$ DO jn=1,jptra |
|---|
| 493 | c$$$ IF (ctrcnm(jn).NE.'NO3') THEN |
|---|
| 494 | c$$$ DO ji=2,jpim1 |
|---|
| 495 | c$$$ ztra = remdmp(jk,jn) * trn(ji,jj,jk,jn) |
|---|
| 496 | c$$$ tra(ji,jj,jk,jn) = tra(ji,jj,jk,jn) - ztra |
|---|
| 497 | c$$$ ztot(ji) = ztot(ji) + ztra |
|---|
| 498 | c$$$ END DO |
|---|
| 499 | c$$$ ENDIF |
|---|
| 500 | c$$$ END DO |
|---|
| 501 | c$$$ DO jn=1,jptra |
|---|
| 502 | c$$$ IF (ctrcnm(jn).EQ.'NO3') THEN |
|---|
| 503 | c$$$ DO ji=2,jpim1 |
|---|
| 504 | c$$$ tra(ji,jj,jk,jn) = tra(ji,jj,jk,jn) + ztot(ji) |
|---|
| 505 | c$$$ END DO |
|---|
| 506 | c$$$#if defined key_trc_diabio |
|---|
| 507 | c$$$ trbio(ji,jj,jk,1)=ztot(ji) |
|---|
| 508 | c$$$#endif |
|---|
| 509 | c$$$ ENDIF |
|---|
| 510 | c$$$ END DO |
|---|
| 511 | c$$$ END DO |
|---|
| 512 | |
|---|
| 513 | C |
|---|
| 514 | C |
|---|
| 515 | C END of slab |
|---|
| 516 | C =========== |
|---|
| 517 | C |
|---|
| 518 | 1000 CONTINUE |
|---|
| 519 | |
|---|
| 520 | #if defined key_trc_diaadd |
|---|
| 521 | |
|---|
| 522 | C Lateral boundary conditions on trc2d |
|---|
| 523 | DO jn=1,jpdia2d |
|---|
| 524 | CALL lbc_lnk(trc2d(:,:,jn),'T',1. ) |
|---|
| 525 | END DO |
|---|
| 526 | |
|---|
| 527 | C Lateral boundary conditions on trc3d |
|---|
| 528 | DO jn=1,jpdia3d |
|---|
| 529 | CALL lbc_lnk(trc3d(:,:,1,jn),'T',1. ) |
|---|
| 530 | END DO |
|---|
| 531 | |
|---|
| 532 | #endif |
|---|
| 533 | |
|---|
| 534 | #if defined key_trc_diabio |
|---|
| 535 | C Lateral boundary conditions on trcbio |
|---|
| 536 | DO jn=1,jpdiabio |
|---|
| 537 | CALL lbc_lnk(trbio(:,:,1,jn),'T',1. ) |
|---|
| 538 | END DO |
|---|
| 539 | #endif |
|---|
| 540 | |
|---|
| 541 | # else |
|---|
| 542 | C |
|---|
| 543 | C no biological model |
|---|
| 544 | C |
|---|
| 545 | # endif |
|---|
| 546 | |
|---|
| 547 | C |
|---|
| 548 | C |
|---|
| 549 | RETURN |
|---|
| 550 | END |
|---|
