[699] | 1 | CC---------------------------------------------------------------------- |
---|
| 2 | CC TOP 1.0 , LOCEAN-IPSL (2005) |
---|
| 3 | CC $Id$ |
---|
| 4 | CC This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt |
---|
| 5 | CC---------------------------------------------------------------------- |
---|
| 6 | |
---|
[618] | 7 | CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
| 8 | CCC p4zbio : PISCES MODEL - Standard parameterization |
---|
| 9 | CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
| 10 | |
---|
| 11 | DO jk = 1,jpkm1 |
---|
| 12 | DO jj = 1,jpj |
---|
| 13 | DO ji = 1,jpi |
---|
| 14 | C |
---|
| 15 | C Evolution of DOC |
---|
| 16 | C ---------------- |
---|
| 17 | C |
---|
| 18 | trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc) |
---|
| 19 | & +orem(ji,jj,jk)+excret2*prorca2(ji,jj,jk) |
---|
| 20 | & +excret*prorca(ji,jj,jk)-olimi(ji,jj,jk)-denitr(ji,jj,jk) |
---|
| 21 | & +grarem(ji,jj,jk)*(1.-sigma1)+grarem2(ji,jj,jk) |
---|
| 22 | & *(1.-sigma2)-xaggdoc(ji,jj,jk)-xaggdoc2(ji,jj,jk) |
---|
| 23 | END DO |
---|
| 24 | END DO |
---|
| 25 | END DO |
---|
| 26 | |
---|
| 27 | DO jk = 1,jpkm1 |
---|
| 28 | DO jj = 1,jpj |
---|
| 29 | DO ji = 1,jpi |
---|
| 30 | C |
---|
| 31 | C Evolution of Detritus |
---|
| 32 | C --------------------- |
---|
| 33 | C |
---|
| 34 | trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc) |
---|
| 35 | & -grazpoc(ji,jj,jk)+grapoc(ji,jj,jk)-grazm(ji,jj,jk) |
---|
| 36 | & +respz(ji,jj,jk)-xagg(ji,jj,jk)+xaggdoc(ji,jj,jk) |
---|
| 37 | & +(1.-0.5*zfracal(ji,jj,jk))*(tortp(ji,jj,jk) |
---|
| 38 | & +respp(ji,jj,jk))+0.5*tortp2(ji,jj,jk) |
---|
| 39 | & +orem2(ji,jj,jk)+tortz(ji,jj,jk)-orem(ji,jj,jk) |
---|
| 40 | C |
---|
| 41 | C Evolution of rapid Detritus |
---|
| 42 | C --------------------- |
---|
| 43 | C |
---|
| 44 | trn(ji,jj,jk,jpgoc) = trn(ji,jj,jk,jpgoc) |
---|
| 45 | & +grapoc2(ji,jj,jk)+respp2(ji,jj,jk)+xagg(ji,jj,jk) |
---|
| 46 | & +tortz2(ji,jj,jk)+respz2(ji,jj,jk)-orem2(ji,jj,jk) |
---|
| 47 | & +0.5*zfracal(ji,jj,jk)*(respp(ji,jj,jk)+tortp(ji,jj,jk)) |
---|
| 48 | & +0.5*tortp2(ji,jj,jk)+xaggdoc2(ji,jj,jk)-grazffe(ji,jj,jk) |
---|
| 49 | C |
---|
| 50 | END DO |
---|
| 51 | END DO |
---|
| 52 | END DO |
---|
| 53 | |
---|
| 54 | |
---|
| 55 | |
---|
| 56 | DO jk = 1,jpkm1 |
---|
| 57 | DO jj = 1,jpj |
---|
| 58 | DO ji = 1,jpi |
---|
| 59 | C |
---|
| 60 | C Evolution of small biogenic Iron |
---|
| 61 | C -------------------------- |
---|
| 62 | C |
---|
| 63 | trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe) |
---|
| 64 | & +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk)) |
---|
| 65 | & -grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk) |
---|
| 66 | & +(1.-0.5*zfracal(ji,jj,jk))*(tortnf(ji,jj,jk) |
---|
| 67 | & +respnf(ji,jj,jk))+0.5*tortdf(ji,jj,jk)+ferat3* |
---|
| 68 | & (tortz(ji,jj,jk)+respz(ji,jj,jk))-ofer(ji,jj,jk) |
---|
| 69 | & +ofer2(ji,jj,jk)-xaggfe(ji,jj,jk) |
---|
| 70 | & +xscave(ji,jj,jk)*zdenom1(ji,jj,jk) |
---|
| 71 | C |
---|
| 72 | C Evolution of big biogenic Iron |
---|
| 73 | C -------------------------- |
---|
| 74 | C |
---|
| 75 | trn(ji,jj,jk,jpbfe) = trn(ji,jj,jk,jpbfe) |
---|
| 76 | & +unass2*(graznf(ji,jj,jk)+grazf(ji,jj,jk)+grazfff(ji,jj,jk) |
---|
| 77 | & +grazpof(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3* |
---|
| 78 | & (tortz2(ji,jj,jk)+respz2(ji,jj,jk))-ofer2(ji,jj,jk) |
---|
| 79 | & +0.5*zfracal(ji,jj,jk)*(respnf(ji,jj,jk)+tortnf(ji,jj,jk)) |
---|
| 80 | & +0.5*tortdf(ji,jj,jk)+respdf(ji,jj,jk)+xaggfe(ji,jj,jk) |
---|
| 81 | & +xbactfer(ji,jj,jk)-grazfff(ji,jj,jk)+xscave(ji,jj,jk) |
---|
| 82 | & *zdenom2(ji,jj,jk) |
---|
| 83 | END DO |
---|
| 84 | END DO |
---|
| 85 | END DO |
---|
| 86 | |
---|