Changeset 14276 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmort.F90

Timestamp:

2021-01-07T23:09:56+01:00 (3 years ago)

Author:

aumont

Message:

numerous updates to PISCES, PISCES-QUOTA and the sediment module

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmort.F90 (modified) (8 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmort.F90

@@ -28 +28 @@
    REAL(wp), PUBLIC :: wchlp   !: Quadratic mortality rate of picophytoplankton
    REAL(wp), PUBLIC :: wchld   !: Quadratic mortality rate of diatoms
-   REAL(wp), PUBLIC :: wchldm  !: Maximum quadratic mortality rate of diatoms
    REAL(wp), PUBLIC :: mpratn  !: Linear mortality rate of nanophytoplankton
    REAL(wp), PUBLIC :: mpratp  !: Linear mortality rate of picophytoplankton
@@ -69 +68 @@
       !!---------------------------------------------------------------------
       INTEGER  :: ji, jj, jk
-      REAL(wp) :: zcompaph
+      REAL(wp) :: zcompaph, zlim1, zlim2
       REAL(wp) :: zfactfe, zfactch, zfactn, zfactp, zprcaca
       REAL(wp) :: ztortp , zrespp , zmortp
@@ -86 +85 @@
                ! blooms (Doney et al. 1996)
                ! -----------------------------------------------------
-               zrespp = wchln * 1.e6 * xstep * xdiss(ji,jj,jk) * zcompaph * trb(ji,jj,jk,jpphy)
+               zlim2   = xlimphy(ji,jj,jk) * xlimphy(ji,jj,jk)
+               zlim1   = 0.25 * ( 1. - zlim2 ) / ( 0.25 + zlim2 ) * trb(ji,jj,jk,jpphy)
+               zrespp = wchln * 1.e6 * xstep * zlim1 * xdiss(ji,jj,jk) * zcompaph
 
                ! Phytoplankton linear mortality
@@ -142 +143 @@
       !!---------------------------------------------------------------------
       INTEGER  :: ji, jj, jk
-      REAL(wp) :: zcompaph
+      REAL(wp) :: zcompaph, zlim1, zlim2
       REAL(wp) :: zfactfe, zfactch, zfactn, zfactp
       REAL(wp) :: ztortp , zrespp , zmortp 
@@ -158 +159 @@
                ! blooms (Doney et al. 1996)
                ! -----------------------------------------------------
-               zrespp = wchlp * 1.e6 * xstep * xdiss(ji,jj,jk) * zcompaph * trb(ji,jj,jk,jppic)
+               zlim2   = xlimpic(ji,jj,jk) * xlimpic(ji,jj,jk)
+               zlim1   = 0.25 * ( 1. - zlim2 ) / ( 0.25 + zlim2 ) * trb(ji,jj,jk,jppic)
+               zrespp = wchlp * 1.e6 * xstep * zlim1 * xdiss(ji,jj,jk) * zcompaph
 
                ! Phytoplankton linear mortality
@@ -229 +232 @@
                zlim2   = xlimdia(ji,jj,jk) * xlimdia(ji,jj,jk)
                zlim1   = 0.25 * ( 1. - zlim2 ) / ( 0.25 + zlim2 ) 
-               zrespp2 = 1.e6 * xstep * (  wchld + wchldm * zlim1 ) * xdiss(ji,jj,jk) * zcompadi * trb(ji,jj,jk,jpdia)
+               zrespp2 = 1.e6 * xstep * wchld * zlim1 * xdiss(ji,jj,jk) * zcompadi * trb(ji,jj,jk,jpdia)
 
                ! Phytoplankton linear mortality
@@ -291 +294 @@
       INTEGER :: ios   ! Local integer output status for namelist read
       !!
-      NAMELIST/namp5zmort/ wchln, wchlp, wchld, wchldm, mpratn, mpratp, mpratd
+      NAMELIST/namp5zmort/ wchln, wchlp, wchld, mpratn, mpratp, mpratd
       !!----------------------------------------------------------------------
 
@@ -310 +313 @@
          WRITE(numout,*) '    quadratic mortality of picophyto.         wchlp     =', wchlp
          WRITE(numout,*) '    quadratic mortality of diatoms            wchld     =', wchld
-         WRITE(numout,*) '    Additional quadratic mortality of diatoms wchldm    =', wchldm
          WRITE(numout,*) '    nanophyto. mortality rate                 mpratn    =', mpratn
          WRITE(numout,*) '    picophyto. mortality rate                 mpratp    =', mpratp