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.
p4zint.F90 in NEMO/trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z – NEMO

source: NEMO/trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zint.F90 @ 9594

Last change on this file since 9594 was 9124, checked in by gm, 6 years ago

dev_merge_2017: ln_timing instead of nn_timing + restricted timing to nemo_init and routine called by step in OPA_SRC

  • Property svn:keywords set to Id
File size: 2.9 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 3.3 , NEMO Consortium (2010)
23   !! $Id$
24   !! Software governed by the CeCILL licence     (NEMOGCM/NEMO_CeCILL.txt)
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.