Changeset 5288 for branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmeso.F90

Timestamp:

2015-05-22T18:22:18+02:00 (9 years ago)

Author:

aumont

Message:

various bug fixes and updates of PISCES quota

File:

• branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmeso.F90 (modified) (8 diffs)

branches/CNRS/dev_r4826_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p5zmeso.F90

@@ -7 +7 @@
    !!             2.0  !  2007-12  (C. Ethe, G. Madec)  F90
    !!             3.4  !  2011-06  (O. Aumont, C. Ethe) Quota model for iron
+   !!             3.6  !  2015-05  (O. Aumont) PISCES quota
    !!----------------------------------------------------------------------
 #if defined key_pisces_quota
@@ -47 +48 @@
    REAL(wp), PUBLIC ::  mzrat2       !: microzooplankton mortality rate 
    REAL(wp), PUBLIC ::  grazrat2     !: maximal mesozoo grazing rate
-   REAL(wp), PUBLIC ::  xkgraz2      !: non assimilated fraction of P by mesozoo 
-   REAL(wp), PUBLIC ::  unass2c       !: Efficicency of mesozoo growth 
-   REAL(wp), PUBLIC ::  unass2n       !: Efficicency of mesozoo growth 
-   REAL(wp), PUBLIC ::  unass2p       !: Efficicency of mesozoo growth 
+   REAL(wp), PUBLIC ::  xkgraz2      !: non assimilated fraction of P by mesozoo
+   REAL(wp), PUBLIC ::  unass2c      !: Efficicency of mesozoo growth 
+   REAL(wp), PUBLIC ::  unass2n      !: Efficicency of mesozoo growth 
+   REAL(wp), PUBLIC ::  unass2p      !: Efficicency of mesozoo growth 
    REAL(wp), PUBLIC ::  epsher2      !: half sturation constant for grazing 2
+   REAL(wp), PUBLIC ::  ssigma2      !: Fraction excreted as semi-labile DOM
+   REAL(wp), PUBLIC ::  srespir2     !: Active respiration
    REAL(wp), PUBLIC ::  grazflux     !: mesozoo flux feeding rate
 
@@ -79 +82 @@
       REAL(wp) :: zepsherf, zepshert, zepsherv, zrespirc, zrespirn, zrespirp, zbasresb, zbasresi
       REAL(wp) :: zgraztotc, zgraztotn, zgraztotp, zgraztotf
-      REAL(wp) :: zexcess, zgradoc, zgradon, zgradop
+      REAL(wp) :: zgradoc, zgradon, zgradop, zgratmp, zgradoct, zgradont, zgrafert, zgradopt
       REAL(wp) :: zgrapoc, zgrapon, zgrapop, zgrapof, zprcaca, zmortz
-      REAL(wp) :: zgrarem, zgraren, zgrarep, zgrafer
+      REAL(wp) :: zexcess, zgrarem, zgraren, zgrarep, zgrafer
       REAL(wp) :: zbeta, zrespz, ztortz, zgrasratp, zgrasratn, zgrasratf
       REAL(wp) :: ztmp1, ztmp2, ztmp3, ztmp4, ztmp5, ztmptot
@@ -265 +268 @@
                zbasresb = MAX(0., zrespz - zexcess)
                zbasresi = zexcess + MIN(0., zrespz - zexcess)
-               zrespirc = 0.2 * zepsherv * zgraztotc + zbasresb
+               zrespirc = srespir2 * zepsherv * zgraztotc + zbasresb
 
                !   Voiding of the excessive elements as organic matter
                !   --------------------------------------------------------
