source: branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4znano.F @ 772

Last change on this file since 772 was 772, checked in by gm, 13 years ago

dev_001_GM - change the name of cpp key to key_top, key_lobster, key_pisces, key_kriest and the corresponding lk_

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 3.0 KB
Line 
1
2CCC $Header$ 
3CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
4C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
5C ---------------------------------------------------------------------------
6CDIR$ LIST
7      SUBROUTINE p4znano
8#if defined key_top && defined key_pisces
9CCC---------------------------------------------------------------------
10CCC
11CCC          ROUTINE p4znano : PISCES MODEL
12CCC          ******************************
13CCC
14CCC  PURPOSE :
15CCC  ---------
16CCC         Compute the mortality terms for nanophytoplankton
17CCC
18CC   INPUT :
19CC   -----
20CC      argument
21CC              None
22CC      common
23CC              all the common defined in opa
24CC
25CC
26CC   OUTPUT :                   : no
27CC   ------
28CC
29CC   EXTERNAL :
30CC   --------
31CC          None
32CC
33CC   MODIFICATIONS:
34CC   --------------
35CC      original  : O. Aumont (2002)
36CC----------------------------------------------------------------------
37CC parameters and commons
38CC ======================
39CDIR$ NOLIST
40      USE oce_trc
41      USE trp_trc
42      USE sms
43      IMPLICIT NONE
44CDIR$ LIST
45CC----------------------------------------------------------------------
46CC local declarations
47CC ==================
48      INTEGER ji, jj, jk
49      REAL zfact,zstep,compaph
50C
51C      Time step duration for biology
52C      ------------------------------
53C
54        zstep=rfact2/rjjss
55C
56        DO jk = 1,jpkm1
57          DO jj = 1,jpj
58            DO ji = 1,jpi
59C
60        compaph = max((trn(ji,jj,jk,jpphy)-1E-8),0.)
61        zfact=1./(trn(ji,jj,jk,jpphy)+rtrn)
62C
63C     Squared mortality of Phyto similar to a sedimentation term during
64C     blooms (Doney et al. 1996)
65C     -----------------------------------------------------------------
66C
67        respp(ji,jj,jk) = wchl*1e6*zstep*zdiss(ji,jj,jk)
68     &    *compaph*trn(ji,jj,jk,jpphy)
69#    if defined key_off_degrad
70     &    *facvol(ji,jj,jk)
71#    endif
72                                                                               
73        respnf(ji,jj,jk) = respp(ji,jj,jk)
74     &    *trn(ji,jj,jk,jpnfe)*zfact
75                                                                               
76        respnch(ji,jj,jk) = respp(ji,jj,jk)
77     &    *trn(ji,jj,jk,jpnch)*zfact
78C
79C     Phytoplankton mortality. This mortality loss is slightly
80C     increased when nutrients are limiting phytoplankton growth
81C     as observed for instance in case of iron limitation.
82C     ----------------------------------------------------------
83C
84        tortp(ji,jj,jk) = mprat*zstep*trn(ji,jj,jk,jpphy)
85     $    /(xkmort+trn(ji,jj,jk,jpphy))*compaph
86#    if defined key_off_degrad
87     &    *facvol(ji,jj,jk)
88#    endif
89                                                                               
90        tortnf(ji,jj,jk)=tortp(ji,jj,jk)
91     &    *trn(ji,jj,jk,jpnfe)*zfact
92                                                                               
93        tortnch(ji,jj,jk)=tortp(ji,jj,jk)
94     &    *trn(ji,jj,jk,jpnch)*zfact
95C
96            END DO
97          END DO
98        END DO
99C
100#endif
101      RETURN
102      END
Note: See TracBrowser for help on using the repository browser.