New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
p4zbio.kriest.h in trunk/NEMO/TOP_SRC/SMS – NEMO

source: trunk/NEMO/TOP_SRC/SMS/p4zbio.kriest.h @ 699

Last change on this file since 699 was 699, checked in by smasson, 17 years ago

insert revision Id
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 3.2 KB
Line 
1CC----------------------------------------------------------------------
2CC  TOP 1.0 , LOCEAN-IPSL (2005) 
3CC $Id$
4CC This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt
5CC----------------------------------------------------------------------
6
7CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
8CCC p4zbio : PISCES MODEL  - Kriest parameterization
9CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
10
11      DO jk = 1,jpkm1
12        DO jj = 1,jpj
13          DO ji = 1,jpi
14C   
15C     Evolution of DOC
16C     ----------------
17C     
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)
23          END DO
24        END DO
25      END DO
26
27
28      DO jk = 1,jpkm1
29        DO jj = 1,jpj
30          DO ji = 1,jpi
31C     
32C     Evolution of Detritus
33C     ---------------------
34C 
35            znumpoc=trn(ji,jj,jk,jpnum)/(trn(ji,jj,jk,jppoc)+rtrn)
36            znumdoc=znumpoc
37
38            trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc)
39     &          -grazpoc(ji,jj,jk)
40     &          +grapoc(ji,jj,jk)+grapoc2(ji,jj,jk)
41     &          -grazm(ji,jj,jk)+tortz2(ji,jj,jk)
42     &          +respz(ji,jj,jk)+respz2(ji,jj,jk)
43     &          +respp(ji,jj,jk)+respp2(ji,jj,jk)
44     &          +tortp2(ji,jj,jk)+tortz(ji,jj,jk)
45     &          +tortp(ji,jj,jk)-orem(ji,jj,jk)
46     &          +xaggdoc(ji,jj,jk)-grazffe(ji,jj,jk)
47C   
48C     Evolution of number of aggregates
49C     ---------------------------------
50C
51            trn(ji,jj,jk,jpnum) = trn(ji,jj,jk,jpnum)
52     &          -xagg(ji,jj,jk)-(orem(ji,jj,jk)
53     &          +grazpoc(ji,jj,jk))*znumpoc
54     &          +(tortp(ji,jj,jk))*xkr_nnano
55     &          +(tortp2(ji,jj,jk)+respp(ji,jj,jk)+tortz(ji,jj,jk)
56     &          +grapoc(ji,jj,jk)-grazm(ji,jj,jk)
57     &          +respz(ji,jj,jk))*xkr_ndiat
58     &          +(grapoc2(ji,jj,jk)+tortz2(ji,jj,jk)
59     &          +respz2(ji,jj,jk))*xkr_nmeso
60     &          +respp2(ji,jj,jk)*xkr_naggr
61     &          +xaggdoc(ji,jj,jk)*znumdoc
62     &          -grazffe(ji,jj,jk)*znumpoc*wsbio4(ji,jj,jk)
63     &          /(wsbio3(ji,jj,jk)+rtrn)
64C
65            trn(ji,jj,jk,jpnum) = max(trn(ji,jj,jk,jpnum),
66     &          trn(ji,jj,jk,jppoc)/xkr_massp/xnumm(jk))
67C
68            trn(ji,jj,jk,jpnum) = min(trn(ji,jj,jk,jpnum),
69     &          trn(ji,jj,jk,jppoc)/xkr_massp/1.1)
70C
71
72          END DO
73        END DO
74      END DO
75
76
77
78      DO jk = 1,jpkm1
79        DO jj = 1,jpj
80          DO ji = 1,jpi
81C
82C     Evolution of biogenic Iron
83C     --------------------------
84C
85          trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe)
86     &     +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk))
87     &     -(1.-unass2)*grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk)
88     &     -(1.-unass2)*grazfff(ji,jj,jk)+unass2*(graznf(ji,jj,jk)
89     &     +grazf(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3
90     &     *(tortz2(ji,jj,jk)+respz2(ji,jj,jk)+tortz(ji,jj,jk)
91     &     +respz(ji,jj,jk))-ofer(ji,jj,jk)+(respnf(ji,jj,jk)
92     &     +tortnf(ji,jj,jk))+tortdf(ji,jj,jk)+respdf(ji,jj,jk)
93     &     +xbactfer(ji,jj,jk)+xscave(ji,jj,jk)*zdenom1(ji,jj,jk)
94          END DO
95        END DO
96      END DO
97
Note: See TracBrowser for help on using the repository browser.