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

Timestamp:

2020-07-02T20:34:16+02:00 (4 years ago)

Author:

aumont

Message:

update of the PISCES comments

File:

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

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

@@ -48 +48 @@
       !!---------------------------------------------------------------------
       !
-      CALL p4z_nano            ! nanophytoplankton
-      CALL p4z_diat            ! diatoms
+      CALL p4z_mort_nano            ! nanophytoplankton
+      CALL p4z_mort_diat            ! diatoms
       !
    END SUBROUTINE p4z_mort
 
 
-   SUBROUTINE p4z_nano
-      !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE p4z_nano  ***
+   SUBROUTINE p4z_mort_nano
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_mort_nano  ***
       !!
       !! ** Purpose :   Compute the mortality terms for nanophytoplankton
       !!
-      !! ** Method  : - ???
+      !! ** Method  :   Both quadratic and simili linear mortality terms
       !!---------------------------------------------------------------------
       INTEGER  ::   ji, jj, jk
@@ -69 +69 @@
       !!---------------------------------------------------------------------
       !
-      IF( ln_timing )   CALL timing_start('p4z_nano')
+      IF( ln_timing )   CALL timing_start('p4z_mort_nano')
       !
       prodcal(:,:,:) = 0._wp   ! calcite production variable set to zero
@@ -101 +101 @@
                tra(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) - zmortp * zfactch
                tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zmortp * zfactfe
+
                ! Production PIC particles due to mortality
                zprcaca = xfracal(ji,jj,jk) * zmortp
                prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca  ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo)
+
+               ! POC associated with the shell is supposed to be routed to 
+               ! big particles because of the ballasting effect
                zfracal = 0.5 * xfracal(ji,jj,jk)
                tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprcaca
@@ -112 +116 @@
                prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + ( 1. - zfracal ) * zmortp
                prodgoc(ji,jj,jk) = prodgoc(ji,jj,jk) + zfracal * zmortp
+
+               ! Update the arrays TRA which contains the biological sources and sinks
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ( 1. - zfracal ) * zmortp * zfactfe
                tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zfracal * zmortp * zfactfe
@@ -124 +130 @@
        ENDIF
       !
-      IF( ln_timing )   CALL timing_stop('p4z_nano')
-      !
-   END SUBROUTINE p4z_nano
-
-
-   SUBROUTINE p4z_diat
-      !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE p4z_diat  ***
+      IF( ln_timing )   CALL timing_stop('p4z_mort_nano')
+      !
+   END SUBROUTINE p4z_mort_nano
+
+
+   SUBROUTINE p4z_mort_diat
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p4z_mort_diat  ***
       !!
       !! ** Purpose :   Compute the mortality terms for diatoms
       !!
-      !! ** Method  : - ???
+      !! ** Method  : - Both quadratic and simili linear mortality terms
       !!---------------------------------------------------------------------
       INTEGER  ::   ji, jj, jk
@@ -144 +150 @@
       !!---------------------------------------------------------------------
       !
-      IF( ln_timing )   CALL timing_start('p4z_diat')
+      IF( ln_timing )   CALL timing_start('p4z_mort_diat')
       !
       ! Aggregation term for diatoms is increased in case of nutrient
       ! stress as observed in reality. The stressed cells become more
       ! sticky and coagulate to sink quickly out of the euphotic zone
+      ! This is due to the production of EPS by stressed cells
       ! -------------------------------------------------------------
 
@@ -183 +190 @@
                tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zmortp2 * zfactsi
                tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zmortp2 * zfactsi
+
+               ! Half of the linear mortality term is routed to big particles
+               ! becaue of the ballasting effect
                tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zrespp2 + 0.5 * ztortp2
                tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) + 0.5 * ztortp2
@@ -199 +209 @@
       ENDIF
       !
-      IF( ln_timing )   CALL timing_stop('p4z_diat')
-      !
-   END SUBROUTINE p4z_diat
+      IF( ln_timing )   CALL timing_stop('p4z_mort_diat')
+      !
+   END SUBROUTINE p4z_mort_diat