1 | |
---|
2 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
3 | !CC p4zbio : PISCES MODEL - Kriest parameterization |
---|
4 | !CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC |
---|
5 | |
---|
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 |
---|
23 | |
---|
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) |
---|
32 | |
---|
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 ) |
---|
47 | |
---|
48 | trn(ji,jj,jk,jpnum) = MAX( trn(ji,jj,jk,jpnum), trn(ji,jj,jk,jppoc) / xkr_massp / xnumm(jk) ) |
---|
49 | |
---|
50 | trn(ji,jj,jk,jpnum) = MIN( trn(ji,jj,jk,jpnum), trn(ji,jj,jk,jppoc) / xkr_massp / 1.1 ) |
---|
51 | |
---|
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 |
---|
64 | END DO |
---|