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