Changeset 3446 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

Timestamp:

2012-08-10T13:13:55+02:00 (12 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB : 2nd step new PISCES updates from Olivier, see ticket #972

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90 (modified) (5 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zmort.F90

@@ -29 +29 @@
    REAL(wp), PUBLIC :: wchl   = 0.001_wp  !:
    REAL(wp), PUBLIC :: wchld  = 0.02_wp   !:
+   REAL(wp), PUBLIC :: wchldm = 0.05_wp   !:
    REAL(wp), PUBLIC :: mprat  = 0.01_wp   !:
    REAL(wp), PUBLIC :: mprat2 = 0.01_wp   !:
@@ -150 +151 @@
       REAL(wp) ::  zfactfe,zfactsi,zfactch, zcompadi
       REAL(wp) ::  zrespp2, ztortp2, zmortp2, zstep
+      REAL(wp) ::  zlim2, zlim1
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -177 +179 @@
                !  Phytoplankton respiration 
                !     ------------------------
-!               zrespp2  = 1.e6 * zstep * (  wchl + wchld * max(0.,2. * ( 0.5 - xlimdia(ji,jj,jk) ) )  )    &
-               zrespp2  = 1.e6 * zstep * (  wchl + wchld * ( 1.0 - xlimdia(ji,jj,jk) ) )      &
-                  &                   * xdiss(ji,jj,jk) * zcompadi * trn(ji,jj,jk,jpdia)
+               zlim2   = xlimdia(ji,jj,jk) * xlimdia(ji,jj,jk)
+               zlim1   = 0.25 * ( 1. - zlim2 ) / ( 0.25 + zlim2 ) 
+               zrespp2 = 1.e6 * zstep * (  wchld + wchldm * zlim1 ) * xdiss(ji,jj,jk) * zcompadi * trn(ji,jj,jk,jpdia)
+              ! zlim1   = 1.0 - xlimdia(ji,jj,jk) 
+              ! zrespp2 = 1.e6 * zstep * (  wchl + wchld * zlim1 ) * xdiss(ji,jj,jk) * zcompadi * trn(ji,jj,jk,jpdia)
 
                !     Phytoplankton mortality. 
                !     ------------------------
-               ztortp2  = mprat2 * zstep * trn(ji,jj,jk,jpdia)  / ( xkmort + trn(ji,jj,jk,jpdia) ) * zcompadi 
+               ztortp2 = mprat2 * zstep * trn(ji,jj,jk,jpdia)  / ( xkmort + trn(ji,jj,jk,jpdia) ) * zcompadi 
 
                zmortp2 = zrespp2 + ztortp2
@@ -235 +239 @@
       !!----------------------------------------------------------------------
 
-      NAMELIST/nampismort/ wchl, wchld, mprat, mprat2, mpratm
+      NAMELIST/nampismort/ wchl, wchld, wchldm, mprat, mprat2, mpratm
 
       REWIND( numnatp )                     ! read numnatp
@@ -246 +250 @@
          WRITE(numout,*) '    quadratic mortality of phytoplankton      wchl      =', wchl
          WRITE(numout,*) '    maximum quadratic mortality of diatoms    wchld     =', wchld
+         WRITE(numout,*) '    maximum quadratic mortality of diatoms    wchld     =', wchldm
          WRITE(numout,*) '    phytoplankton mortality rate              mprat     =', mprat
          WRITE(numout,*) '    Diatoms mortality rate                    mprat2    =', mprat2