-               zgradoc = (1. - unass2c - zepsherv) * zgraztotc - zbasresi
-               zgradon = (1. - unass2n) * zgraztotn - zepsherv * no3rat3 * zgraztotc
-               zgradop = (1. - unass2p) * zgraztotp - zepsherv * po4rat3 * zgraztotc
-               zgrafer = (1. - unass2c) * zgraztotf - zepsherv * ferat3 * zgraztotc
+               zgradoct = (1. - unass2c - zepsherv) * zgraztotc - zbasresi
+               zgradont = (1. - unass2n) * zgraztotn - zepsherv * no3rat3 * zgraztotc
+               zgradopt = (1. - unass2p) * zgraztotp - zepsherv * po4rat3 * zgraztotc
+               zgrafert = (1. - unass2c) * zgraztotf - zepsherv * ferat3 * zgraztotc
                ztmp1   = ( 1. - epsher2 - unass2c ) /( 1. - 0.8 * epsher2 ) * ztortz
-               zgradoc = zgradoc + ztmp1
-               zgradon = zgradon + no3rat3 * ztmp1
-               zgradop = zgradop + po4rat3 * ztmp1
-               zgrafer = zgrafer + ferat3  * ztmp1
-               zgrarem = 0.2 * epsher2 /( 1. - 0.8 * epsher2 ) * ztortz
-               zgraren = no3rat3 * zgrarem
-               zgrarep = po4rat3 * zgrarem
-               zgrafer = zgrafer + ferat3 * zgrarem
+               zgradoc = (zgradoct + ztmp1) * ssigma2
+               zgradon = (zgradont + no3rat3 * ztmp1) * ssigma2
+               zgradop = (zgradopt + po4rat3 * ztmp1) * ssigma2
+               zgratmp = 0.2 * epsher2 /( 1. - 0.8 * epsher2 ) * ztortz
+               zgrarem = zgratmp + ( zgradoct + ztmp1 ) * (1.0 - ssigma2)
+               zgraren = no3rat3 * zgratmp + ( zgradont + no3rat3 * ztmp1 ) * (1.0 - ssigma2)
+               zgrarep = po4rat3 * zgratmp + ( zgradopt + po4rat3 * ztmp1 ) * (1.0 - ssigma2)
+               zgrafer = zgrafert + ferat3 * ( ztmp1 + zgratmp )
 
                !   Defecation as a result of non assimilated products
@@ -302 +305 @@
                tra(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + zgradon
                tra(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + zgradop
-               tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * (zgraren + zrespirn)
+               tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * (zgrarem + zrespirc)
                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgrafer
-               tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zrespirc
+               tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zrespirc + zgrarem
                tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * (zgraren + zrespirn)
                tra(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) + zepsherv * zgraztotc - zrespirc   &
@@ -349 +352 @@
                zprcaca = part2 * zprcaca
                tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrazcal - zprcaca
-               tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2. * ( zgrazcal + zprcaca )
+               tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + 2. * ( zgrazcal - zprcaca )
                tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zgrazcal + zprcaca
             END DO
@@ -389 +392 @@
       NAMELIST/nampismes/ part2, grazrat2, resrat2, mzrat2, xpref2c, xpref2p, xpref2z, xpref2m,   &
          &                xpref2d, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, &
-         &                xthresh2mes, xthresh2, xkgraz2, epsher2, unass2c, &
-      &                   unass2n, unass2p, grazflux
+         &                xthresh2mes, xthresh2, xkgraz2, epsher2, ssigma2, unass2c, &
+      &                   unass2n, unass2p, srespir2, grazflux
 
       INTEGER :: ios                 ! Local integer output status for namelist read
@@ -428 +431 @@
          WRITE(numout,*) '    C egested fraction of food by mesozoo          unass2p       =', unass2p
          WRITE(numout,*) '    Efficicency of Mesozoo growth                  epsher2      =', epsher2
-         WRITE(numout,*) '    half sturation constant for grazing 2          xkgraz2      =', xkgraz2
+         WRITE(numout,*) '    Fraction excreted as semi-labile DOM           ssigma2       =', ssigma2
+         WRITE(numout,*) '    Active respiration                             srespir2      =', srespir2
+         WRITE(numout,*) '    half sturation constant for grazing 2          xkgraz2       =', xkgraz2
       ENDIF