Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90

Timestamp:

2010-12-27T18:33:53+01:00 (13 years ago)

Author:

rblod

Message:

Update NEMOGCM from branch nemo_v3_3_beta

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90 (modified) (10 diffs, 1 prop)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zmort.F90

@@ -25 +25 @@
 
    PUBLIC   p4z_mort    
+   PUBLIC   p4z_mort_init    
 
 
@@ -35 +36 @@
      mpratm = 0.01_wp           !:
 
-   !! * Module variables
-   REAL(wp) :: zstep
-
-
 
    !!* Substitution
 #  include "top_substitute.h90"
    !!----------------------------------------------------------------------
-   !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) 
+   !! NEMO/TOP 3.3 , NEMO Consortium (2010)
    !! $Id$ 
-   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
+   !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
 CONTAINS
 
-   SUBROUTINE p4z_mort( kt, jnt )
+   SUBROUTINE p4z_mort( kt )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE p4z_mort  ***
@@ -59 +56 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt, jnt ! ocean time step
-      !!---------------------------------------------------------------------
-
-      IF( ( kt * jnt ) == nittrc000  )   CALL p4z_mort_init      ! Initialization (first time-step only)
-
-      zstep = rfact2 / rday      ! Time step duration for biology
+      INTEGER, INTENT(in) ::   kt ! ocean time step
+      !!---------------------------------------------------------------------
 
       CALL p4z_nano            ! nanophytoplankton
@@ -83 +76 @@
       INTEGER  :: ji, jj, jk
       REAL(wp) :: zcompaph
-      REAL(wp) :: zfactfe,zfactch,zprcaca,zfracal
-      REAL(wp) :: ztortp,zrespp,zmortp
+      REAL(wp) :: zfactfe, zfactch, zprcaca, zfracal
+      REAL(wp) :: ztortp , zrespp , zmortp , zstep
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
 
 
-#if defined key_trc_dia3d
+#if defined key_diatrc
      prodcal(:,:,:) = 0.  !: Initialisation of calcite production variable
 #endif
@@ -99 +92 @@
                zcompaph = MAX( ( trn(ji,jj,jk,jpphy) - 1e-8 ), 0.e0 )
 
-!     Squared mortality of Phyto similar to a sedimentation term during
-!     blooms (Doney et al. 1996)
-!     -----------------------------------------------------------------
-               zrespp = wchl * 1.e6 * zstep * xdiss(ji,jj,jk)   &
-# if defined key_off_degrad
-                  &        * facvol(ji,jj,jk)     &
+# if defined key_degrad
+               zstep =  xstep * facvol(ji,jj,jk)  
+# else
+               zstep =  xstep  
 # endif
-                  &        * zcompaph * trn(ji,jj,jk,jpphy)
-
-!     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 * zstep * trn(ji,jj,jk,jpphy)          &
-# if defined key_off_degrad
-                  &          * facvol(ji,jj,jk)     &
-# endif
-                  &   / ( xkmort + trn(ji,jj,jk,jpphy) ) * zcompaph
-
+               !     Squared mortality of Phyto similar to a sedimentation term during
+               !     blooms (Doney et al. 1996)
+               zrespp = wchl * 1.e6 * zstep * xdiss(ji,jj,jk) * zcompaph * trn(ji,jj,jk,jpphy) 
+
+               !     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 * trn(ji,jj,jk,jpphy) / ( xkmort + trn(ji,jj,jk,jpphy) ) * zcompaph
 
                zmortp = zrespp + ztortp
@@ -130 +117 @@
                tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zmortp * zfactfe
                zprcaca = xfracal(ji,jj,jk) * zmortp
-#if defined key_trc_dia3d
+#if defined key_diatrc
                prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca  ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo)
 #endif
@@ -169 +156 @@
       INTEGER  ::  ji, jj, jk
       REAL(wp) ::  zfactfe,zfactsi,zfactch, zcompadi
-      REAL(wp) ::  zrespp2, ztortp2, zmortp2
+      REAL(wp) ::  zrespp2, ztortp2, zmortp2, zstep
       CHARACTER (len=25) :: charout
  
@@ -175 +162 @@
 
 
-!    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
-!     ------------------------------------------------------------
+      !    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
+      !     ------------------------------------------------------------
 
       DO jk = 1, jpkm1
@@ -186 +173 @@
                zcompadi = MAX( ( trn(ji,jj,jk,jpdia) - 1e-8), 0. )
 
-!    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
-!     ------------------------------------------------------------
-
+               !    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
+               !     ------------------------------------------------------------
+
+# if defined key_degrad
+               zstep =  xstep * facvol(ji,jj,jk)  
+# else
+               zstep =  xstep  
+# endif
+               !  Phytoplankton respiration 
+               !     ------------------------
                zrespp2  = 1.e6 * zstep * (  wchl + wchld * ( 1.- xlimdia(ji,jj,jk) )  )    &
-# if defined key_off_degrad
-                  &       * facvol(ji,jj,jk)       &
-# endif
                   &       * xdiss(ji,jj,jk) * zcompadi * trn(ji,jj,jk,jpdia)
-                                                                               
-
-!     Phytoplankton mortality. 
-!     ------------------------
-               ztortp2  = mprat2 * zstep * trn(ji,jj,jk,jpdia)     &
-# if defined key_off_degrad
-                  &        * facvol(ji,jj,jk)       &
-# endif
-                  &      / ( xkmort + trn(ji,jj,jk,jpdia) ) * zcompadi
-
-                zmortp2 = zrespp2 + ztortp2
-
-!   Update the arrays tra which contains the biological sources and sinks
-!   ---------------------------------------------------------------------
+
+               !     Phytoplankton mortality. 
+               !     ------------------------
+               ztortp2  = mprat2 * zstep * trn(ji,jj,jk,jpdia)  / ( xkmort + trn(ji,jj,jk,jpdia) ) * zcompadi 
+
+               zmortp2 = zrespp2 + ztortp2
+
+               !   Update the arrays tra which contains the biological sources and sinks
+               !   ---------------------------------------------------------------------
                zfactch = trn(ji,jj,jk,jpdch) / ( trn(ji,jj,jk,jpdia) + rtrn )
                zfactfe = trn(ji,jj,jk,jpdfe) / ( trn(ji,jj,jk,jpdia) + rtrn )
@@ -249 +235 @@
       !!
       !! ** Method  :   Read the nampismort namelist and check the parameters
-      !!      called at the first timestep (nittrc000)
+      !!      called at the first timestep
       !!
       !! ** input   :   Namelist nampismort