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.std.h in trunk/NEMO/TOP_SRC/SMS – NEMO

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

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

insert revision Id

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 2.9 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
7CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
8CCC p4zbio : PISCES MODEL  - Standard parameterization
9CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
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)-xaggdoc2(ji,jj,jk)
23          END DO
24        END DO
25      END DO
26
27      DO jk = 1,jpkm1
28        DO jj = 1,jpj
29          DO ji = 1,jpi
30C     
31C     Evolution of Detritus
32C     ---------------------
33C     
34          trn(ji,jj,jk,jppoc) = trn(ji,jj,jk,jppoc)
35     &     -grazpoc(ji,jj,jk)+grapoc(ji,jj,jk)-grazm(ji,jj,jk)
36     &     +respz(ji,jj,jk)-xagg(ji,jj,jk)+xaggdoc(ji,jj,jk)
37     &     +(1.-0.5*zfracal(ji,jj,jk))*(tortp(ji,jj,jk)
38     &     +respp(ji,jj,jk))+0.5*tortp2(ji,jj,jk)
39     &     +orem2(ji,jj,jk)+tortz(ji,jj,jk)-orem(ji,jj,jk)
40C   
41C     Evolution of rapid Detritus
42C     ---------------------
43C   
44          trn(ji,jj,jk,jpgoc) = trn(ji,jj,jk,jpgoc)
45     &     +grapoc2(ji,jj,jk)+respp2(ji,jj,jk)+xagg(ji,jj,jk)
46     &     +tortz2(ji,jj,jk)+respz2(ji,jj,jk)-orem2(ji,jj,jk)
47     &     +0.5*zfracal(ji,jj,jk)*(respp(ji,jj,jk)+tortp(ji,jj,jk))
48     &     +0.5*tortp2(ji,jj,jk)+xaggdoc2(ji,jj,jk)-grazffe(ji,jj,jk)
49C
50          END DO
51        END DO
52      END DO
53
54
55
56      DO jk = 1,jpkm1
57        DO jj = 1,jpj
58          DO ji = 1,jpi
59C
60C     Evolution of small biogenic Iron
61C     --------------------------
62C
63          trn(ji,jj,jk,jpsfe) = trn(ji,jj,jk,jpsfe)
64     &     +unass*(grazpf(ji,jj,jk)+grazsf(ji,jj,jk))
65     &     -grazpof(ji,jj,jk)-(1.-unass)*grazmf(ji,jj,jk)
66     &     +(1.-0.5*zfracal(ji,jj,jk))*(tortnf(ji,jj,jk)
67     &     +respnf(ji,jj,jk))+0.5*tortdf(ji,jj,jk)+ferat3*
68     &     (tortz(ji,jj,jk)+respz(ji,jj,jk))-ofer(ji,jj,jk)
69     &     +ofer2(ji,jj,jk)-xaggfe(ji,jj,jk)
70     &     +xscave(ji,jj,jk)*zdenom1(ji,jj,jk)
71C
72C     Evolution of big biogenic Iron
73C     --------------------------
74C
75          trn(ji,jj,jk,jpbfe) = trn(ji,jj,jk,jpbfe)
76     &     +unass2*(graznf(ji,jj,jk)+grazf(ji,jj,jk)+grazfff(ji,jj,jk)
77     &     +grazpof(ji,jj,jk)+ferat3*grazz(ji,jj,jk))+ferat3*
78     &     (tortz2(ji,jj,jk)+respz2(ji,jj,jk))-ofer2(ji,jj,jk)
79     &     +0.5*zfracal(ji,jj,jk)*(respnf(ji,jj,jk)+tortnf(ji,jj,jk))
80     &     +0.5*tortdf(ji,jj,jk)+respdf(ji,jj,jk)+xaggfe(ji,jj,jk)
81     &     +xbactfer(ji,jj,jk)-grazfff(ji,jj,jk)+xscave(ji,jj,jk)
82     &     *zdenom2(ji,jj,jk)
83          END DO
84        END DO
85      END DO
86
Note: See TracBrowser for help on using the repository browser.