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.
p5zmicro.F90 in branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z – NEMO

source: branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zmicro.F90 @ 7180

Last change on this file since 7180 was 7180, checked in by aumont, 8 years ago

various bug fixes on iron chemistry

  • Property svn:executable set to *
File size: 23.0 KB
Line 
1MODULE p5zmicro
2   !!======================================================================
3   !!                         ***  MODULE p5zmicro  ***
4   !! TOP :   PISCES Compute the sources/sinks for microzooplankton
5   !!======================================================================
6   !! History :   1.0  !  2004     (O. Aumont) Original code
7   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
8   !!             3.4  !  2011-06  (O. Aumont, C. Ethe) Quota model for iron
9   !!             3.6  !  2015-05  (O. Aumont) PISCES quota
10   !!----------------------------------------------------------------------
11#if defined key_pisces_quota
12   !!----------------------------------------------------------------------
13   !!   'key_pisces_quota'                                 PISCES bio-model
14   !!----------------------------------------------------------------------
15   !!   p5z_micro       :   Compute the sources/sinks for microzooplankton
16   !!   p5z_micro_init  :   Initialize and read the appropriate namelist
17   !!----------------------------------------------------------------------
18   USE oce_trc         !  shared variables between ocean and passive tracers
19   USE trc             !  passive tracers common variables
20   USE sms_pisces      !  PISCES Source Minus Sink variables
21   USE p5zsink         !  vertical flux of particulate matter due to sinking
22   USE p5zlim          !  Phytoplankton limitation terms
23   USE iom             !  I/O manager
24   USE prtctl_trc      !  print control for debugging
25
26   IMPLICIT NONE
27   PRIVATE
28
29   PUBLIC   p5z_micro         ! called in p5zbio.F90
30   PUBLIC   p5z_micro_init    ! called in trcsms_pisces.F90
31
32   !! * Shared module variables
33   REAL(wp), PUBLIC ::  part        !: part of calcite not dissolved in microzoo guts
34   REAL(wp), PUBLIC ::  xprefc     !: microzoo preference for POC
35   REAL(wp), PUBLIC ::  xprefn     !: microzoo preference for nanophyto
36   REAL(wp), PUBLIC ::  xprefp     !: microzoo preference for picophyto
37   REAL(wp), PUBLIC ::  xprefd     !: microzoo preference for diatoms
38   REAL(wp), PUBLIC ::  xprefz     !: microzoo preference for microzoo
39   REAL(wp), PUBLIC ::  xthreshdia  !: diatoms feeding threshold for microzooplankton
40   REAL(wp), PUBLIC ::  xthreshpic  !: picophyto feeding threshold for microzooplankton
41   REAL(wp), PUBLIC ::  xthreshphy  !: nanophyto threshold for microzooplankton
42   REAL(wp), PUBLIC ::  xthreshzoo  !: microzoo threshold for microzooplankton
43   REAL(wp), PUBLIC ::  xthreshpoc  !: poc threshold for microzooplankton
44   REAL(wp), PUBLIC ::  xthresh     !: feeding threshold for microzooplankton
45   REAL(wp), PUBLIC ::  resrat      !: exsudation rate of microzooplankton
46   REAL(wp), PUBLIC ::  mzrat       !: microzooplankton mortality rate
47   REAL(wp), PUBLIC ::  grazrat     !: maximal microzoo grazing rate
48   REAL(wp), PUBLIC ::  xkgraz      !: non assimilated fraction of P by microzoo
49   REAL(wp), PUBLIC ::  unassc      !: Efficicency of microzoo growth
50   REAL(wp), PUBLIC ::  unassn      !: Efficicency of microzoo growth
51   REAL(wp), PUBLIC ::  unassp      !: Efficicency of microzoo growth
52   REAL(wp), PUBLIC ::  epsher      !: half sturation constant for grazing 1
53   REAL(wp), PUBLIC ::  srespir     !: half sturation constant for grazing 1
54   REAL(wp), PUBLIC ::  ssigma      !: Fraction excreted as semi-labile DOM
55   LOGICAL,  PUBLIC ::  bmetexc     !: Use of excess carbon for respiration
56
57
58   !!* Substitution
59#  include "top_substitute.h90"
60   !!----------------------------------------------------------------------
61   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
62   !! $Id: p4zmicro.F90 3160 2011-11-20 14:27:18Z cetlod $
63   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
64   !!----------------------------------------------------------------------
65
66CONTAINS
67
68   SUBROUTINE p5z_micro( kt, knt )
69      !!---------------------------------------------------------------------
70      !!                     ***  ROUTINE p5z_micro  ***
71      !!
72      !! ** Purpose :   Compute the sources/sinks for microzooplankton
73      !!
74      !! ** Method  : - ???
75      !!---------------------------------------------------------------------
76      INTEGER, INTENT(in) ::  kt  ! ocean time step
77      INTEGER, INTENT(in) ::  knt 
78      !
79      INTEGER  :: ji, jj, jk
80      REAL(wp) :: zcompadi, zcompaz , zcompaph, zcompapoc, zcompapon, zcompapop
81      REAL(wp) :: zcompapi, zgraze  , zdenom, zfact, zstep, zfood, zfoodlim
82      REAL(wp) :: ztmp1, ztmp2, ztmp3, ztmp4, ztmp5, ztmptot
83      REAL(wp) :: zepsherf, zepshert, zepsherv, zrespirc, zrespirn, zrespirp, zbasresb, zbasresi
84      REAL(wp) :: zgraztotc, zgraztotn, zgraztotp, zgraztotf, zbasresn, zbasresp, zbasresf
85      REAL(wp) :: zgradoc, zgradon, zgradop, zgraref, zgradoct, zgradont, zgradopt, zgrareft
86      REAL(wp) :: zexcess, zgraren, zgrarep, zgrarem
87      REAL(wp) :: zgrapoc, zgrapon, zgrapop, zgrapof, zprcaca, zmortz
88      REAL(wp) :: zrespz, ztortz, zgrasratf, zgrasratn, zgrasratp
89      REAL(wp) :: zgraznc, zgraznn, zgraznp, zgrazpoc, zgrazpon, zgrazpop, zgrazpof
90      REAL(wp) :: zgrazdc, zgrazdn, zgrazdp, zgrazdf, zgraznf, zgrazz
91      REAL(wp) :: zgrazpc, zgrazpn, zgrazpp, zgrazpf, zbeta, zrfact2, zmetexcess
92      REAL(wp), POINTER, DIMENSION(:,:,:) :: zgrazing, zw3d
93      CHARACTER (len=25) :: charout
94      !!---------------------------------------------------------------------
95      !
96      IF( nn_timing == 1 )  CALL timing_start('p5z_micro')
97      !
98      IF( lk_iomput )  CALL wrk_alloc( jpi, jpj, jpk, zgrazing )
99      !
100      zmetexcess = 0.0
101      IF ( bmetexc ) zmetexcess = 1.0
102      !
103      DO jk = 1, jpkm1
104         DO jj = 1, jpj
105            DO ji = 1, jpi
106               zcompaz = MAX( ( trb(ji,jj,jk,jpzoo) - 1.e-9 ), 0.e0 )
107               zstep   = xstep
108# if defined key_degrad
109               zstep = zstep * facvol(ji,jj,jk)
110# endif
111               zfact   = zstep * tgfunc2(ji,jj,jk) * zcompaz
112
113               !   Michaelis-Menten mortality rates of microzooplankton
114               !   -----------------------------------------------------
115               zrespz = resrat * zfact * ( trb(ji,jj,jk,jpzoo) / ( xkmort + trb(ji,jj,jk,jpzoo) )  &
116               &        + 3. * nitrfac(ji,jj,jk) )
117
118               !   Zooplankton mortality. A square function has been selected with
119               !   no real reason except that it seems to be more stable and may mimic predation.
120               !   ------------------------------------------------------------------------------
121               ztortz = mzrat * 1.e6 * zfact * trb(ji,jj,jk,jpzoo)
122
123               !   Computation of the abundance of the preys
124               !   A threshold can be specified in the namelist
125               !   --------------------------------------------
126               zcompadi  = MIN( MAX( ( trb(ji,jj,jk,jpdia) - xthreshdia ), 0.e0 ), xsizedia )
127               zcompaph  = MAX( ( trb(ji,jj,jk,jpphy) - xthreshphy ), 0.e0 )
128               zcompaz   = MAX( ( trb(ji,jj,jk,jpzoo) - xthreshzoo ), 0.e0 )
129               zcompapi  = MAX( ( trb(ji,jj,jk,jppic) - xthreshpic ), 0.e0 )
130               zcompapoc = MAX( ( trb(ji,jj,jk,jppoc) - xthreshpoc ), 0.e0 )
131               
132               !   Microzooplankton grazing
133               !   ------------------------
134               zfood     = xprefn * zcompaph + xprefc * zcompapoc + xprefd * zcompadi   &
135               &           + xprefz * zcompaz + xprefp * zcompapi
136               zfoodlim  = MAX( 0. , zfood - min(xthresh,0.5*zfood) )
137               zdenom    = zfoodlim / ( xkgraz + zfoodlim )
138               zgraze    = grazrat * zstep * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpzoo) 
139
140               !   An active switching parameterization is used here.
141               !   We don't use the KTW parameterization proposed by
142               !   Vallina et al. because it tends to produce to steady biomass
143               !   composition and the variance of Chl is too low as it grazes
144               !   too strongly on winning organisms. Thus, instead of a square
145               !   a 1.5 power value is used which decreases the pressure on the
146               !   most abundant species
147               !   ------------------------------------------------------------ 
148               ztmp1 = xprefn * zcompaph**1.5
149               ztmp2 = xprefp * zcompapi**1.5
150               ztmp3 = xprefc * zcompapoc**1.5
151               ztmp4 = xprefd * zcompadi**1.5
152               ztmp5 = xprefz * zcompaz**1.5
153               ztmptot = ztmp1 + ztmp2 + ztmp3 + ztmp4 + ztmp5 + rtrn
154               ztmp1 = ztmp1 / ztmptot
155               ztmp2 = ztmp2 / ztmptot
156               ztmp3 = ztmp3 / ztmptot
157               ztmp4 = ztmp4 / ztmptot
158               ztmp5 = ztmp5 / ztmptot
159
160               !   Microzooplankton regular grazing on the different preys
161               !   -------------------------------------------------------
162               zgraznc   = zgraze  * ztmp1  * zdenom
163               zgraznn   = zgraznc * trb(ji,jj,jk,jpnph) / (trb(ji,jj,jk,jpphy) + rtrn)
164               zgraznp   = zgraznc * trb(ji,jj,jk,jppph) / (trb(ji,jj,jk,jpphy) + rtrn)
165               zgraznf   = zgraznc * trb(ji,jj,jk,jpnfe) / (trb(ji,jj,jk,jpphy) + rtrn)
166               zgrazpc   = zgraze  * ztmp2  * zdenom
167               zgrazpn   = zgrazpc * trb(ji,jj,jk,jpnpi) / (trb(ji,jj,jk,jppic) + rtrn)
168               zgrazpp   = zgrazpc * trb(ji,jj,jk,jpppi) / (trb(ji,jj,jk,jppic) + rtrn)
169               zgrazpf   = zgrazpc * trb(ji,jj,jk,jppfe) / (trb(ji,jj,jk,jppic) + rtrn)
170               zgrazz    = zgraze  * ztmp5   * zdenom
171               zgrazpoc  = zgraze  * ztmp3   * zdenom
172               zgrazpon  = zgrazpoc * trb(ji,jj,jk,jppon) / ( trb(ji,jj,jk,jppoc) + rtrn )
173               zgrazpop  = zgrazpoc * trb(ji,jj,jk,jppop) / ( trb(ji,jj,jk,jppoc) + rtrn )
174               zgrazpof  = zgrazpoc* trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
175               zgrazdc   = zgraze  * ztmp4  * zdenom
176               zgrazdn   = zgrazdc * trb(ji,jj,jk,jpndi) / (trb(ji,jj,jk,jpdia) + rtrn)
177               zgrazdp   = zgrazdc * trb(ji,jj,jk,jppdi) / (trb(ji,jj,jk,jpdia) + rtrn)
178               zgrazdf   = zgrazdc * trb(ji,jj,jk,jpdfe) / (trb(ji,jj,jk,jpdia) + rtrn)
179               !
180               zgraztotc = zgraznc + zgrazpoc + zgrazdc + zgrazz + zgrazpc
181               zgraztotn = zgraznn + zgrazpn + zgrazpon + zgrazdn + zgrazz * no3rat3
182               zgraztotp = zgraznp + zgrazpp + zgrazpop + zgrazdp + zgrazz * po4rat3
183               zgraztotf = zgraznf + zgrazpf + zgrazpof + zgrazdf + zgrazz * ferat3
184               !
185               ! Grazing by microzooplankton
186               IF( ln_diatrc .AND. lk_iomput )  zgrazing(ji,jj,jk) = zgraztotc
187
188               !   Stoichiometruc ratios of the food ingested by zooplanton
189               !   --------------------------------------------------------
190               zgrasratf =  (zgraztotf + rtrn) / ( zgraztotc + rtrn )
191               zgrasratn =  (zgraztotn + rtrn) / ( zgraztotc + rtrn )
192               zgrasratp =  (zgraztotp + rtrn) / ( zgraztotc + rtrn )
193
194               !   Growth efficiency is made a function of the quality
195               !   and the quantity of the preys
196               !   ---------------------------------------------------
197               zepshert  = MIN( 1., zgrasratn/ no3rat3, zgrasratp/ po4rat3, zgrasratf / ferat3)
198               zbeta = 1./ (epsher - 0.2)
199               zepsherf = 0.2 + 1./ (zbeta + 0.04 * 12. * zfood * 1E6 )
200               zepsherv  = zepsherf * zepshert
201
202               !   Respiration of microzooplankton
203               !   Excess carbon in the food is used preferentially
204               !   ------------------------------------------------
205               zexcess  = zgraztotc * zepsherf * (1.0 - zepshert) * zmetexcess
206               zbasresb = MAX(0., zrespz - zexcess)
207               zbasresi = zexcess + MIN(0., zrespz - zexcess) 
208               zrespirc = srespir * zepsherv * zgraztotc + zbasresb
209               
210               !   When excess carbon is used, the other elements in excess
211               !   are also used proportionally to their abundance
212               !   --------------------------------------------------------
213               zexcess  = ( zgrasratn/ no3rat3 - zepshert ) / ( 1.0 - zepshert + rtrn)
214               zbasresn = zbasresi * zexcess * zgrasratn 
215               zexcess  = ( zgrasratp/ po4rat3 - zepshert ) / ( 1.0 - zepshert + rtrn)
216               zbasresp = zbasresi * zexcess * zgrasratp
217               zexcess  = ( zgrasratf/ ferat3 - zepshert ) / ( 1.0 - zepshert + rtrn)
218               zbasresf = zbasresi * zexcess * zgrasratf
219
220               !   Voiding of the excessive elements as DOM
221               !   ----------------------------------------
222               zgradoct   = (1. - unassc - zepsherv) * zgraztotc - zbasresi 
223               zgradont   = (1. - unassn) * zgraztotn - zepsherv * no3rat3 * zgraztotc - zbasresn
224               zgradopt   = (1. - unassp) * zgraztotp - zepsherv * po4rat3 * zgraztotc - zbasresp
225               zgrareft   = (1. - unassc) * zgraztotf - zepsherv * ferat3 * zgraztotc - zbasresf
226
227               !  Since only semilabile DOM is represented in PISCES
228               !  part of DOM is in fact labile and is then released
229               !  as dissolved inorganic compounds (ssigma)
230               !  --------------------------------------------------
231               zgradoc =  zgradoct * ssigma
232               zgradon =  zgradont * ssigma
233               zgradop =  zgradopt * ssigma
234               zgrarem = (1.0 - ssigma) * zgradoct
235               zgraren = (1.0 - ssigma) * zgradont
236               zgrarep = (1.0 - ssigma) * zgradopt
237               zgraref = zgrareft
238
239               !   Defecation as a result of non assimilated products
240               !   --------------------------------------------------
241               zgrapoc   = zgraztotc * unassc
242               zgrapon   = zgraztotn * unassn
243               zgrapop   = zgraztotp * unassp
244               zgrapof   = zgraztotf * unassc
245
246               !  Addition of respiration to the release of inorganic nutrients
247               !  -------------------------------------------------------------
248               zgrarem = zgrarem + zbasresi + zrespirc
249               zgraren = zgraren + zbasresn + zrespirc * no3rat3
250               zgrarep = zgrarep + zbasresp + zrespirc * po4rat3
251               zgraref = zgraref + zbasresf + zrespirc * ferat3
252
253               !   Update of the TRA arrays
254               !   ------------------------
255               tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) + zgrarep
256               tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zgraren
257               tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zgradoc
258#if defined key_ligand
259               tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zgradoc * ldocz
260#endif
261               tra(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + zgradon
262               tra(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zgradop
263               tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * zgrarem 
264               tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgraref
265               tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zepsherv * zgraztotc - zrespirc   &
266               &                     - ztortz - zgrazz
267               tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) - zgraznc
268               tra(ji,jj,jk,jpnph) = tra(ji,jj,jk,jpnph) - zgraznn
269               tra(ji,jj,jk,jppph) = tra(ji,jj,jk,jppph) - zgraznp
270               tra(ji,jj,jk,jppic) = tra(ji,jj,jk,jppic) - zgrazpc
271               tra(ji,jj,jk,jpnpi) = tra(ji,jj,jk,jpnpi) - zgrazpn
272               tra(ji,jj,jk,jpppi) = tra(ji,jj,jk,jpppi) - zgrazpp
273               tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazdc
274               tra(ji,jj,jk,jpndi) = tra(ji,jj,jk,jpndi) - zgrazdn
275               tra(ji,jj,jk,jppdi) = tra(ji,jj,jk,jppdi) - zgrazdp
276               tra(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) - zgraznc * trb(ji,jj,jk,jpnch)/(trb(ji,jj,jk,jpphy)+rtrn)
277               tra(ji,jj,jk,jppch) = tra(ji,jj,jk,jppch) - zgrazpc * trb(ji,jj,jk,jppch)/(trb(ji,jj,jk,jppic)+rtrn)
278               tra(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) - zgrazdc * trb(ji,jj,jk,jpdch)/(trb(ji,jj,jk,jpdia)+rtrn)
279               tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zgrazdc * trb(ji,jj,jk,jpdsi)/(trb(ji,jj,jk,jpdia)+rtrn)
280               tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zgrazdc * trb(ji,jj,jk,jpdsi)/(trb(ji,jj,jk,jpdia)+rtrn)
281               tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zgraznf
282               tra(ji,jj,jk,jppfe) = tra(ji,jj,jk,jppfe) - zgrazpf
283               tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) - zgrazdf
284               tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + ztortz + zgrapoc - zgrazpoc 
285               prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + ztortz + zgrapoc
286               conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zgrazpoc
287               tra(ji,jj,jk,jppon) = tra(ji,jj,jk,jppon) + no3rat3 * ztortz + zgrapon - zgrazpon
288               tra(ji,jj,jk,jppop) = tra(ji,jj,jk,jppop) + po4rat3 * ztortz + zgrapop - zgrazpop
289               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ferat3 * ztortz  + zgrapof - zgrazpof
290#if defined key_kriest
291               tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + zgrapoc * xkr_dmicro
292#endif
293               !
294               ! calcite production
295               zprcaca = xfracal(ji,jj,jk) * zgraznc
296               prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca  ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo)
297               !
298               zprcaca = part * zprcaca
299               tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrarem - zprcaca
300               tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * zprcaca     &
301               &                     + rno3 * zgraren
302               tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) + zprcaca
303#if defined key_kriest
304               tra(ji,jj,jk,jpnum) = tra(ji,jj,jk,jpnum) + ztortz * xkr_dmicro &
305                                                         - zgrazpoc * trb(ji,jj,jk,jpnum) / ( trb(ji,jj,jk,jppoc) + rtrn )
306#endif
307            END DO
308         END DO
309      END DO
310      !
311      IF( lk_iomput .AND. knt == nrdttrc ) THEN
312         CALL wrk_alloc( jpi, jpj, jpk, zw3d )
313         IF( iom_use( "GRAZ1" ) ) THEN
314            zw3d(:,:,:) = zgrazing(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:)  !  Total grazing of phyto by zooplankton
315            CALL iom_put( "GRAZ1", zw3d )
316         ENDIF
317         CALL wrk_dealloc( jpi, jpj, jpk, zw3d )
318      ENDIF
319      !
320      IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
321         WRITE(charout, FMT="('micro')")
322         CALL prt_ctl_trc_info(charout)
323         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
324      ENDIF
325      !
326      IF( lk_iomput )  CALL wrk_dealloc( jpi, jpj, jpk, zgrazing )
327      !
328      IF( nn_timing == 1 )  CALL timing_stop('p5z_micro')
329      !
330   END SUBROUTINE p5z_micro
331
332
333   SUBROUTINE p5z_micro_init
334
335      !!----------------------------------------------------------------------
336      !!                  ***  ROUTINE p5z_micro_init  ***
337      !!
338      !! ** Purpose :   Initialization of microzooplankton parameters
339      !!
340      !! ** Method  :   Read the nampiszoo namelist and check the parameters
341      !!                called at the first timestep (nittrc000)
342      !!
343      !! ** input   :   Namelist nampiszoo
344      !!
345      !!----------------------------------------------------------------------
346
347      NAMELIST/nampiszoo/ part, grazrat, bmetexc, resrat, mzrat, xprefc, xprefn, &
348         &                xprefp, xprefd, xprefz, xthreshdia, xthreshphy, &
349         &                xthreshpic, xthreshpoc, xthreshzoo, xthresh, xkgraz, &
350         &                epsher, ssigma, srespir, unassc, unassn, unassp
351
352      INTEGER :: ios                 ! Local integer output status for namelist read
353
354      REWIND( numnatp_ref )              ! Namelist nampiszoo in reference namelist : Pisces microzooplankton
355      READ  ( numnatp_ref, nampiszoo, IOSTAT = ios, ERR = 901)
356901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiszoo in reference namelist', lwp )
357
358      REWIND( numnatp_cfg )              ! Namelist nampiszoo in configuration namelist : Pisces microzooplankton
359      READ  ( numnatp_cfg, nampiszoo, IOSTAT = ios, ERR = 902 )
360902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampiszoo in configuration namelist', lwp )
361      IF(lwm) WRITE ( numonp, nampiszoo )
362
363      IF(lwp) THEN                         ! control print
364         WRITE(numout,*) ' '
365         WRITE(numout,*) ' Namelist parameters for microzooplankton, nampiszoo'
366         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~'
367         WRITE(numout,*) '    part of calcite not dissolved in microzoo guts  part        =', part
368         WRITE(numout,*) '    microzoo preference for POC                     xprefc     =', xprefc
369         WRITE(numout,*) '    microzoo preference for nano                    xprefn     =', xprefn
370         WRITE(numout,*) '    microzoo preference for pico                    xprefp     =', xprefp
371         WRITE(numout,*) '    microzoo preference for diatoms                 xprefd     =', xprefd
372         WRITE(numout,*) '    microzoo preference for microzoo                xprefz     =', xprefz
373         WRITE(numout,*) '    diatoms feeding threshold  for microzoo         xthreshdia  =', xthreshdia
374         WRITE(numout,*) '    nanophyto feeding threshold for microzoo        xthreshphy  =', xthreshphy
375         WRITE(numout,*) '    picophyto feeding threshold for microzoo        xthreshpic  =', xthreshpic
376         WRITE(numout,*) '    poc feeding threshold for microzoo              xthreshpoc  =', xthreshpoc
377         WRITE(numout,*) '    microzoo feeding threshold for microzoo         xthreshzoo  =', xthreshzoo
378         WRITE(numout,*) '    feeding threshold for microzooplankton          xthresh     =', xthresh
379         WRITE(numout,*) '    exsudation rate of microzooplankton             resrat      =', resrat
380         WRITE(numout,*) '    microzooplankton mortality rate                 mzrat       =', mzrat
381         WRITE(numout,*) '    maximal microzoo grazing rate                   grazrat     =', grazrat
382         WRITE(numout,*) '    C egested fraction of fodd by microzoo          unassc      =', unassc
383         WRITE(numout,*) '    N egested fraction of fodd by microzoo          unassn      =', unassn
384         WRITE(numout,*) '    P egested fraction of fodd by microzoo          unassp      =', unassp
385         WRITE(numout,*) '    Efficicency of microzoo growth                  epsher      =', epsher
386         WRITE(numout,*) '    Fraction excreted as semi-labile DOM            ssigma      =', ssigma
387         WRITE(numout,*) '    Active respiration                              srespir     =', srespir
388         WRITE(numout,*) '    half sturation constant for grazing 1           xkgraz      =', xkgraz
389         WRITE(numout,*) '    Use of excess carbon for respiration            bmetexc     =', bmetexc
390      ENDIF
391
392   END SUBROUTINE p5z_micro_init
393
394#else
395   !!======================================================================
396   !!  Dummy module :                                   No PISCES bio-model
397   !!======================================================================
398CONTAINS
399   SUBROUTINE p5z_micro                    ! Empty routine
400   END SUBROUTINE p5z_micro
401#endif 
402
403   !!======================================================================
404END MODULE  p5zmicro
Note: See TracBrowser for help on using the repository browser.