Changeset 3446 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

Timestamp:

2012-08-10T13:13:55+02:00 (12 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB : 2nd step new PISCES updates from Olivier, see ticket #972

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 (modified) (8 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

@@ -77 +77 @@
       !
       INTEGER  ::   ji, jj, jk
-      REAL(wp) ::   zsilfac, zfact, znanotot, zdiattot, zconctemp, zconctemp2
+      REAL(wp) ::   zsilfac, znanotot, zdiattot, zconctemp, zconctemp2
       REAL(wp) ::   zratio, zmax, zsilim, ztn, zadap
       REAL(wp) ::   zlim, zsilfac2, zsiborn, zprod, zproreg, zproreg2
@@ -135 +135 @@
                DO ji = 1, jpi
                   IF( etot(ji,jj,jk) > 1.E-3 ) THEN
-                  zval = MAX( 1., zstrn(ji,jj) )
-                  zval = 1.5 * zval / ( 12. + zval )
-                  zprbio(ji,jj,jk) = prmax(ji,jj,jk) * zval
-                  zprdia(ji,jj,jk) = zprbio(ji,jj,jk)
+                     zval = MAX( 1., zstrn(ji,jj) )
+                     zval = 1.5 * zval / ( 12. + zval )
+                     zprbio(ji,jj,jk) = prmax(ji,jj,jk) * zval
+                     zprdia(ji,jj,jk) = zprbio(ji,jj,jk)
                   ENDIF
                END DO
@@ -156 +156 @@
 !CDIR NOVERRCHK
                DO ji = 1, jpi
-
                   ! Computation of the P-I slope for nanos and diatoms
                   IF( etot(ji,jj,jk) > 1.E-3 ) THEN
-                      ztn    = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. )
-                      zadap  = ztn / ( 2.+ ztn )
-
+                      ztn         = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. )
+                      zadap       = ztn / ( 2.+ ztn )
                       zconctemp   = MAX( 0.e0 , trn(ji,jj,jk,jpdia) - xsizedia )
                       zconctemp2  = trn(ji,jj,jk,jpdia) - zconctemp
-
-                      znanotot = enano(ji,jj,jk) * zstrn(ji,jj)
-                      zdiattot = ediat(ji,jj,jk) * zstrn(ji,jj)
-
-                      zfact  = EXP( -0.21 * znanotot )
-                      zpislopead (ji,jj,jk) = pislope  * ( 1.+ zadap  * zfact )  &
+                      znanotot    = enano(ji,jj,jk) * zstrn(ji,jj)
+                      zdiattot    = ediat(ji,jj,jk) * zstrn(ji,jj)
+                      !
+                      zpislopead (ji,jj,jk) = pislope * ( 1.+ zadap  * EXP( -0.21 * znanotot ) )  &
                          &                   * trn(ji,jj,jk,jpnch) /( trn(ji,jj,jk,jpphy) * 12. + rtrn)
-
+                      !
                       zpislopead2(ji,jj,jk) = (pislope * zconctemp2 + pislope2 * zconctemp) / ( trn(ji,jj,jk,jpdia) + rtrn )   &
                          &                   * trn(ji,jj,jk,jpdch) /( trn(ji,jj,jk,jpdia) * 12. + rtrn)
@@ -201 +197 @@
                   ! Computation of the P-I slope for nanos and diatoms
                   IF( etot(ji,jj,jk) > 1.E-3 ) THEN
-                      ztn    = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. )
-                      zadap  = ztn / ( 2.+ ztn )
-
+                      ztn         = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. )
+                      zadap       = ztn / ( 2.+ ztn )
                       zconctemp   = MAX( 0.e0 , trn(ji,jj,jk,jpdia) - xsizedia )
                       zconctemp2  = trn(ji,jj,jk,jpdia) - zconctemp
-
-                      zfact  = EXP( -0.21 * enano(ji,jj,jk) )
-                      zpislopead (ji,jj,jk) = pislope  * ( 1.+ zadap  * zfact )
+                      !
+                      zpislopead (ji,jj,jk) = pislope  * ( 1.+ zadap  * EXP( -0.21 * enano(ji,jj,jk) ) )
                       zpislopead2(ji,jj,jk) = (pislope * zconctemp2 + pislope2 * zconctemp)
 
@@ -264 +258 @@
                   zsilim = MIN( zprdia(ji,jj,jk) / ( prmax(ji,jj,jk) + rtrn ), xlimsi(ji,jj,jk) )
                   zsilfac = 4.4 * EXP( -4.23 * zsilim ) * MAX( 0.e0, MIN( 1., 2.2 * ( zlim - 0.5 ) )  ) + 1.e0
-                  zsiborn = MAX( 0.e0, ( trn(ji,jj,jk,jpsil) - 15.e-6 ) )
-                  zsilfac2 = 1.+ 2.* zsiborn / ( zsiborn + xksi2 )
-                  zsilfac = MIN( 5.4, zsilfac * zsilfac2)
-                  zysopt(ji,jj,jk) = grosip * zlim * zsilfac
+                  zsiborn = trn(ji,jj,jk,jpsil) * trn(ji,jj,jk,jpsil) * trn(ji,jj,jk,jpsil)
+                  IF (gphit(ji,jj) < -30 ) THEN
+                    zsilfac2 = 1. + 2. * zsiborn / ( zsiborn + xksi2**3 )
+                  ELSE
+                    zsilfac2 = 1. +      zsiborn / ( zsiborn + xksi2**3 )
+                  ENDIF
+                  zysopt(ji,jj,jk) = grosip * zlim * zsilfac * zsilfac2
               ENDIF
             END DO
