[618] | 1 | |
---|
[719] | 2 | |
---|
[618] | 3 | CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
| 4 | CCC p4zbio : PISCES MODEL - Kriest parameterization |
---|
| 5 | CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
| 6 | |
---|
| 7 | DO jk = 1,jpkm1 |
---|
| 8 | DO jj = 1,jpj |
---|
| 9 | DO ji = 1,jpi |
---|
| 10 | C |
---|
| 11 | C Evolution of DOC |
---|
| 12 | C ---------------- |
---|
| 13 | C |
---|
| 14 | trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc) |
---|
| 15 | & +orem(ji,jj,jk)+excret2*prorca2(ji,jj,jk) |
---|
| 16 | & +excret*prorca(ji,jj,jk)-olimi(ji,jj,jk)-denitr(ji,jj,jk) |
---|
| 17 | & +grarem(ji,jj,jk)*(1.-sigma1)+grarem2(ji,jj,jk) |
---|
| 18 | & *(1.-sigma2)-xaggdoc(ji,jj,jk) |
---|
| 19 | END DO |
---|
| 20 | END DO |
---|
| 21 | END DO |
---|
| 22 | |
---|
| 23 | |
---|
| 24 | DO jk = 1,jpkm1 |
---|
| 25 | DO jj = 1,jpj |
---|
| 26 | DO ji = 1,jpi |
---|
| 27 | C |
---|
| 28 | C Evolution of Detritus |
---|
| 29 | C --------------------- |
---|
| 30 | C |
---|
| 31 | znumpoc=trn(ji,jj,jk,jpnum)/(trn(ji,jj,jk,jppoc)+rtrn) |
---|
| 32 | znumdoc=znumpoc |
---|
| 33 | |
---|
| 34 | trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc) |
---|
| 35 | & -grazpoc(ji,jj,jk) |
---|
| 36 | & +grapoc(ji,jj,jk)+grapoc2(ji,jj,jk) |
---|
| 37 | & -grazm(ji,jj,jk)+tortz2(ji,jj,jk) |
---|
| 38 | & +respz(ji,jj,jk)+respz2(ji,jj,jk) |
---|
| 39 | & +respp(ji,jj,jk)+respp2(ji,jj,jk) |
---|
| 40 | & +tortp2(ji,jj,jk)+tortz(ji,jj,jk) |
---|
| 41 | & +tortp(ji,jj,jk)-orem(ji,jj,jk) |
---|
| 42 | & +xaggdoc(ji,jj,jk)-grazffe(ji,jj,jk) |
---|
| 43 | C |
---|
| 44 | C Evolution of number of aggregates |
---|
| 45 | C --------------------------------- |
---|
| 46 | C |
---|
| 47 | trn(ji,jj,jk,jpnum) = trn(ji,jj,jk,jpnum) |
---|
| 48 | & -xagg(ji,jj,jk)-(orem(ji,jj,jk) |
---|
| 49 | & +grazpoc(ji,jj,jk))*znumpoc |
---|
| 50 | & +(tortp(ji,jj,jk))*xkr_nnano |
---|
| 51 | & +(tortp2(ji,jj,jk)+respp(ji,jj,jk)+tortz(ji,jj,jk) |
---|
| 52 | & +grapoc(ji,jj,jk)-grazm(ji,jj,jk) |
---|
| 53 | & +respz(ji,jj,jk))*xkr_ndiat |
---|
| 54 | & +(grapoc2(ji,jj,jk)+tortz2(ji,jj,jk) |
---|
| 55 | & +respz2(ji,jj,jk))*xkr_nmeso |
---|
| 56 | & +respp2(ji,jj,jk)*xkr_naggr |
---|
| 57 | & +xaggdoc(ji,jj,jk)*znumdoc |
---|
| 58 | & -grazffe(ji,jj,jk)*znumpoc*wsbio4(ji,jj,jk) |
---|
| 59 | & /(wsbio3(ji,jj,jk)+rtrn) |
---|
| 60 | C |
---|
| 61 | trn(ji,jj,jk,jpnum) = max(trn(ji,jj,jk,jpnum), |
---|
| 62 | & trn(ji,jj,jk,jppoc)/xkr_massp/xnumm(jk)) |
---|
| 63 | C |
---|
| 64 | trn(ji,jj,jk,jpnum) = min(trn(ji,jj,jk,jpnum), |
---|
| 65 | & trn(ji,jj,jk,jppoc)/xkr_massp/1.1) |
---|
| 66 | C |
---|
| 67 | |
---|
| 68 | END DO |
---|
| 69 | END DO |
---|
| 70 | END DO |
---|
| 71 | |
---|
| 72 | |
---|
| 73 | |
---|
| 74 | DO jk = 1,jpkm1 |
---|
| 75 | DO jj = 1,jpj |
---|
| 76 | DO ji = 1,jpi |
---|
| 77 | C |
---|
| 78 | C Evolution of biogenic Iron |
---|
| 79 | C -------------------------- |
---|
| 80 | C |
---|
| 81 | trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe) |
---|
| 82 | & +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk)) |
---|
| 83 | & -(1.-unass2)*grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk) |
---|
| 84 | & -(1.-unass2)*grazfff(ji,jj,jk)+unass2*(graznf(ji,jj,jk) |
---|
| 85 | & +grazf(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3 |
---|
| 86 | & *(tortz2(ji,jj,jk)+respz2(ji,jj,jk)+tortz(ji,jj,jk) |
---|
| 87 | & +respz(ji,jj,jk))-ofer(ji,jj,jk)+(respnf(ji,jj,jk) |
---|
| 88 | & +tortnf(ji,jj,jk))+tortdf(ji,jj,jk)+respdf(ji,jj,jk) |
---|
| 89 | & +xbactfer(ji,jj,jk)+xscave(ji,jj,jk)*zdenom1(ji,jj,jk) |
---|
| 90 | END DO |
---|
| 91 | END DO |
---|
| 92 | END DO |
---|
| 93 | |
---|