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.
p2zbio.F90 in branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z – NEMO

source: branches/2013/dev_r3853_CNRS9_ConfSetting/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90 @ 3875

Last change on this file since 3875 was 3875, checked in by clevy, 11 years ago

Configuration Setting/Step? 1, see ticket:#1074
File size: 29.3 KB
Line 
1MODULE p2zbio
2   !!======================================================================
3   !!                         ***  MODULE p2zbio  ***
4   !! TOP :   LOBSTER
5   !!======================================================================
6   !! History :    -   !  1999-07  (M. Levy) Original code
7   !!              -   !  2000-12  (E. Kestenare) assign a parameter to name individual tracers
8   !!              -   !  2001-03  (M. Levy)  LNO3 + dia2d
9   !!             2.0  !  2007-12  (C. Deltel, G. Madec)  F90
10   !!----------------------------------------------------------------------
11#if defined key_pisces_reduced
12   !!----------------------------------------------------------------------
13   !!   'key_pisces_reduced'                                     LOBSTER bio-model
14   !!----------------------------------------------------------------------
15   !!   p2z_bio        : 
16   !!----------------------------------------------------------------------
17   USE oce_trc         !
18   USE trc             !
19   USE sms_pisces
20   USE p2zopt
21   USE lbclnk          !
22   USE prtctl_trc      ! Print control for debbuging
23   USE trdmod_oce
24   USE trdmod_trc
25   USE iom
26   
27   IMPLICIT NONE
28   PRIVATE
29
30   PUBLIC   p2z_bio         ! called in ???
31   PUBLIC   p2z_bio_init    ! called in ???
32
33   REAL(wp) ::   tmumax     = 1.21e-5    ! maximal phytoplankton growth rate            [s-1]
34   REAL(wp) ::   rgamma     =  0.05      ! phytoplankton exudation fraction             [%]
35   REAL(wp) ::   fphylab    =  0.75      ! NH4 fraction of phytoplankton exsudation
36   REAL(wp) ::   tmminp     =  5.8e-7    ! minimal phytoplancton mortality rate         [0.05/86400 s-1=20 days]
37   REAL(wp) ::   aki        =  33.       ! light photosynthesis half saturation constant[W/m2]
38   !
39   REAL(wp) ::   akno3      =  0.7       ! nitrate limitation half-saturation value     [mmol/m3]
40   REAL(wp) ::   aknh4      =  0.001     ! ammonium limitation half-saturation value    [mmol/m3]
41   REAL(wp) ::   taunn      =  5.80e-7   ! nitrification rate                           [s-1]
42   REAL(wp) ::   psinut     =  3.        ! inhibition of nitrate uptake by ammonium
43   !
44   REAL(wp) ::   taudn      = 5.80e-7    ! detritus breakdown rate                        [0.1/86400 s-1=10 days]
45   REAL(wp) ::   fdetlab    = 0.         ! NH4 fraction of detritus dissolution
46   !
47   REAL(wp) ::   taudomn    = 6.43e-8    ! DOM breakdown rate                             [s-1]
48   !                                     ! slow remineralization rate of semi-labile dom to nh4 (1 month)
49   !
50   REAL(wp) ::   rppz       = 0.         ! ivlev coeff for zoo mortality
51   REAL(wp) ::   taus       = 9.26E-6    ! specific zooplankton maximal grazing rate              [s-1]
52   !                                     ! 0.75/86400 s-1=8.680555E-6    1/86400 = 1.15e-5
53   REAL(wp) ::   aks        = 1.         ! half-saturation constant for total zooplankton grazing [mmolN.m-3]
54   REAL(wp) ::   rpnaz      = 0.3        ! non-assimilated phytoplankton by zooplancton           [%]
55   REAL(wp) ::   rdnaz      = 0.3        ! non-assimilated detritus by zooplankton                [%]
56   REAL(wp) ::   tauzn      = 8.1e-7     ! zooplancton specific excretion rate                    [0.1/86400 s-1=10 days]
57   REAL(wp) ::   tmminz     = 2.31e-6    ! minimal zooplankton mortality rate                     [(mmolN/m3)-1 d-1]
58   REAL(wp) ::   fzoolab    = 0.5        ! NH4 fraction of zooplankton excretion
59   REAL(wp) ::   fdbod      = 0.5        ! zooplankton mortality fraction that goes to detritus
60
61   !!* Substitution
62#  include "top_substitute.h90"
63   !!----------------------------------------------------------------------
64   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
65   !! $Id: p2zbio.F90 3294 2012-01-28 16:44:18Z rblod $
66   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
67   !!----------------------------------------------------------------------
68
69CONTAINS
70
71   SUBROUTINE p2z_bio( kt )
72      !!---------------------------------------------------------------------
73      !!                     ***  ROUTINE p2z_bio  ***
74      !!
75      !! ** Purpose :   compute the now trend due to biogeochemical processes
76      !!              and add it to the general trend of passive tracers equations
77      !!
78      !! ** Method  :   each now biological flux is calculated in function of now
79      !!              concentrations of tracers.
80      !!              depending on the tracer, these fluxes are sources or sinks.
81      !!              the total of the sources and sinks for each tracer
82      !!              is added to the general trend.
83      !!       
84      !!                      tra = tra + zf...tra - zftra...
85      !!                                     |         |
86      !!                                     |         |
87      !!                                  source      sink
88      !!       
89      !!              IF 'key_diabio' defined , the biogeochemical trends
90      !!              for passive tracers are saved for futher diagnostics.
91      !!---------------------------------------------------------------------
92      !!
93      INTEGER, INTENT( in ) ::   kt      ! ocean time-step index     
94      !!
95      INTEGER  ::   ji, jj, jk, jl
96      REAL(wp) ::   zdet, zzoo, zphy, zno3, znh4, zdom      ! now concentrations
97      REAL(wp) ::   zlno3, zlnh4, zle, zlt                  ! limitation terms for phyto
98      REAL(wp) ::   zno3phy, znh4phy, zphynh4, zphydom
99      REAL(wp) ::   zphydet, zphyzoo, zdetzoo
100      REAL(wp) ::   zzoonh4, zzoodom, zzoodet, zdetnh4, zdetdom
101      REAL(wp) ::   znh4no3, zdomnh4, zppz, zpdz, zpppz, zppdz, zfood
102      REAL(wp) ::   zfilpz, zfildz, zphya, zzooa, zno3a
103      REAL(wp) ::   znh4a, zdeta, zdoma, zzoobod, zboddet, zdomaju
104      REAL(wp) ::   ze3t
105      REAL(wp), POINTER,   DIMENSION(:,:,:) :: zw2d
106      REAL(wp), POINTER, DIMENSION(:,:,:,:) :: zw3d
107      CHARACTER (len=25) :: charout
108      !!---------------------------------------------------------------------
109      !
110      IF( nn_timing == 1 )  CALL timing_start('p2z_bio')
111      !
112      IF( ln_diatrc ) THEN
113         CALL wrk_alloc( jpi, jpj,     17, zw2d )
114         CALL wrk_alloc( jpi, jpj, jpk, 3, zw3d )
115      ENDIF
116
117      IF( kt == nittrc000 ) THEN
118         IF(lwp) WRITE(numout,*)
119         IF(lwp) WRITE(numout,*) ' p2z_bio: LOBSTER bio-model'
120         IF(lwp) WRITE(numout,*) ' ~~~~~~~'
121      ENDIF
122
123      xksi(:,:) = 0.e0        ! zooplakton closure ( fbod)
124      IF( ln_diatrc ) THEN
125         zw2d  (:,:,:) = 0.e0
126         zw3d(:,:,:,:) = 0.e0
127      ENDIF
128
129      !                                      ! -------------------------- !
130      DO jk = 1, jpkbm1                      !  Upper ocean (bio-layers)  !
131         !                                   ! -------------------------- !
132         DO jj = 2, jpjm1
133            DO ji = fs_2, fs_jpim1 
134               ! trophic variables( det, zoo, phy, no3, nh4, dom)
135               ! ------------------------------------------------
136
137               ! negative trophic variables DO not contribute to the fluxes
138               zdet = MAX( 0.e0, trn(ji,jj,jk,jpdet) )
139               zzoo = MAX( 0.e0, trn(ji,jj,jk,jpzoo) )
140               zphy = MAX( 0.e0, trn(ji,jj,jk,jpphy) )
141               zno3 = MAX( 0.e0, trn(ji,jj,jk,jpno3) )
142               znh4 = MAX( 0.e0, trn(ji,jj,jk,jpnh4) )
143               zdom = MAX( 0.e0, trn(ji,jj,jk,jpdom) )
144
145               ! Limitations
146               zlt   = 1.
147               zle   = 1. - EXP( -etot(ji,jj,jk) / aki / zlt )
148               ! psinut,akno3,aknh4 added by asklod AS Kremeur 2005-03
149               zlno3 = zno3 * EXP( -psinut * znh4 ) / ( akno3 + zno3 )
150               zlnh4 = znh4 / (znh4+aknh4) 
151
152               ! sinks and sources
153               !    phytoplankton production and exsudation
154               zno3phy = tmumax * zle * zlt * zlno3 * zphy
155               znh4phy = tmumax * zle * zlt * zlnh4 * zphy
156
157               !    fphylab added by asklod AS Kremeur 2005-03
158               zphydom = rgamma * (1 - fphylab) * (zno3phy + znh4phy)
159               zphynh4 = rgamma * fphylab * (zno3phy + znh4phy)
160               ! zooplankton production
161               !    preferences
162               zppz = rppz
163               zpdz = 1. - rppz
164               zpppz = ( zppz * zphy ) / ( ( zppz * zphy + zpdz * zdet ) + 1.e-13 )
165               zppdz = ( zpdz * zdet ) / ( ( zppz * zphy + zpdz * zdet ) + 1.e-13 )
166               zfood = zpppz * zphy + zppdz * zdet
167               !    filtration
168               zfilpz = taus * zpppz / (aks + zfood)
169               zfildz = taus * zppdz / (aks + zfood)
170               !    grazing
171               zphyzoo = zfilpz * zphy * zzoo
172               zdetzoo = zfildz * zdet * zzoo
173
174               ! fecal pellets production
175               zzoodet = rpnaz * zphyzoo + rdnaz * zdetzoo
176
177               ! zooplankton liquide excretion
178               zzoonh4 = tauzn * fzoolab * zzoo 
179               zzoodom = tauzn * (1 - fzoolab) * zzoo
180
181               ! mortality
182               !    phytoplankton mortality
183               zphydet = tmminp * zphy
184
185               !    zooplankton mortality
186               !    closure : flux grazing is redistributed below level jpkbio
187               zzoobod = tmminz * zzoo * zzoo
188               xksi(ji,jj) = xksi(ji,jj) + (1-fdbod) * zzoobod * fse3t(ji,jj,jk)
189               zboddet = fdbod * zzoobod
190
191               ! detritus and dom breakdown
192               zdetnh4 = taudn * fdetlab * zdet
193               zdetdom = taudn * (1 - fdetlab) * zdet
194
195               zdomnh4 = taudomn * zdom
196
197               ! flux added to express how the excess of nitrogen from
198               ! PHY, ZOO and DET to DOM goes directly to NH4 (flux of ajustment)
199               zdomaju = (1 - redf/reddom) * (zphydom + zzoodom + zdetdom)
200
201               ! Nitrification
202               znh4no3 = taunn * znh4
203
204               ! determination of trends
205               !    total trend for each biological tracer
206               zphya =   zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo - zphydet
207               zzooa =   zphyzoo + zdetzoo - zzoodet - zzoodom - zzoonh4 - zzoobod
208               zno3a = - zno3phy + znh4no3
209               znh4a = - znh4phy - znh4no3 + zphynh4 + zzoonh4 + zdomnh4 + zdetnh4 + zdomaju
210               zdeta =   zphydet + zzoodet - zdetzoo - zdetnh4 - zdetdom + zboddet
211               zdoma =   zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju
212
213               ! tracer flux at totox-point added to the general trend
214               tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta
215               tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa
216               tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya
217               tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a
218               tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a
219               tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma
220
221
222               IF( ( ln_diabio .AND. .NOT. lk_iomput ) .OR. l_trdtrc ) THEN
223                  trbio(ji,jj,jk,jp_pcs0_trd     ) = zno3phy
224                  trbio(ji,jj,jk,jp_pcs0_trd +  1) = znh4phy
225                  trbio(ji,jj,jk,jp_pcs0_trd +  2) = zphynh4
226                  trbio(ji,jj,jk,jp_pcs0_trd +  3) = zphydom
227                  trbio(ji,jj,jk,jp_pcs0_trd +  4) = zphyzoo
228                  trbio(ji,jj,jk,jp_pcs0_trd +  5) = zphydet
229                  trbio(ji,jj,jk,jp_pcs0_trd +  6) = zdetzoo
230                  !  trend number 8 in p2zsed
231                  trbio(ji,jj,jk,jp_pcs0_trd +  8) = zzoodet
232                  trbio(ji,jj,jk,jp_pcs0_trd +  9) = zzoobod
233                  trbio(ji,jj,jk,jp_pcs0_trd + 10) = zzoonh4
234                  trbio(ji,jj,jk,jp_pcs0_trd + 11) = zzoodom
235                  trbio(ji,jj,jk,jp_pcs0_trd + 12) = znh4no3
236                  trbio(ji,jj,jk,jp_pcs0_trd + 13) = zdomnh4
237                  trbio(ji,jj,jk,jp_pcs0_trd + 14) = zdetnh4
238                  trbio(ji,jj,jk,jp_pcs0_trd + 15) = zdetdom
239                  !  trend number 17 in p2zexp
240                ENDIF
241                IF( ln_diatrc ) THEN
242                  ! convert fluxes in per day
243                  ze3t = fse3t(ji,jj,jk) * 86400.
244                  zw2d(ji,jj,1)  = zw2d(ji,jj,1)  + zno3phy * ze3t
245                  zw2d(ji,jj,2)  = zw2d(ji,jj,2)  + znh4phy * ze3t
246                  zw2d(ji,jj,3)  = zw2d(ji,jj,3)  + zphydom * ze3t
247                  zw2d(ji,jj,4)  = zw2d(ji,jj,4)  + zphynh4 * ze3t
248                  zw2d(ji,jj,5)  = zw2d(ji,jj,5)  + zphyzoo * ze3t
249                  zw2d(ji,jj,6)  = zw2d(ji,jj,6)  + zphydet * ze3t
250                  zw2d(ji,jj,7)  = zw2d(ji,jj,7)  + zdetzoo * ze3t
251                  zw2d(ji,jj,8)  = zw2d(ji,jj,8)  + zzoodet * ze3t
252                  zw2d(ji,jj,9)  = zw2d(ji,jj,9)  + zzoobod * ze3t
253                  zw2d(ji,jj,10) = zw2d(ji,jj,10) + zzoonh4 * ze3t
254                  zw2d(ji,jj,11) = zw2d(ji,jj,11) + zzoodom * ze3t
255                  zw2d(ji,jj,12) = zw2d(ji,jj,12) + znh4no3 * ze3t
256                  zw2d(ji,jj,13) = zw2d(ji,jj,13) + zdomnh4 * ze3t
257                  zw2d(ji,jj,14) = zw2d(ji,jj,14) + zdetnh4 * ze3t
258                  zw2d(ji,jj,15) = zw2d(ji,jj,15) + ( zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo - zphydet ) * ze3t
259                  zw2d(ji,jj,16) = zw2d(ji,jj,16) + ( zphyzoo + zdetzoo - zzoodet - zzoobod - zzoonh4 - zzoodom ) * ze3t
260                  zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t
261                  !   
262                  zw3d(ji,jj,jk,1) = zno3phy * 86400
263                  zw3d(ji,jj,jk,2) = znh4phy * 86400     
264                  zw3d(ji,jj,jk,3) = znh4no3 * 86400   
265                   !
266                ENDIF
267            END DO
268         END DO
269      END DO
270
271      !                                      ! -------------------------- !
272      DO jk = jpkb, jpkm1                    !  Upper ocean (bio-layers)  !
273         !                                   ! -------------------------- !
274         DO jj = 2, jpjm1
275            DO ji = fs_2, fs_jpim1 
276               ! remineralisation of all quantities towards nitrate
277
278               !    trophic variables( det, zoo, phy, no3, nh4, dom)
279               !       negative trophic variables DO not contribute to the fluxes
280               zdet = MAX( 0.e0, trn(ji,jj,jk,jpdet) )
281               zzoo = MAX( 0.e0, trn(ji,jj,jk,jpzoo) )
282               zphy = MAX( 0.e0, trn(ji,jj,jk,jpphy) )
283               zno3 = MAX( 0.e0, trn(ji,jj,jk,jpno3) )
284               znh4 = MAX( 0.e0, trn(ji,jj,jk,jpnh4) )
285               zdom = MAX( 0.e0, trn(ji,jj,jk,jpdom) )
286
287               !    Limitations
288               zlt   = 0.e0
289               zle   = 0.e0
290               zlno3 = 0.e0
291               zlnh4 = 0.e0
292
293               !    sinks and sources
294               !       phytoplankton production and exsudation
295               zno3phy = 0.e0
296               znh4phy = 0.e0
297               zphydom = 0.e0
298               zphynh4 = 0.e0
299
300               !    zooplankton production
301               zphyzoo = 0.e0      ! grazing
302               zdetzoo = 0.e0
303
304               zzoodet = 0.e0      ! fecal pellets production
305
306               zzoonh4 = tauzn * fzoolab * zzoo         ! zooplankton liquide excretion
307               zzoodom = tauzn * (1 - fzoolab) * zzoo
308
309               !    mortality
310               zphydet = tmminp * zphy      ! phytoplankton mortality
311
312               zzoobod = 0.e0               ! zooplankton mortality
313               zboddet = 0.e0               ! closure : flux fbod is redistributed below level jpkbio
314
315               !    detritus and dom breakdown
316               zdetnh4 = taudn * fdetlab * zdet
317               zdetdom = taudn * (1 - fdetlab) * zdet
318
319               zdomnh4 = taudomn * zdom
320               zdomaju = (1 - redf/reddom) * (zphydom + zzoodom + zdetdom)
321
322               !    Nitrification
323               znh4no3 = taunn * znh4
324
325
326               ! determination of trends
327               !     total trend for each biological tracer
328               zphya =   zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo - zphydet
329               zzooa =   zphyzoo + zdetzoo - zzoodet - zzoodom - zzoonh4 - zzoobod
330               zno3a = - zno3phy + znh4no3 
331               znh4a = - znh4phy - znh4no3 + zphynh4 + zzoonh4 + zdomnh4 + zdetnh4 + zdomaju
332               zdeta = zphydet + zzoodet  - zdetzoo - zdetnh4 - zdetdom + zboddet
333               zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju
334
335               ! tracer flux at totox-point added to the general trend
336               tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta
337               tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa
338               tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya
339               tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a
340               tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a
341               tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma
342               !
343               IF( ( ln_diabio .AND. .NOT. lk_iomput ) .OR. l_trdtrc ) THEN
344                  trbio(ji,jj,jk,jp_pcs0_trd     ) = zno3phy
345                  trbio(ji,jj,jk,jp_pcs0_trd +  1) = znh4phy
346                  trbio(ji,jj,jk,jp_pcs0_trd +  2) = zphynh4
347                  trbio(ji,jj,jk,jp_pcs0_trd +  3) = zphydom
348                  trbio(ji,jj,jk,jp_pcs0_trd +  4) = zphyzoo
349                  trbio(ji,jj,jk,jp_pcs0_trd +  5) = zphydet
350                  trbio(ji,jj,jk,jp_pcs0_trd +  6) = zdetzoo
351                  !  trend number 8 in p2zsed
352                  trbio(ji,jj,jk,jp_pcs0_trd +  8) = zzoodet
353                  trbio(ji,jj,jk,jp_pcs0_trd +  9) = zzoobod
354                  trbio(ji,jj,jk,jp_pcs0_trd + 10) = zzoonh4
355                  trbio(ji,jj,jk,jp_pcs0_trd + 11) = zzoodom
356                  trbio(ji,jj,jk,jp_pcs0_trd + 12) = znh4no3
357                  trbio(ji,jj,jk,jp_pcs0_trd + 13) = zdomnh4
358                  trbio(ji,jj,jk,jp_pcs0_trd + 14) = zdetnh4
359                  trbio(ji,jj,jk,jp_pcs0_trd + 15) = zdetdom
360                  !  trend number 17 in p2zexp
361                ENDIF
362                IF( ln_diatrc ) THEN
363                  ! convert fluxes in per day
364                  ze3t = fse3t(ji,jj,jk) * 86400.
365                  zw2d(ji,jj,1)  = zw2d(ji,jj,1)  + zno3phy * ze3t
366                  zw2d(ji,jj,2)  = zw2d(ji,jj,2)  + znh4phy * ze3t
367                  zw2d(ji,jj,3)  = zw2d(ji,jj,3)  + zphydom * ze3t
368                  zw2d(ji,jj,4)  = zw2d(ji,jj,4)  + zphynh4 * ze3t
369                  zw2d(ji,jj,5)  = zw2d(ji,jj,5)  + zphyzoo * ze3t
370                  zw2d(ji,jj,6)  = zw2d(ji,jj,6)  + zphydet * ze3t
371                  zw2d(ji,jj,7)  = zw2d(ji,jj,7)  + zdetzoo * ze3t
372                  zw2d(ji,jj,8)  = zw2d(ji,jj,8)  + zzoodet * ze3t
373                  zw2d(ji,jj,9)  = zw2d(ji,jj,9)  + zzoobod * ze3t
374                  zw2d(ji,jj,10) = zw2d(ji,jj,10) + zzoonh4 * ze3t
375                  zw2d(ji,jj,11) = zw2d(ji,jj,11) + zzoodom * ze3t
376                  zw2d(ji,jj,12) = zw2d(ji,jj,12) + znh4no3 * ze3t
377                  zw2d(ji,jj,13) = zw2d(ji,jj,13) + zdomnh4 * ze3t
378                  zw2d(ji,jj,14) = zw2d(ji,jj,14) + zdetnh4 * ze3t
379                  zw2d(ji,jj,15) = zw2d(ji,jj,15) + ( zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo - zphydet ) * ze3t
380                  zw2d(ji,jj,16) = zw2d(ji,jj,16) + ( zphyzoo + zdetzoo - zzoodet - zzoobod - zzoonh4 - zzoodom ) * ze3t
381                  zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t
382                  !   
383                  zw3d(ji,jj,jk,1) = zno3phy * 86400
384                  zw3d(ji,jj,jk,2) = znh4phy * 86400
385                  zw3d(ji,jj,jk,3) = znh4no3 * 86400
386                   !
387                ENDIF
388            END DO
389         END DO
390      END DO
391
392      IF( ln_diatrc ) THEN
393         !
394         DO jl = 1, 17 
395            CALL lbc_lnk( zw2d(:,:,jl),'T', 1. )
396         END DO
397         DO jl = 1, 3
398            CALL lbc_lnk( zw3d(:,:,:,jl),'T', 1. )
399         END DO
400         IF( lk_iomput ) THEN
401            ! Save diagnostics
402            CALL iom_put( "TNO3PHY", zw2d(:,:,1) )
403            CALL iom_put( "TNH4PHY", zw2d(:,:,2) )
404            CALL iom_put( "TPHYDOM", zw2d(:,:,3) )
405            CALL iom_put( "TPHYNH4", zw2d(:,:,4) )
406            CALL iom_put( "TPHYZOO", zw2d(:,:,5) )
407            CALL iom_put( "TPHYDET", zw2d(:,:,6) )
408            CALL iom_put( "TDETZOO", zw2d(:,:,7) )
409            CALL iom_put( "TZOODET", zw2d(:,:,8) )
410            CALL iom_put( "TZOOBOD", zw2d(:,:,9) )
411            CALL iom_put( "TZOONH4", zw2d(:,:,10) )
412            CALL iom_put( "TZOODOM", zw2d(:,:,11) )
413            CALL iom_put( "TNH4NO3", zw2d(:,:,12) )
414            CALL iom_put( "TDOMNH4", zw2d(:,:,13) )
415            CALL iom_put( "TDETNH4", zw2d(:,:,14) )
416            CALL iom_put( "TPHYTOT", zw2d(:,:,15) )
417            CALL iom_put( "TZOOTOT", zw2d(:,:,16) )
418            !
419            CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) )
420            CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) )
421            CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) )
422            !
423         ELSE
424            !
425            trc2d(:,:,jp_pcs0_2d    ) = zw2d(:,:,1) 
426            trc2d(:,:,jp_pcs0_2d + 1) = zw2d(:,:,2) 
427            trc2d(:,:,jp_pcs0_2d + 2) = zw2d(:,:,3) 
428            trc2d(:,:,jp_pcs0_2d + 3) = zw2d(:,:,4) 
429            trc2d(:,:,jp_pcs0_2d + 4) = zw2d(:,:,5) 
430            trc2d(:,:,jp_pcs0_2d + 5) = zw2d(:,:,6) 
431            trc2d(:,:,jp_pcs0_2d + 6) = zw2d(:,:,7) 
432                     ! trend number 8 is in p2zsed.F
433            trc2d(:,:,jp_pcs0_2d +  8) = zw2d(:,:,8) 
434            trc2d(:,:,jp_pcs0_2d +  9) = zw2d(:,:,9) 
435            trc2d(:,:,jp_pcs0_2d + 10) = zw2d(:,:,10) 
436            trc2d(:,:,jp_pcs0_2d + 11) = zw2d(:,:,11) 
437            trc2d(:,:,jp_pcs0_2d + 12) = zw2d(:,:,12) 
438            trc2d(:,:,jp_pcs0_2d + 13) = zw2d(:,:,13) 
439            trc2d(:,:,jp_pcs0_2d + 14) = zw2d(:,:,14) 
440            trc2d(:,:,jp_pcs0_2d + 15) = zw2d(:,:,15) 
441            trc2d(:,:,jp_pcs0_2d + 16) = zw2d(:,:,16) 
442            trc2d(:,:,jp_pcs0_2d + 17) = zw2d(:,:,17) 
443            ! trend number 19 is in p2zexp.F
444            trc3d(:,:,:,jp_pcs0_3d    ) = zw3d(:,:,:,1) 
445            trc3d(:,:,:,jp_pcs0_3d + 1) = zw3d(:,:,:,2) 
446            trc3d(:,:,:,jp_pcs0_3d + 2) = zw3d(:,:,:,3) 
447         ENDIF
448        !
449      ENDIF
450
451      IF( ln_diabio .AND. .NOT. lk_iomput )  THEN
452         DO jl = jp_pcs0_trd, jp_pcs1_trd
453            CALL lbc_lnk( trbio(:,:,1,jl),'T', 1. )
454         END DO
455      ENDIF
456      !
457      IF( l_trdtrc ) THEN
458         DO jl = jp_pcs0_trd, jp_pcs1_trd
459            CALL trd_mod_trc( trbio(:,:,:,jl), jl, kt )   ! handle the trend
460         END DO
461      ENDIF
462
463      IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
464         WRITE(charout, FMT="('bio')")
465         CALL prt_ctl_trc_info(charout)
466         CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
467      ENDIF
468      !
469      IF( ln_diatrc ) THEN
470         CALL wrk_dealloc( jpi, jpj,     17, zw2d )
471         CALL wrk_dealloc( jpi, jpj, jpk, 3, zw3d )
472      ENDIF
473      !
474      IF( nn_timing == 1 )  CALL timing_stop('p2z_bio')
475      !
476   END SUBROUTINE p2z_bio
477
478   SUBROUTINE p2z_bio_init
479      !!----------------------------------------------------------------------
480      !!                  ***  ROUTINE p2z_bio_init  ***
481      !!
482      !! ** Purpose :  biological parameters
483      !!
484      !! ** Method  :   Read namelist and check the parameters
485      !!
486      !!----------------------------------------------------------------------
487      NAMELIST/namlobphy/ tmumax, rgamma, fphylab, tmminp, aki
488      NAMELIST/namlobnut/ akno3, aknh4, taunn, psinut
489      NAMELIST/namlobzoo/ rppz, taus, aks, rpnaz, rdnaz, tauzn, fzoolab, fdbod, tmminz
490      NAMELIST/namlobdet/  taudn, fdetlab
491      NAMELIST/namlobdom/ taudomn
492      INTEGER :: ios                 ! Local integer output status for namelist read
493      !!----------------------------------------------------------------------
494
495      REWIND( numnatp_ref )              ! Namelist namlobphy in reference namelist : Lobster biological parameters
496      READ  ( numnatp_ref, namlobphy, IOSTAT = ios, ERR = 901)
497901   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobphy in reference namelist', lwp )
498
499      REWIND( numnatp_cfg )              ! Namelist namlobphy in configuration namelist : Lobster biological parameters
500      READ  ( numnatp_cfg, namlobphy, IOSTAT = ios, ERR = 902 )
501902   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobphy in configuration namelist', lwp )
502      WRITE ( numonp, namlobphy )
503
504      IF(lwp) THEN
505          WRITE(numout,*) ' Namelist namlobphy'
506          WRITE(numout,*) '    phyto max growth rate                                tmumax    =', 86400 * tmumax, ' d'
507          WRITE(numout,*) '    phytoplankton exudation fraction                     rgamma    =', rgamma
508          WRITE(numout,*) '    NH4 fraction of phytoplankton exsudation             fphylab   =', fphylab
509          WRITE(numout,*) '    minimal phyto mortality rate                         tmminp    =', 86400 * tmminp
510          WRITE(numout,*) '    light hlaf saturation constant                       aki       =', aki
511          WRITE(numout,*) ' '
512      ENDIF
513
514      REWIND( numnatp_ref )              ! Namelist namlobnut in reference namelist : Lobster nutriments parameters
515      READ  ( numnatp_ref, namlobnut, IOSTAT = ios, ERR = 903)
516903   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobnut in reference namelist', lwp )
517
518      REWIND( numnatp_cfg )              ! Namelist namlobnut in configuration namelist : Lobster nutriments parameters
519      READ  ( numnatp_cfg, namlobnut, IOSTAT = ios, ERR = 904 )
520904   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobnut in configuration namelist', lwp )
521      WRITE ( numonp, namlobnut )
522
523      IF(lwp) THEN
524          WRITE(numout,*) ' Namelist namlobnut'
525          WRITE(numout,*) '    half-saturation nutrient for no3 uptake              akno3     =', akno3
526          WRITE(numout,*) '    half-saturation nutrient for nh4 uptake              aknh4     =', aknh4
527          WRITE(numout,*) '    nitrification rate                                   taunn     =', taunn
528          WRITE(numout,*) '    inhibition of no3 uptake by nh4                      psinut    =', psinut
529          WRITE(numout,*) ' '
530      ENDIF
531
532      REWIND( numnatp_ref )              ! Namelist namlobzoo in reference namelist : Lobster zooplankton parameters
533      READ  ( numnatp_ref, namlobzoo, IOSTAT = ios, ERR = 905)
534905   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobzoo in reference namelist', lwp )
535
536      REWIND( numnatp_cfg )              ! Namelist namlobzoo in configuration namelist : Lobster zooplankton parameters
537      READ  ( numnatp_cfg, namlobzoo, IOSTAT = ios, ERR = 906 )
538906   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobzoo in configuration namelist', lwp )
539      WRITE ( numonp, namlobzoo )
540
541      IF(lwp) THEN
542          WRITE(numout,*) ' Namelist namlobzoo'
543          WRITE(numout,*) '    zoo preference for phyto                             rppz      =', rppz
544          WRITE(numout,*) '    maximal zoo grazing rate                             taus      =', 86400 * taus, ' d'
545          WRITE(numout,*) '    half saturation constant for zoo food                aks       =', aks
546          WRITE(numout,*) '    non-assimilated phyto by zoo                         rpnaz     =', rpnaz
547          WRITE(numout,*) '    non-assimilated detritus by zoo                      rdnaz     =', rdnaz
548          WRITE(numout,*) '    zoo specific excretion rate                          tauzn     =', 86400 * tauzn
549          WRITE(numout,*) '    minimal zoo mortality rate                           tmminz    =', 86400 * tmminz
550          WRITE(numout,*) '    NH4 fraction of zooplankton excretion                fzoolab   =', fzoolab
551          WRITE(numout,*) '    Zooplankton mortality fraction that goes to detritus fdbod     =', fdbod
552          WRITE(numout,*) ' '
553      ENDIF
554
555      REWIND( numnatp_ref )              ! Namelist namlobdet in reference namelist : Lobster detritus parameters
556      READ  ( numnatp_ref, namlobdet, IOSTAT = ios, ERR = 907)
557907   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdet in reference namelist', lwp )
558
559      REWIND( numnatp_cfg )              ! Namelist namlobdet in configuration namelist : Lobster detritus parameters
560      READ  ( numnatp_cfg, namlobdet, IOSTAT = ios, ERR = 908 )
561908   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdet in configuration namelist', lwp )
562      WRITE ( numonp, namlobdet )
563
564      IF(lwp) THEN
565          WRITE(numout,*) ' Namelist namlobdet'
566          WRITE(numout,*) '    detrital breakdown rate                              taudn     =', 86400 * taudn , ' d'
567          WRITE(numout,*) '    NH4 fraction of detritus dissolution                 fdetlab   =', fdetlab
568          WRITE(numout,*) ' '
569      ENDIF
570
571      REWIND( numnatp_ref )              ! Namelist namlobdom in reference namelist : Lobster DOM breakdown rate
572      READ  ( numnatp_ref, namlobdom, IOSTAT = ios, ERR = 909)
573909   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdom in reference namelist', lwp )
574
575      REWIND( numnatp_cfg )              ! Namelist namlobdom in configuration namelist : Lobster DOM breakdown rate
576      READ  ( numnatp_cfg, namlobdom, IOSTAT = ios, ERR = 910 )
577910   IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdom in configuration namelist', lwp )
578      WRITE ( numonp, namlobdom )
579
580      IF(lwp) THEN
581          WRITE(numout,*) ' Namelist namlobdom'
582          WRITE(numout,*) '    DOM breakdown rate                                 taudomn     =', 86400 * taudn , ' d'
583          WRITE(numout,*) ' '
584      ENDIF
585      !
586   END SUBROUTINE p2z_bio_init
587
588#else
589   !!======================================================================
590   !!  Dummy module :                                   No PISCES bio-model
591   !!======================================================================
592CONTAINS
593   SUBROUTINE p2z_bio( kt )                   ! Empty routine
594      INTEGER, INTENT( in ) ::   kt
595      WRITE(*,*) 'p2z_bio: You should not have seen this print! error?', kt
596   END SUBROUTINE p2z_bio
597#endif 
598
599   !!======================================================================
600END MODULE  p2zbio
Note: See TracBrowser for help on using the repository browser.