source: trunk/NEMO/TOP_SRC/SMS/p4zbio.kriest.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: 3.0 KB
Line 
1
2
3CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
4CCC p4zbio : PISCES MODEL  - Kriest parameterization
5CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
6
7      DO jk = 1,jpkm1
8        DO jj = 1,jpj
9          DO ji = 1,jpi
10C   
11C     Evolution of DOC
12C     ----------------
13C     
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
27C     
28C     Evolution of Detritus
29C     ---------------------
30
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)
43C   
44C     Evolution of number of aggregates
45C     ---------------------------------
46C
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)
60C
61            trn(ji,jj,jk,jpnum) = max(trn(ji,jj,jk,jpnum),
62     &          trn(ji,jj,jk,jppoc)/xkr_massp/xnumm(jk))
63C
64            trn(ji,jj,jk,jpnum) = min(trn(ji,jj,jk,jpnum),
65     &          trn(ji,jj,jk,jppoc)/xkr_massp/1.1)
66C
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
77C
78C     Evolution of biogenic Iron
79C     --------------------------
80C
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
Note: See TracBrowser for help on using the repository browser.