source: branches/dev_001_GM/NEMO/TOP_SRC/PISCES_SMS/p4zint.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: 2.6 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 p4zint(kt)
8#if defined key_top && defined key_pisces
9CCC
10CCC
11CCC       ROUTINE p4zint : PISCES MODEL
12CCC       *****************************
13CCC
14CC
15CC   PURPOSE :
16CC   ---------
17CC           *P4ZINT* INTERPOLATION AND COMPUTATION OF 
18CC                    VARIOUS ACCESSORY FIELDS
19CC   INPUT :
20CC   -----
21CC      argument
22CC              kt              : time step
23CC
24CC   EXTERNAL :
25CC   ----------
26CC          NONE
27CC
28CC   MODIFICATIONS:
29CC   --------------
30CC      original      : 2004    O. Aumont
31CC ----------------------------------------------------------------
32CC parameters and commons
33CC ======================
34CDIR$ NOLIST
35      USE oce_trc
36      USE trp_trc
37      USE sms
38      IMPLICIT NONE
39CDIR$ LIST
40CC-----------------------------------------------------------------
41CC------
42CC local declarations
43CC ==================
44C
45      INTEGER kt
46      INTEGER ji, jj
47      INTEGER iman
48      INTEGER nspyr,nvit1t,nvit2t
49      REAL zpdtan, zman, zpdtmo, zdemi
50      REAL zt, zdum
51C
52C
53      zpdtan = raass / rdt
54      nspyr  = nint(zpdtan)
55      zman   = 12.
56      iman   = 12
57      zpdtmo = zpdtan / zman
58      zdemi  = zpdtmo / 2.
59      zt     = ( float ( kt) + zdemi) / zpdtmo
60     
61C  recherche de l'indice des enregistrements
62C  du modele dynamique encadrant le pas de temps kt.
63C  --------------------------------------------------
64C
65      xtvit = zt - float(int ( zt))
66      nvit1t = int (zt)
67      nvit2t = nvit1t+1
68      nvit1t = MOD ( nvit1t, iman)
69      IF ( nvit1t .EQ. 0 ) nvit1t = iman
70      nvit2t = MOD ( nvit2t, iman)
71      IF ( nvit2t .EQ. 0 ) nvit2t = iman
72C
73C Interpolation of dust deposition
74C --------------------------------
75C
76         dust(:,:) = (1.-xtvit)*dustmo(:,:,nvit1t)
77     $            +xtvit*dustmo(:,:,nvit2t)
78C
79C
80C Computation of phyto and zoo metabolic rate
81C -------------------------------------------
82C
83
84         Tgfunc(:,:,:) = exp(0.063913*tn(:,:,:))
85         Tgfunc2(:,:,:) = exp(0.07608*tn(:,:,:))
86C
87C      Computation of the silicon dependant half saturation
88C      constant for silica uptake
89C       ---------------------------------------------------
90C
91        DO ji=1,jpi
92          DO jj=1,jpj
93          zdum=trn(ji,jj,1,jpsil)**2
94          xksimax(ji,jj) = max(xksimax(ji,jj),(1.+7.*zdum
95     &      /(xksi2*xksi2*25.+zdum))*1E-6)
96          END DO
97        END DO
98C
99        IF (nday_year.EQ.365) THEN
100           xksi=xksimax
101           xksimax=0.
102        ENDIF
103C
104#endif
105C
106      RETURN
107      END
Note: See TracBrowser for help on using the repository browser.