| 1 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
|---|
| 2 | !CC p4zbio : PISCES MODEL - Standard parameterization |
|---|
| 3 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
|---|
| 4 | |
|---|
| 5 | DO jk = 1, jpkm1 |
|---|
| 6 | |
|---|
| 7 | ! Evolution of DOC |
|---|
| 8 | ! ---------------- |
|---|
| 9 | trn(:,:,jk,jpdoc) = trn(:,:,jk,jpdoc) + orem(:,:,jk) + excret2 * prorca2(:,:,jk) & |
|---|
| 10 | & + excret * prorca(:,:,jk) - olimi(:,:,jk) - denitr(:,:,jk) & |
|---|
| 11 | & + grarem(:,:,jk) * (1.-sigma1) + grarem2(:,:,jk) * (1.-sigma2) & |
|---|
| 12 | & - xaggdoc(:,:,jk) - xaggdoc2(:,:,jk) |
|---|
| 13 | |
|---|
| 14 | ! Evolution of Detritus |
|---|
| 15 | ! --------------------- |
|---|
| 16 | trn(:,:,jk,jppoc) = trn(:,:,jk,jppoc) - grazpoc(:,:,jk) + grapoc(:,:,jk) - grazm (:,:,jk) & |
|---|
| 17 | & + respz (:,:,jk) - xagg (:,:,jk) + xaggdoc(:,:,jk) & |
|---|
| 18 | & + ( 1.-0.5 * zfracal(:,:,jk) ) * ( tortp(:,:,jk) + respp(:,:,jk) ) & |
|---|
| 19 | & + 0.5 * tortp2 (:,:,jk) & |
|---|
| 20 | & + orem2(:,:,jk) + tortz(:,:,jk) - orem(:,:,jk) |
|---|
| 21 | |
|---|
| 22 | ! Evolution of rapid Detritus |
|---|
| 23 | ! --------------------------- |
|---|
| 24 | trn(:,:,jk,jpgoc) = trn(:,:,jk,jpgoc) + grapoc2(:,:,jk) + respp2(:,:,jk) + xagg (:,:,jk) & |
|---|
| 25 | & + tortz2 (:,:,jk) + respz2(:,:,jk) - orem2(:,:,jk) & |
|---|
| 26 | & + 0.5*zfracal(:,:,jk)*(respp(:,:,jk)+tortp(:,:,jk)) & |
|---|
| 27 | & + 0.5*tortp2(:,:,jk)+xaggdoc2(:,:,jk)-grazffe(:,:,jk) |
|---|
| 28 | |
|---|
| 29 | ! Evolution of small biogenic Iron |
|---|
| 30 | ! -------------------------------- |
|---|
| 31 | trn(:,:,jk,jpsfe) = trn(:,:,jk,jpsfe) + unass * ( grazpf(:,:,jk) + grazsf(:,:,jk) ) & |
|---|
| 32 | & - grazpof(:,:,jk) - ( 1.- unass ) * grazmf(:,:,jk) & |
|---|
| 33 | & + ( 1.- 0.5 * zfracal(:,:,jk) ) * ( tortnf(:,:,jk) & |
|---|
| 34 | & + respnf(:,:,jk) ) + 0.5 * tortdf(:,:,jk) + ferat3 * & |
|---|
| 35 | & ( tortz(:,:,jk) + respz(:,:,jk) ) - ofer(:,:,jk) & |
|---|
| 36 | & + ofer2(:,:,jk) - xaggfe(:,:,jk) & |
|---|
| 37 | & + xscave(:,:,jk) * zdenom1(:,:,jk) |
|---|
| 38 | |
|---|
| 39 | ! Evolution of big biogenic Iron |
|---|
| 40 | ! ------------------------------ |
|---|
| 41 | trn(:,:,jk,jpbfe) = trn(:,:,jk,jpbfe) + unass2 * ( graznf (:,:,jk) + grazf(:,:,jk) + grazfff(:,:,jk) & |
|---|
| 42 | & + grazpof(:,:,jk) + ferat3 * grazz(:,:,jk) ) & |
|---|
| 43 | & + ferat3 * ( tortz2 (:,:,jk) + respz2(:,:,jk) ) & |
|---|
| 44 | & - ofer2(:,:,jk) & |
|---|
| 45 | & + 0.5 * zfracal(:,:,jk) * ( respnf(:,:,jk) + tortnf(:,:,jk) ) & |
|---|
| 46 | & + 0.5 * tortdf (:,:,jk) + respdf(:,:,jk) + xaggfe(:,:,jk) & |
|---|
| 47 | & + xbactfer(:,:,jk) - grazfff(:,:,jk) + xscave(:,:,jk) * zdenom2(:,:,jk) |
|---|
| 48 | ! |
|---|
| 49 | END DO |
|---|
