Changeset 2528 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zlim.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/p4zlim.F90 (modified) (9 diffs)

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

@@ -23 +23 @@
 
    PUBLIC p4z_lim    
+   PUBLIC p4z_lim_init    
 
    !! * Shared module variables
@@ -43 +44 @@
 #  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_lim( kt, jnt )
+   SUBROUTINE p4z_lim( kt )
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE p4z_lim  ***
@@ -59 +60 @@
       !! ** Method  : - ???
       !!---------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt, jnt ! ocean time step
+      INTEGER, INTENT(in)  :: kt
       INTEGER  ::   ji, jj, jk
       REAL(wp) ::   zlim1, zlim2, zlim3, zlim4, zno3, zferlim
@@ -67 +68 @@
 
 
-      IF( ( kt * jnt ) == nittrc000  )   CALL p4z_lim_init      ! Initialization (first time-step only)
-
-
-!  Tuning of the iron concentration to a minimum
-!  level that is set to the detection limit
-!  -------------------------------------
+      !  Tuning of the iron concentration to a minimum
+      !  level that is set to the detection limit
+      !  -------------------------------------
 
       DO jk = 1, jpkm1
@@ -85 +83 @@
       END DO
 
-!  Computation of a variable Ks for iron on diatoms
-!  taking into account that increasing biomass is
-!  made of generally bigger cells
-!  ------------------------------------------------
+      !  Computation of a variable Ks for iron on diatoms taking into account
+      !  that increasing biomass is made of generally bigger cells
+      !  ------------------------------------------------
 
       DO jk = 1, jpkm1
@@ -107 +104 @@
       END DO
 
-      DO jk = 1, jpkm1
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-    
-!      Michaelis-Menten Limitation term for nutrients
-!      Small flagellates
-!      -----------------------------------------------
+     !  Michaelis-Menten Limitation term for nutrients Small flagellates
+     !      -----------------------------------------------
+      DO jk = 1, jpkm1
+         DO jj = 1, jpj
+            DO ji = 1, jpi
               zdenom = 1. / &
                   & ( conc0 * concnnh4 + concnnh4 * trn(ji,jj,jk,jpno3) + conc0 * trn(ji,jj,jk,jpnh4) )
@@ -132 +127 @@
       END DO
 
-      DO jk = 1, jpkm1
-         DO jj = 1, jpj
-            DO ji = 1, jpi
-
-!   Michaelis-Menten Limitation term for nutrients Diatoms
-!   ----------------------------------------------
+      !   Michaelis-Menten Limitation term for nutrients Diatoms
+      !   ----------------------------------------------
+      DO jk = 1, jpkm1
+         DO jj = 1, jpj
+            DO ji = 1, jpi
               zdenom = 1. / &
                   & ( conc1 * concdnh4 + concdnh4 * trn(ji,jj,jk,jpno3) + conc1 * trn(ji,jj,jk,jpnh4) )
@@ -161 +155 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               ztemp = MAX( 0., tn(ji,jj,jk) )
+               ztemp = MAX( 0., tsn(ji,jj,jk,jp_tem) )
                xfracal(ji,jj,jk) = caco3r * xlimphy(ji,jj,jk)   &
                   &                       * MAX( 0.0001, ztemp / ( 2.+ ztemp ) )   &
@@ -181 +175 @@
       !!
       !! ** Method  :   Read the nampislim namelist and check the parameters
-      !!      called at the first timestep (nittrc000)
+      !!      called at the first timestep (nit000)
       !!
       !! ** input   :   Namelist nampislim