Changeset 7180 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

Timestamp:

2016-11-03T16:41:10+01:00 (8 years ago)

Author:

aumont

Message:

various bug fixes on iron chemistry

File:

• branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90 (modified) (9 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zprod.F90

@@ -81 +81 @@
       REAL(wp) ::   zprod, zproreg, zproreg2, zprochln, zprochld
       REAL(wp) ::   zmaxday, zdocprod, zpislopen, zpisloped
-      REAL(wp) ::   zrum, zcodel, zargu, zval, zfeup
+      REAL(wp) ::   zmxltst, zmxlday
+      REAL(wp) ::   zrum, zcodel, zargu, zval, zfeup, chlcnm_n, chlcdm_n
       REAL(wp) ::   zfact
       CHARACTER (len=25) :: charout
-      REAL(wp), POINTER, DIMENSION(:,:  ) :: zstrn, zw2d
+      REAL(wp), POINTER, DIMENSION(:,:  ) :: zstrn, zw2d, zmixnano, zmixdiat
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zpislopeadn, zpislopeadd, zprdia, zprbio, zprdch, zprnch, zysopt, zw3d   
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zprorcan, zprorcad, zprofed, zprofen, zpronewn, zpronewd
@@ -93 +94 @@
       !
       !  Allocate temporary workspace
-      CALL wrk_alloc( jpi, jpj,      zstrn )
+      CALL wrk_alloc( jpi, jpj,      zmixnano, zmixdiat, zstrn )
       CALL wrk_alloc( jpi, jpj, jpk, zpislopeadn, zpislopeadd, zprdia, zprbio, zprdch, zprnch, zysopt ) 
       CALL wrk_alloc( jpi, jpj, jpk, zmxl_fac, zmxl_chl )
@@ -104 +105 @@
 
       ! Computation of the optimal production
-      prmax(:,:,:) = 0.6_wp * r1_rday * tgfunc(:,:,:) 
+      prmax(:,:,:) = 0.8_wp * r1_rday * tgfunc(:,:,:) 
       IF( lk_degrad )  prmax(:,:,:) = prmax(:,:,:) * facvol(:,:,:) 
 
@@ -127 +128 @@
                IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN
                   zval = MAX( 1., zstrn(ji,jj) )
-                  zmxl_fac(ji,jj,jk) = zval
+                  IF( fsdept(ji,jj,jk) <= hmld(ji,jj) ) THEN
+                     zval = zval * MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn ))
+                  ENDIF
                   zmxl_chl(ji,jj,jk) = zval / 24.
-                  IF( fsdepw(ji,jj,jk+1) <= hmld(ji,jj) ) THEN
-                     zval = MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn ))
-                     zmxl_fac(ji,jj,jk) = zmxl_fac(ji,jj,jk) * zval
-                     zmxl_chl(ji,jj,jk) = zmxl_chl(ji,jj,jk) * zval
-                  ENDIF
-                  zmxl_fac(ji,jj,jk) = ( 1. - exp( -0.2 * zmxl_fac(ji,jj,jk) ) )
-                  zmxl_chl(ji,jj,jk) = zmxl_chl(ji,jj,jk)
+                  zmxl_fac(ji,jj,jk) = 1.5 * zval / ( 12. + zval )
                ENDIF
             END DO
@@ -143 +140 @@
       zprbio(:,:,:) = prmax(:,:,:) * zmxl_fac(:,:,:)
       zprdia(:,:,:) = zprbio(:,:,:)
+
+      ! Maximum light intensity
+      WHERE( zstrn(:,:) < 1.e0 ) zstrn(:,:) = 24.
 
       ! Computation of the P-I slope for nanos and diatoms
@@ -253 +253 @@
       END DO
 
+      !  Mixed-layer effect on production 
       !  Sea-ice effect on production
 
       DO jk = 1, jpkm1
-         zprbio(:,:,jk) = zprbio(:,:,jk) * ( 1. - fr_i(:,:) )
-         zprdia(:,:,jk) = zprdia(:,:,jk) * ( 1. - fr_i(:,:) )
-      END DO
-
+         DO jj = 1, jpj
+            DO ji = 1, jpi
+               zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * ( 1. - fr_i(ji,jj) )
+               zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * ( 1. - fr_i(ji,jj) )
+            END DO
+         END DO
+      END DO
 
       ! Computation of the various production terms 
@@ -303 +307 @@
                   zprod    = rday * zprorcan(ji,jj,jk) * zprnch(ji,jj,jk) * xlimphy(ji,jj,jk)
                   zprochln = chlcmin * 12. * zprorcan (ji,jj,jk)
-                  zprochln = zprochln + (chlcnm-chlcmin) * 12. * zprod / &
+                  chlcnm_n   = MIN ( chlcnm, ( chlcnm / (1. - 1.14 / 43.4 *tsn(ji,jj,jk,jp_tem))) * (1. - 1.14 / 43.4 * 20.))
+                  zprochln = zprochln + (chlcnm_n-chlcmin) * 12. * zprod / &
                                         & (  zpislopeadn(ji,jj,jk) * znanotot +rtrn)
                   !  production terms for diatoms ( chlorophyll )
@@ -309 +314 @@
                   zprod    = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * xlimdia(ji,jj,jk)
                   zprochld = chlcmin * 12. * zprorcad(ji,jj,jk)
-                  zprochld = zprochld + (chlcdm-chlcmin) * 12. * zprod / &
+                  chlcdm_n   = MIN ( chlcdm, ( chlcdm / (1. - 1.14 / 43.4 * tsn(ji,jj,jk,jp_tem))) * (1. - 1.14 / 43.4 * 20.))
+                  zprochld = zprochld + (chlcdm_n-chlcmin) * 12. * zprod / &
                                         & ( zpislopeadd(ji,jj,jk) * zdiattot +rtrn )
                   !   Update the arrays TRA which contain the Chla sources and sinks
@@ -485 +491 @@
      ENDIF
      !
-     CALL wrk_dealloc( jpi, jpj,      zstrn )
+     CALL wrk_dealloc( jpi, jpj,  zmixnano, zmixdiat,    zstrn )
      CALL wrk_dealloc( jpi, jpj, jpk, zpislopeadn, zpislopeadd, zprdia, zprbio, zprdch, zprnch, zysopt ) 
      CALL wrk_dealloc( jpi, jpj, jpk, zmxl_fac, zmxl_chl )