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.
p4znano.F in branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS – NEMO

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, 16 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
RevLine 
[341]1
[719]2CCC $Header$ 
[341]3CCC  TOP 1.0 , LOCEAN-IPSL (2005) 
4C This software is governed by CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
5C ---------------------------------------------------------------------------
[186]6CDIR$ LIST
7      SUBROUTINE p4znano
[772]8#if defined key_top && defined key_pisces
[186]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
[339]49      REAL zfact,zstep,compaph
[186]50C
[339]51C      Time step duration for biology
52C      ------------------------------
53C
54        zstep=rfact2/rjjss
55C
[186]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.)
[339]61        zfact=1./(trn(ji,jj,jk,jpphy)+rtrn)
[186]62C
63C     Squared mortality of Phyto similar to a sedimentation term during
64C     blooms (Doney et al. 1996)
65C     -----------------------------------------------------------------
66C
[339]67        respp(ji,jj,jk) = wchl*1e6*zstep*zdiss(ji,jj,jk)
68     &    *compaph*trn(ji,jj,jk,jpphy)
[186]69#    if defined key_off_degrad
70     &    *facvol(ji,jj,jk)
71#    endif
[339]72                                                                               
[186]73        respnf(ji,jj,jk) = respp(ji,jj,jk)
[339]74     &    *trn(ji,jj,jk,jpnfe)*zfact
75                                                                               
[186]76        respnch(ji,jj,jk) = respp(ji,jj,jk)
[339]77     &    *trn(ji,jj,jk,jpnch)*zfact
78C
[186]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
[339]84        tortp(ji,jj,jk) = mprat*zstep*trn(ji,jj,jk,jpphy)
85     $    /(xkmort+trn(ji,jj,jk,jpphy))*compaph
[186]86#    if defined key_off_degrad
87     &    *facvol(ji,jj,jk)
88#    endif
[339]89                                                                               
[186]90        tortnf(ji,jj,jk)=tortp(ji,jj,jk)
[339]91     &    *trn(ji,jj,jk,jpnfe)*zfact
92                                                                               
[186]93        tortnch(ji,jj,jk)=tortp(ji,jj,jk)
[339]94     &    *trn(ji,jj,jk,jpnch)*zfact
[186]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.