source: trunk/NEMO/TOP_SRC/PISCES/trcini_pisces.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

File size: 3.8 KB
Line 
1MODULE trcini_pisces
2   !!======================================================================
3   !!                         ***  MODULE trcini_pisces  ***
4   !! TOP :   initialisation of the PISCES biochemical model
5   !!======================================================================
6   !! History :    -   !  1988-07  (E. Maier-Reiner) Original code
7   !!              -   !  1999-10  (O. Aumont, C. Le Quere)
8   !!              -   !  2002     (O. Aumont)  PISCES
9   !!             1.0  !  2005-03  (O. Aumont, A. El Moussaoui) F90
10   !!             2.0  !  2007-12  (C. Ethe, G. Madec) from trcini.pisces.h90
11   !!----------------------------------------------------------------------
12#if defined key_pisces
13   !!----------------------------------------------------------------------
14   !!   'key_pisces'                                       PISCES bio-model
15   !!----------------------------------------------------------------------
16   !! trc_ini_pisces   : PISCES biochemical model initialisation
17   !!----------------------------------------------------------------------
18   USE par_trc         ! TOP parameters
19   USE sms             ! Source Minus Sink variables
20   USE oce_trc         ! ocean variables
21   USE trp_trc         !
22   USE p4zche 
23
24
25   IMPLICIT NONE
26   PRIVATE
27
28   PUBLIC   trc_ini_pisces   ! called by trcini.F90 module
29
30
31#  include "domzgr_substitute.h90"
32#  include "top_substitute.h90"
33   !!----------------------------------------------------------------------
34   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
35   !! $Id: trcini_pisces.F90 776 2007-12-19 14:10:14Z gm $
36   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
37   !!----------------------------------------------------------------------
38
39CONTAINS
40
41   SUBROUTINE trc_ini_pisces
42      !!----------------------------------------------------------------------
43      !!                   ***  ROUTINE trc_ini_pisces ***
44      !!
45      !! ** Purpose :   Initialisation of the PISCES biochemical model
46      !!----------------------------------------------------------------------
47      INTEGER ::   jk
48
49      !!----------------------------------------------------------------------
50
51
52      IF(lwp) WRITE(numout,*)
53      IF(lwp) WRITE(numout,*) ' trc_ini_pisces :   PISCES biochemical model initialisation'
54      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
55
56
57      !                                            ! Time-step
58      rfact   = rdttra(1) * float(ndttrc)          ! ---------
59      rfactr  = 1. / rfact
60      rfact2  = rfact / float(nrdttrc)
61      rfact2r = 1. / rfact2
62
63      IF(lwp) WRITE(numout,*) '    Tracer  time step    rfact  = ', rfact, ' rdt = ', rdt
64      IF(lwp) write(numout,*) '    Biology time step    rfact2 = ', rfact2
65
66
67      !----------------------------------------------------------------------
68      ! Initialize biological variables
69      !----------------------------------------------------------------------
70      ! Set biological ratios
71      ! ---------------------
72      rno3   = (16.+2.) / 122.
73      po4r   =   1.e0   / 122.
74      o2nit  =  32.     / 122.
75      rdenit =  97.6    /  16.
76      o2ut   = 140.     / 122.
77
78      CALL p4z_che        ! initialize the chemical constants
79
80      ndayflxtr = 0      !  Initialize a counter for the computation of chemistry
81
82      IF(lwp) WRITE(numout,*) ' Initialisation of PISCES done'
83      !
84   END SUBROUTINE trc_ini_pisces
85   
86#else
87   !!----------------------------------------------------------------------
88   !!   Dummy module                            No PISCES biochemical model
89   !!----------------------------------------------------------------------
90CONTAINS
91   SUBROUTINE trc_ini_pisces             ! Empty routine
92   END SUBROUTINE trc_ini_pisces
93#endif
94
95   !!======================================================================
96END MODULE trcini_pisces
Note: See TracBrowser for help on using the repository browser.