Changeset 4529 for trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

Timestamp:

2014-03-15T12:00:04+01:00 (10 years ago)

Author:

cetlod

Message:

Bugfix and minor improvment on PISCES, see ticket #1258

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 (modified) (7 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

@@ -35 +35 @@
    REAL(wp), PUBLIC ::  pislope         !:
    REAL(wp), PUBLIC ::  pislope2        !:
+   REAL(wp), PUBLIC ::  xadap           !:
    REAL(wp), PUBLIC ::  excret          !:
    REAL(wp), PUBLIC ::  excret2         !:
@@ -158 +159 @@
                   IF( etot(ji,jj,jk) > 1.E-3 ) THEN
                       ztn         = MAX( 0., tsn(ji,jj,jk,jp_tem) - 15. )
-                      zadap       = ztn / ( 2.+ ztn )
+                      zadap       = xadap * ztn / ( 2.+ ztn )
                       zconctemp   = MAX( 0.e0 , trn(ji,jj,jk,jpdia) - xsizedia )
                       zconctemp2  = trn(ji,jj,jk,jpdia) - zconctemp
@@ -227 +228 @@
       ENDIF
 
+
       !  Computation of a proxy of the N/C ratio
       !  ---------------------------------------
@@ -235 +237 @@
 !CDIR NOVERRCHK
             DO ji = 1, jpi
-                zval = ( xnanonh4(ji,jj,jk) + xnanono3(ji,jj,jk) ) * prmax(ji,jj,jk) / ( zprbio(ji,jj,jk) + rtrn )
-                quotan(ji,jj,jk) = MIN( 1., 0.5 + 0.5 * zval )
-                zval = ( xdiatnh4(ji,jj,jk) + xdiatno3(ji,jj,jk) ) * prmax(ji,jj,jk) / ( zprdia(ji,jj,jk) + rtrn )
-                quotad(ji,jj,jk) = MIN( 1., 0.5 + 0.5 * zval )
+                zval = MIN( xnanopo4(ji,jj,jk), ( xnanonh4(ji,jj,jk) + xnanono3(ji,jj,jk) ) )   &
+                &      * prmax(ji,jj,jk) / ( zprbio(ji,jj,jk) + rtrn )
+                quotan(ji,jj,jk) = MIN( 1., 0.2 + 0.8 * zval )
+                zval = MIN( xdiatpo4(ji,jj,jk), ( xdiatnh4(ji,jj,jk) + xdiatno3(ji,jj,jk) ) )   &
+                &      * prmax(ji,jj,jk) / ( zprdia(ji,jj,jk) + rtrn )
+                quotad(ji,jj,jk) = MIN( 1., 0.2 + 0.8 * zval )
             END DO
          END DO
@@ -274 +278 @@
             zmxltst = MAX( 0.e0, hmld(ji,jj) - heup(ji,jj) )
             zmxlday = zmxltst * zmxltst * r1_rday
-            zmixnano(ji,jj) = 1. - zmxlday / ( 3. + zmxlday )
+            zmixnano(ji,jj) = 1. - zmxlday / ( 2. + zmxlday )
             zmixdiat(ji,jj) = 1. - zmxlday / ( 4. + zmxlday )
          END DO
@@ -470 +474 @@
       !!----------------------------------------------------------------------
       !
-      NAMELIST/nampisprod/ pislope, pislope2, ln_newprod, bresp, excret, excret2,  &
+      NAMELIST/nampisprod/ pislope, pislope2, xadap, ln_newprod, bresp, excret, excret2,  &
          &                 chlcnm, chlcdm, chlcmin, fecnm, fecdm, grosip
       INTEGER :: ios                 ! Local integer output status for namelist read
@@ -491 +495 @@
          WRITE(numout,*) '    mean Si/C ratio                           grosip       =', grosip
          WRITE(numout,*) '    P-I slope                                 pislope      =', pislope
+         WRITE(numout,*) '    Acclimation factor to low light           xadap       =', xadap
          WRITE(numout,*) '    excretion ratio of nanophytoplankton      excret       =', excret
          WRITE(numout,*) '    excretion ratio of diatoms                excret2      =', excret2