source: NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zint.F90 @ 10368

Last change on this file since 10368 was 10068, checked in by nicolasmartin, 3 years ago

First part of modifications to have a common default header : fix typos and SVN keywords properties

  • Property svn:keywords set to Id
File size: 2.8 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   !!   p4z_int        :  interpolation and computation of various accessory fields
10   !!----------------------------------------------------------------------
11   USE oce_trc         !  shared variables between ocean and passive tracers
12   USE trc             !  passive tracers common variables
13   USE sms_pisces      !  PISCES Source Minus Sink variables
14
15   IMPLICIT NONE
16   PRIVATE
17
18   PUBLIC   p4z_int 
19   REAL(wp) ::   xksilim = 16.5e-6_wp   ! Half-saturation constant for the Si half-saturation constant computation
20
21   !!----------------------------------------------------------------------
22   !! NEMO/TOP 4.0 , NEMO Consortium (2018)
23   !! $Id$
24   !! Software governed by the CeCILL license (see ./LICENSE)
25   !!----------------------------------------------------------------------
26CONTAINS
27
28   SUBROUTINE p4z_int( kt )
29      !!---------------------------------------------------------------------
30      !!                     ***  ROUTINE p4z_int  ***
31      !!
32      !! ** Purpose :   interpolation and computation of various accessory fields
33      !!
34      !!---------------------------------------------------------------------
35      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index
36      !
37      INTEGER  :: ji, jj                 ! dummy loop indices
38      REAL(wp) :: zvar                   ! local variable
39      !!---------------------------------------------------------------------
40      !
41      IF( ln_timing )   CALL timing_start('p4z_int')
42      !
43      ! Computation of phyto and zoo metabolic rate
44      ! -------------------------------------------
45      tgfunc (:,:,:) = EXP( 0.063913 * tsn(:,:,:,jp_tem) )
46      tgfunc2(:,:,:) = EXP( 0.07608  * tsn(:,:,:,jp_tem) )
47
48      ! Computation of the silicon dependant half saturation  constant for silica uptake
49      ! ---------------------------------------------------
50      DO ji = 1, jpi
51         DO jj = 1, jpj
52            zvar = trb(ji,jj,1,jpsil) * trb(ji,jj,1,jpsil)
53            xksimax(ji,jj) = MAX( xksimax(ji,jj), ( 1.+ 7.* zvar / ( xksilim * xksilim + zvar ) ) * 1e-6 )
54         END DO
55      END DO
56      !
57      IF( nday_year == nyear_len(1) ) THEN
58         xksi   (:,:) = xksimax(:,:)
59         xksimax(:,:) = 0._wp
60      ENDIF
61      !
62      IF( ln_timing )   CALL timing_stop('p4z_int')
63      !
64   END SUBROUTINE p4z_int
65
66   !!======================================================================
67END MODULE p4zint
Note: See TracBrowser for help on using the repository browser.