source: trunk/NEMO/TOP_SRC/PISCES/p4zint.F90 @ 935

Last change on this file since 935 was 935, checked in by cetlod, 13 years ago

adding modules for PISCES SMS model, see ticket 141

  • Property svn:executable set to *
File size: 3.4 KB
Line 
1MODULE p4zint
2   !!======================================================================
3   !!                         ***  MODULE p4zint  ***
4   !! TOP :   PISCES interpolation and computation of various accessory fields
5   !!======================================================================
6   !! History :   1.0  !  2004-03 (O. Aumont) Original code
7   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
8   !!----------------------------------------------------------------------
9#if defined key_pisces
10   !!----------------------------------------------------------------------
11   !!   'key_pisces'                                       PISCES bio-model
12   !!----------------------------------------------------------------------
13   !!   p4z_int        :  interpolation and computation of various accessory fields
14   !!----------------------------------------------------------------------
15   USE oce_trc         !
16   USE trp_trc
17   USE sms
18
19   IMPLICIT NONE
20   PRIVATE
21
22   PUBLIC   p4z_int    ! called in p4zprg.F90
23
24   REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) ::   &   !:
25      Tgfunc,            &  !:  Temp. dependancy of various biological rates
26      Tgfunc2               !:  Temp. dependancy of mesozooplankton rates
27
28   !! * Module variables
29   REAL(wp) :: &
30      xksilim = 16.5E-6        ! Half-saturation constant for the computation of the Si half-saturation constant
31
32
33   !!----------------------------------------------------------------------
34   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
35   !! $Header:$
36   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
37   !!----------------------------------------------------------------------
38
39CONTAINS
40
41   SUBROUTINE p4z_int
42      !!---------------------------------------------------------------------
43      !!                     ***  ROUTINE p4z_int  ***
44      !!
45      !! ** Purpose :   interpolation and computation of various accessory fields
46      !!
47      !! ** Method  : - ???
48      !!---------------------------------------------------------------------
49      !!
50      INTEGER  ::   ji, jj
51      REAL(wp) ::   zdum
52      !!---------------------------------------------------------------------
53
54      ! Computation of phyto and zoo metabolic rate
55      ! -------------------------------------------
56
57      Tgfunc (:,:,:) = EXP( 0.063913 * tn(:,:,:) )
58      Tgfunc2(:,:,:) = EXP( 0.07608  * tn(:,:,:) )
59
60      ! Computation of the silicon dependant half saturation
61      ! constant for silica uptake
62      ! ---------------------------------------------------
63
64      DO ji = 1, jpi
65         DO jj = 1, jpj
66            zdum = trn(ji,jj,1,jpsil) * trn(ji,jj,1,jpsil)
67            xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zdum / ( xksilim * xksilim + zdum ) ) * 1e-6 )
68         END DO
69      END DO
70
71      IF( nday_year == 365 ) THEN
72         xksi    = xksimax
73         xksimax = 0.e0
74      ENDIF
75      !
76   END SUBROUTINE p4z_int
77
78#else
79   !!======================================================================
80   !!  Dummy module :                                   No PISCES bio-model
81   !!======================================================================
82CONTAINS
83   SUBROUTINE p4z_int                   ! Empty routine
84      WRITE(*,*) 'p4z_int: You should not have seen this print! error?'
85   END SUBROUTINE p4z_int
86#endif 
87
88   !!======================================================================
89END MODULE  p4zint
Note: See TracBrowser for help on using the repository browser.