Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmort.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/p5zmort.F90 (modified) (8 diffs)

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

@@ -53 +53 @@
       !!---------------------------------------------------------------------
 
-      CALL p5z_nano            ! nanophytoplankton
-      CALL p5z_pico            ! picophytoplankton
-      CALL p5z_diat            ! diatoms
+      CALL p5z_mort_nano            ! nanophytoplankton
+      CALL p5z_mort_pico            ! picophytoplankton
+      CALL p5z_mort_diat            ! diatoms
 
    END SUBROUTINE p5z_mort
 
 
-   SUBROUTINE p5z_nano
-      !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE p5z_nano  ***
+   SUBROUTINE p5z_mort_nano
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p5z_mort_nano  ***
       !!
       !! ** Purpose :   Compute the mortality terms for nanophytoplankton
       !!
-      !! ** Method  : - ???
+      !! ** Method  : - Both quadratic and simili linear mortality terms
       !!---------------------------------------------------------------------
       INTEGER  :: ji, jj, jk
@@ -75 +75 @@
       !!---------------------------------------------------------------------
       !
-      IF( ln_timing )   CALL timing_start('p5z_nano')
+      IF( ln_timing )   CALL timing_start('p5z_mort_nano')
       !
       prodcal(:,:,:) = 0.  !: calcite production variable set to zero
@@ -82 +82 @@
             DO ji = 1, jpi
                zcompaph = MAX( ( trb(ji,jj,jk,jpphy) - 1e-9 ), 0.e0 )
+
                ! Quadratic mortality of nano due to aggregation during
                ! blooms (Doney et al. 1996)
@@ -127 +128 @@
        ENDIF
       !
-      IF( ln_timing )   CALL timing_stop('p5z_nano')
-      !
-   END SUBROUTINE p5z_nano
-
-
-   SUBROUTINE p5z_pico
-      !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE p5z_pico  ***
+      IF( ln_timing )   CALL timing_stop('p5z_mort_nano')
+      !
+   END SUBROUTINE p5z_mort_nano
+
+
+   SUBROUTINE p5z_mort_pico
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p5z_mort_pico  ***
       !!
       !! ** Purpose :   Compute the mortality terms for picophytoplankton
       !!
-      !! ** Method  : - ???
+      !! ** Method  : - Both quadratic and semilininear terms are used
       !!---------------------------------------------------------------------
       INTEGER  :: ji, jj, jk
@@ -147 +148 @@
       !!---------------------------------------------------------------------
       !
-      IF( ln_timing )   CALL timing_start('p5z_pico')
+      IF( ln_timing )   CALL timing_start('p5z_mort_pico')
       !
       DO jk = 1, jpkm1
@@ -191 +192 @@
        ENDIF
       !
-      IF( ln_timing )   CALL timing_stop('p5z_pico')
-      !
-   END SUBROUTINE p5z_pico
-
-
-   SUBROUTINE p5z_diat
-      !!---------------------------------------------------------------------
-      !!                     ***  ROUTINE p5z_diat  ***
+      IF( ln_timing )   CALL timing_stop('p5z_mort_pico')
+      !
+   END SUBROUTINE p5z_mort_pico
+
+
+   SUBROUTINE p5z_mort_diat
+      !!---------------------------------------------------------------------
+      !!                     ***  ROUTINE p5z_mort_diat  ***
       !!
       !! ** Purpose :   Compute the mortality terms for diatoms
@@ -211 +212 @@
       !!---------------------------------------------------------------------
       !
-      IF( ln_timing )   CALL timing_start('p5z_diat')
+      IF( ln_timing )   CALL timing_start('p5z_mort_diat')
       !
 
@@ -271 +272 @@
       ENDIF
       !
-      IF( ln_timing )   CALL timing_stop('p5z_diat')
-      !
-   END SUBROUTINE p5z_diat
+      IF( ln_timing )   CALL timing_stop('p5z_mort_diat')
+      !
+   END SUBROUTINE p5z_mort_diat