source: trunk/NEMO/TOP_SRC/SMS/p4zbio.std.h @ 719

Last change on this file since 719 was 719, checked in by ctlod, 14 years ago

get back to the nemo_v2_3 version for trunk

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.7 KB
Line 
1CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2CCC p4zbio : PISCES MODEL  - Standard parameterization
3CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
4
5      DO jk = 1,jpkm1
6        DO jj = 1,jpj
7          DO ji = 1,jpi
8C   
9C     Evolution of DOC
10C     ----------------
11C     
12          trn(ji,jj,jk,jpdoc) = trn(ji,jj,jk,jpdoc)
13     &      +orem(ji,jj,jk)+excret2*prorca2(ji,jj,jk)
14     &      +excret*prorca(ji,jj,jk)-olimi(ji,jj,jk)-denitr(ji,jj,jk)
15     &      +grarem(ji,jj,jk)*(1.-sigma1)+grarem2(ji,jj,jk)
16     &      *(1.-sigma2)-xaggdoc(ji,jj,jk)-xaggdoc2(ji,jj,jk)
17          END DO
18        END DO
19      END DO
20
21      DO jk = 1,jpkm1
22        DO jj = 1,jpj
23          DO ji = 1,jpi
24C     
25C     Evolution of Detritus
26C     ---------------------
27C     
28          trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc)
29     &     -grazpoc(ji,jj,jk)+grapoc(ji,jj,jk)-grazm(ji,jj,jk)
30     &     +respz(ji,jj,jk)-xagg(ji,jj,jk)+xaggdoc(ji,jj,jk)
31     &     +(1.-0.5*zfracal(ji,jj,jk))*(tortp(ji,jj,jk)
32     &     +respp(ji,jj,jk))+0.5*tortp2(ji,jj,jk)
33     &     +orem2(ji,jj,jk)+tortz(ji,jj,jk)-orem(ji,jj,jk)
34C   
35C     Evolution of rapid Detritus
36C     ---------------------
37C   
38          trn(ji,jj,jk,jpgoc) = trn(ji,jj,jk,jpgoc)
39     &     +grapoc2(ji,jj,jk)+respp2(ji,jj,jk)+xagg(ji,jj,jk)
40     &     +tortz2(ji,jj,jk)+respz2(ji,jj,jk)-orem2(ji,jj,jk)
41     &     +0.5*zfracal(ji,jj,jk)*(respp(ji,jj,jk)+tortp(ji,jj,jk))
42     &     +0.5*tortp2(ji,jj,jk)+xaggdoc2(ji,jj,jk)-grazffe(ji,jj,jk)
43C
44          END DO
45        END DO
46      END DO
47
48
49
50      DO jk = 1,jpkm1
51        DO jj = 1,jpj
52          DO ji = 1,jpi
53C
54C     Evolution of small biogenic Iron
55C     --------------------------
56C
57          trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe)
58     &     +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk))
59     &     -grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk)
60     &     +(1.-0.5*zfracal(ji,jj,jk))*(tortnf(ji,jj,jk)
61     &     +respnf(ji,jj,jk))+0.5*tortdf(ji,jj,jk)+ferat3*
62     &     (tortz(ji,jj,jk)+respz(ji,jj,jk))-ofer(ji,jj,jk)
63     &     +ofer2(ji,jj,jk)-xaggfe(ji,jj,jk)
64     &     +xscave(ji,jj,jk)*zdenom1(ji,jj,jk)
65C
66C     Evolution of big biogenic Iron
67C     --------------------------
68C
69          trn(ji,jj,jk,jpbfe) = trn(ji,jj,jk,jpbfe)
70     &     +unass2*(graznf(ji,jj,jk)+grazf(ji,jj,jk)+grazfff(ji,jj,jk)
71     &     +grazpof(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3*
72     &     (tortz2(ji,jj,jk)+respz2(ji,jj,jk))-ofer2(ji,jj,jk)
73     &     +0.5*zfracal(ji,jj,jk)*(respnf(ji,jj,jk)+tortnf(ji,jj,jk))
74     &     +0.5*tortdf(ji,jj,jk)+respdf(ji,jj,jk)+xaggfe(ji,jj,jk)
75     &     +xbactfer(ji,jj,jk)-grazfff(ji,jj,jk)+xscave(ji,jj,jk)
76     &     *zdenom2(ji,jj,jk)
77          END DO
78        END DO
79      END DO
80
Note: See TracBrowser for help on using the repository browser.