@@ -312 +309 @@
                   zprofen(ji,jj,jk) = fecnm * prmax(ji,jj,jk)  &
                   &             * ( 4. - 4.5 * xlimnfe(ji,jj,jk) / ( xlimnfe(ji,jj,jk) + 0.5 ) )    &
-                  &             * trn(ji,jj,jk,jpfer) / ( trn(ji,jj,jk,jpfer) + concnfe(ji,jj,jk) )  &
+                  &             * biron(ji,jj,jk) / ( biron(ji,jj,jk) + concnfe(ji,jj,jk) )  &
                   &             * zmax * trn(ji,jj,jk,jpphy) * rfact2
                   !  production terms for diatomees
@@ -323 +320 @@
                   zprofed(ji,jj,jk) = fecdm * prmax(ji,jj,jk)  &
                   &             * ( 4. - 4.5 * xlimdfe(ji,jj,jk) / ( xlimdfe(ji,jj,jk) + 0.5 ) )    &
-                  &             * trn(ji,jj,jk,jpfer) / ( trn(ji,jj,jk,jpfer) + concdfe(ji,jj,jk) )  &
+                  &             * biron(ji,jj,jk) / ( biron(ji,jj,jk) + concdfe(ji,jj,jk) )  &
                   &             * zmax * trn(ji,jj,jk,jpdia) * rfact2
                ENDIF
@@ -426 +423 @@
          IF( lk_iomput ) THEN
            IF( jnt == nrdttrc ) THEN
-              CALL iom_put( "PPPHY" , zprorca (:,:,:) * zrfact2 * tmask(:,:,:) )  ! primary production by nanophyto
-              CALL iom_put( "PPPHY2", zprorcad(:,:,:) * zrfact2 * tmask(:,:,:) )  ! primary production by diatom
-              CALL iom_put( "PPNEWN", zpronew (:,:,:) * zrfact2 * tmask(:,:,:) )  ! new primary production by nanophyto
-              CALL iom_put( "PPNEWD", zpronewd(:,:,:) * zrfact2 * tmask(:,:,:) )  ! new primary production by diatom
-              CALL iom_put( "PBSi"  , zprorcad(:,:,:) * zrfact2 * tmask(:,:,:) * zysopt(:,:,:) ) ! biogenic silica production
-              CALL iom_put( "PFeD"  , zprofed (:,:,:) * zrfact2 * tmask(:,:,:) )  ! biogenic iron production by diatom
-              CALL iom_put( "PFeN"  , zprofen (:,:,:) * zrfact2 * tmask(:,:,:) )  ! biogenic iron production by nanophyto
+              CALL iom_put( "PPPHY"  , zprorca (:,:,:) * zrfact2 * tmask(:,:,:) )  ! primary production by nanophyto
+              CALL iom_put( "PPPHY2" , zprorcad(:,:,:) * zrfact2 * tmask(:,:,:) )  ! primary production by diatom
+              CALL iom_put( "PPNEWN" , zpronew (:,:,:) * zrfact2 * tmask(:,:,:) )  ! new primary production by nanophyto
+              CALL iom_put( "PPNEWD" , zpronewd(:,:,:) * zrfact2 * tmask(:,:,:) )  ! new primary production by diatom
+              CALL iom_put( "PBSi"   , zprorcad(:,:,:) * zrfact2 * tmask(:,:,:) * zysopt(:,:,:) ) ! biogenic silica production
+              CALL iom_put( "PFeD"   , zprofed (:,:,:) * zrfact2 * tmask(:,:,:) )  ! biogenic iron production by diatom
+              CALL iom_put( "PFeN"   , zprofen (:,:,:) * zrfact2 * tmask(:,:,:) )  ! biogenic iron production by nanophyto
               CALL iom_put( "Mumax"  , prmax(:,:,:) * tmask(:,:,:) )  ! Maximum growth rate
-              CALL iom_put( "MuN"  , zprbio(:,:,:) * xlimphy(:,:,:) * tmask(:,:,:) )  ! Realized growth rate for nanophyto
-              CALL iom_put( "MuD"  , zprdia(:,:,:) * xlimdia(:,:,:) * tmask(:,:,:) )  ! Realized growth rate for diatoms
-              CALL iom_put( "LNnut"  , xlimphy (:,:,:) * tmask(:,:,:) )  ! Nutrient limitation term
-              CALL iom_put( "LDnut"  , xlimdia (:,:,:) * tmask(:,:,:) )  ! Nutrient limitation term
-              CALL iom_put( "LNFe"  , xlimnfe (:,:,:) * tmask(:,:,:) )  ! Iron limitation term
-              CALL iom_put( "LDFe"  , xlimdfe (:,:,:) * tmask(:,:,:) )  ! Iron limitation term
-              CALL iom_put( "LNlight"  , zprbio (:,:,:) / (prmax(:,:,:) + rtrn) * tmask(:,:,:) )  ! light limitation term
-              CALL iom_put( "LDlight"  , zprdia (:,:,:) / (prmax(:,:,:) + rtrn) * tmask(:,:,:) )  ! light limitation term
+              CALL iom_put( "MuN"    , zprbio(:,:,:) * xlimphy(:,:,:) * tmask(:,:,:) )  ! Realized growth rate for nanophyto
+              CALL iom_put( "MuD"    , zprdia(:,:,:) * xlimdia(:,:,:) * tmask(:,:,:) )  ! Realized growth rate for diatoms
+              CALL iom_put( "LNlight", zprbio (:,:,:) / (prmax(:,:,:) + rtrn) * tmask(:,:,:) )  ! light limitation term
+              CALL iom_put( "LDlight", zprdia (:,:,:) / (prmax(:,:,:) + rtrn) * tmask(:,:,:) )  ! light limitation term
            ENDIF
          ELSE