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