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 | |
---|
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 | |
---|