source: trunk/NEMO/TOP_SRC/PISCES/trclsm_pisces.F90 @ 1007

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

Update PISCES modules to take into account the re-organization of TOP initialization phase, see ticket 170

  • Property svn:executable set to *
File size: 5.1 KB
Line 
1MODULE trclsm_pisces
2   !!======================================================================
3   !!                      ***  MODULE trclsm_lobster  ***
4   !! TOP :   initialisation of some run parameters for PISCES bio-model
5   !!======================================================================
6   !! History :    -   !  1999-10 (M.A. Foujols, M. Levy) original code
7   !!              -   !  2000-01 (L. Bopp) hamocc3, p3zd
8   !!             1.0  !  2003-08 (C. Ethe)  module F90
9   !!             2.0  !  2007-12  (C. Ethe, G. Madec) from trclsm.pisces.h90
10   !!----------------------------------------------------------------------
11#if defined key_pisces
12   !!----------------------------------------------------------------------
13   !!   'key_pisces'   :                                   PISCES bio-model
14   !!----------------------------------------------------------------------
15   !! trc_lsm_pisces       : PISCES model namelist read
16   !!----------------------------------------------------------------------
17   USE oce_trc         ! Ocean variables
18   USE par_trc         ! TOP parameters
19   USE trc             ! TOP variables
20   USE sms             ! sms trends
21
22
23   IMPLICIT NONE
24   PRIVATE
25
26   PUBLIC   trc_lsm_pisces   ! called by trclsm.F90 module
27
28
29   !!----------------------------------------------------------------------
30   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007)
31   !! $Id: trclsm_pisces.F90 773 2007-12-17 16:55:32Z gm $
32   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
33   !!----------------------------------------------------------------------
34
35CONTAINS
36
37   SUBROUTINE trc_lsm_pisces
38      !!----------------------------------------------------------------------
39      !!                     ***  trc_lsm_pisces  *** 
40      !!
41      !! ** Purpose :   read PISCES namelist
42      !!
43      !! ** input   :   file 'namelist.trc.sms' containing the following
44      !!             namelist: natext, natbio, natsms
45      !!                       natkriest ("key_kriest")
46      !!----------------------------------------------------------------------
47      CHARACTER (len=32) ::   clname
48      !!
49      NAMELIST/natbio/ part, nrdttrc, wsbio, xkmort, ferat3, wsbio2
50
51#if defined key_kriest
52      NAMELIST/natkrpar/xkr_eta, xkr_zeta, xkr_mass_min, xkr_mass_max
53#endif
54      !!----------------------------------------------------------------------
55
56      IF(lwp) WRITE(numout,*)
57      IF(lwp) WRITE(numout,*) ' trc_lsm_pisces : read PISCES namelists'
58      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~~~~'
59
60
61      !                               ! Open the namelist file
62      !                               ! ----------------------
63      clname ='namelist.sms.pisces'
64      CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   &
65         &           1, numout, .FALSE., 1 )
66
67      !                               ! natbio : biological parameters
68      !                               ! ------------------------------
69      REWIND( numnat )                     ! read natbio
70      READ  ( numnat, natbio )
71
72      IF(lwp) THEN                         ! control print
73         WRITE(numout,*) ' Namelist : natbio'
74         WRITE(numout,*) '    part of calcite not dissolved in guts     part      =', part
75         WRITE(numout,*) '    frequence pour la biologie                nrdttrc   =', nrdttrc
76         WRITE(numout,*) '    POC sinking speed                         wsbio     =', wsbio
77         WRITE(numout,*) '    half saturation constant for mortality    xkmort    =', xkmort
78         WRITE(numout,*) '    Fe/C in zooplankton                       ferat3    =', ferat3
79         WRITE(numout,*) '    Big particles sinking speed               wsbio2    =', wsbio2
80      ENDIF
81
82#if defined key_kriest
83
84      !                               ! natkrpar : kriest parameters
85      !                               ! -----------------------------
86      xkr_eta      = 0.62       
87      xkr_zeta     = 1.62       
88      xkr_mass_min = 0.0002     
89      xkr_mass_max = 1.     
90
91      REWIND( numnat )                     ! read natkriest
92      READ  ( numnat, natkrpar )
93
94      IF(lwp) THEN
95         WRITE(numout,*)
96         WRITE(numout,*) ' Namelist : natkrpar'
97         WRITE(numout,*) '    Sinking  exponent                        xkr_eta      = ', xkr_eta
98         WRITE(numout,*) '    N content exponent                       xkr_zeta     = ', xkr_zeta
99         WRITE(numout,*) '    Minimum mass for Aggregates              xkr_mass_min = ', xkr_mass_min
100         WRITE(numout,*) '    Maximum mass for Aggregates              xkr_mass_max = ', xkr_mass_max
101         WRITE(numout,*)
102     ENDIF
103
104
105     ! Computation of some variables
106     xkr_massp = 5.7E-6 * 7.6 * xkr_mass_min**xkr_zeta
107
108#endif
109      !
110   END SUBROUTINE trc_lsm_pisces
111
112#else
113   !!----------------------------------------------------------------------
114   !!  Dummy module :                                   No PISCES bio-model
115   !!----------------------------------------------------------------------
116CONTAINS
117   SUBROUTINE trc_lsm_pisces                      ! Empty routine
118   END  SUBROUTINE  trc_lsm_pisces
119#endif 
120
121   !!======================================================================
122END MODULE trclsm_pisces
Note: See TracBrowser for help on using the repository browser.