Changeset 12496 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zmort.F90

Timestamp:

2020-03-02T12:33:55+01:00 (4 years ago)

Author:

aumont

Message:

uptake of ligand by phytoplankton + bug correction in the river supply (h_rnf) + cosmetic changes

File:

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

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

@@ -23 +23 @@
    PUBLIC   p4z_mort_init    
 
-   REAL(wp), PUBLIC ::   wchl     !:
+   REAL(wp), PUBLIC ::   wchln    !:
    REAL(wp), PUBLIC ::   wchld    !:
    REAL(wp), PUBLIC ::   wchldm   !:
-   REAL(wp), PUBLIC ::   mprat    !:
-   REAL(wp), PUBLIC ::   mprat2   !:
+   REAL(wp), PUBLIC ::   mpratn   !:
+   REAL(wp), PUBLIC ::   mpratd   !:
 
    !!----------------------------------------------------------------------
@@ -84 +84 @@
                !     Squared mortality of Phyto similar to a sedimentation term during
                !     blooms (Doney et al. 1996)
-               zrespp = wchl * 1.e6 * xstep * xdiss(ji,jj,jk) * zcompaph * zsizerat 
+               zrespp = wchln * 1.e6 * xstep * xdiss(ji,jj,jk) * zcompaph * zsizerat 
 
                !     Phytoplankton mortality. This mortality loss is slightly
                !     increased when nutrients are limiting phytoplankton growth
                !     as observed for instance in case of iron limitation.
-               ztortp = mprat * xstep * zcompaph / ( xkmort + trb(ji,jj,jk,jpphy) ) * zsizerat
+               ztortp = mpratn * xstep * zcompaph / ( xkmort + trb(ji,jj,jk,jpphy) ) * zsizerat
 
                zmortp = zrespp + ztortp
@@ -169 +169 @@
                !     Phytoplankton mortality. 
                !     ------------------------
-               ztortp2 = mprat2 * xstep * trb(ji,jj,jk,jpdia)  / ( xkmort + trb(ji,jj,jk,jpdia) ) * zcompadi 
+               ztortp2 = mpratd * xstep * trb(ji,jj,jk,jpdia)  / ( xkmort + trb(ji,jj,jk,jpdia) ) * zcompadi 
 
                zmortp2 = zrespp2 + ztortp2
@@ -218 +218 @@
       INTEGER ::   ios   ! Local integer
       !
-      NAMELIST/namp4zmort/ wchl, wchld, wchldm, mprat, mprat2
+      NAMELIST/namp4zmort/ wchln, wchld, wchldm, mpratn, mpratd
       !!----------------------------------------------------------------------
       !
@@ -237 +237 @@
       IF(lwp) THEN                         ! control print
          WRITE(numout,*) '   Namelist : namp4zmort'
-         WRITE(numout,*) '      quadratic mortality of phytoplankton        wchl   =', wchl
+         WRITE(numout,*) '      quadratic mortality of phytoplankton        wchln  =', wchln
          WRITE(numout,*) '      maximum quadratic mortality of diatoms      wchld  =', wchld
          WRITE(numout,*) '      maximum quadratic mortality of diatoms      wchldm =', wchldm
-         WRITE(numout,*) '      phytoplankton mortality rate                mprat  =', mprat
-         WRITE(numout,*) '      Diatoms mortality rate                      mprat2 =', mprat2
+         WRITE(numout,*) '      phytoplankton mortality rate                mpratn =', mpratn
+         WRITE(numout,*) '      Diatoms mortality rate                      mpratd =', mpratd
       ENDIF
       !