source: branches/UKMO/r6232_tracer_advection/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zint.F90 @ 9295

Last change on this file since 9295 was 9295, checked in by jcastill, 3 years ago

Remove svn keywords

File size: 3.5 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         !  shared variables between ocean and passive tracers
16   USE trc             !  passive tracers common variables
17   USE sms_pisces      !  PISCES Source Minus Sink variables
18   USE iom
19
20   IMPLICIT NONE
21   PRIVATE
22
23   PUBLIC   p4z_int 
24   REAL(wp) ::   xksilim = 16.5e-6_wp   ! Half-saturation constant for the Si half-saturation constant computation
25
26   !!----------------------------------------------------------------------
27   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
28   !! $Id$
29   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
30   !!----------------------------------------------------------------------
31CONTAINS
32
33   SUBROUTINE p4z_int( kt )
34      !!---------------------------------------------------------------------
35      !!                     ***  ROUTINE p4z_int  ***
36      !!
37      !! ** Purpose :   interpolation and computation of various accessory fields
38      !!
39      !!---------------------------------------------------------------------
40      !
41      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index
42      !
43      INTEGER  :: ji, jj                 ! dummy loop indices
44      REAL(wp) :: zvar                   ! local variable
45      !!---------------------------------------------------------------------
46      !
47      IF( nn_timing == 1 )  CALL timing_start('p4z_int')
48      !
49      ! Computation of phyto and zoo metabolic rate
50      ! -------------------------------------------
51      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) )
52      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) )
53
54      ! Computation of the silicon dependant half saturation  constant for silica uptake
55      ! ---------------------------------------------------
56      DO ji = 1, jpi
57         DO jj = 1, jpj
58            zvar = trb(ji,jj,1,jpsil) * trb(ji,jj,1,jpsil)
59            xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zvar / ( xksilim * xksilim + zvar ) ) * 1e-6 )
60         END DO
61      END DO
62      !
63      IF( nday_year == nyear_len(1) ) THEN
64         xksi   (:,:) = xksimax(:,:)
65         xksimax(:,:) = 0._wp
66      ENDIF
67      !
68      IF( nn_timing == 1 )  CALL timing_stop('p4z_int')
69      !
70   END SUBROUTINE p4z_int
71
72#else
73   !!======================================================================
74   !!  Dummy module :                                   No PISCES bio-model
75   !!======================================================================
76CONTAINS
77   SUBROUTINE p4z_int                   ! Empty routine
78      WRITE(*,*) 'p4z_int: You should not have seen this print! error?'
79   END SUBROUTINE p4z_int
80#endif 
81
82   !!======================================================================
83END MODULE p4zint
Note: See TracBrowser for help on using the repository browser.