Changeset 7627 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z

Timestamp:

2017-02-01T11:22:41+01:00 (7 years ago)

Author:

aumont

Message:

change in the bacterial consumption of Fe + some more diags

Location:

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z

Files:

• p5zpoc.F90 (modified) (3 diffs)
• p5zrem.F90 (modified) (5 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zpoc.F90

@@ -214 +214 @@
                zopon2 = xremipn / xremipc * zremig * trb(ji,jj,jk,jpgon)
                zopop2 = xremipp / xremipc * zremig * trb(ji,jj,jk,jpgop)
-               zofer2 = xremipn / xremipc * zremig * trb(ji,jj,jk,jpbfe)
+               zofer2 = zremig * trb(ji,jj,jk,jpbfe)
 
                ! Update the appropriate tracers trends
@@ -396 +396 @@
                zopon  = xremipn / xremipc * zremip * trb(ji,jj,jk,jppon)
                zopop  = xremipp / xremipc * zremip * trb(ji,jj,jk,jppop)
-               zofer  = xremipn / xremipc * zremip * trb(ji,jj,jk,jpsfe)
+               zofer  = zremip * trb(ji,jj,jk,jpsfe)
 #if defined key_kriest
                zopoc2 = xremipc * zremip * trb(ji,jj,jk,jpnum)
@@ -422 +422 @@
           CALL iom_put( "REMINP" , zremipoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate
           CALL iom_put( "REMING" , zremigoc(:,:,:) * tmask(:,:,:) )  ! Remineralisation rate
-          CALL iom_put( "REMINF" , zfolimi(:,:,:) * tmask(:,:,:)  * 1.e+3 * rfact2r )  ! Remineralisation rate
+          CALL iom_put( "REMINF" , zfolimi(:,:,:) * tmask(:,:,:)  * 1.E9 * zrfact2 )  ! Remineralisation rate
       ENDIF
       !

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P5Z/p5zrem.F90

@@ -80 +80 @@
       CHARACTER (len=25) :: charout
       REAL(wp), POINTER, DIMENSION(:,:  ) :: ztempbac 
-      REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zwork1, zdepprod, zfacsi, zfacsib, znitr
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zdepprod, zdepeff, zfacsi, zfacsib
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: znitr, zfebact, zwork1
       !!---------------------------------------------------------------------
       !
@@ -87 +88 @@
       ! Allocate temporary workspace
       CALL wrk_alloc( jpi, jpj,      ztempbac                  )
-      CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zwork1, zfacsi, zfacsib, znitr )
+      CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib)
+      CALL wrk_alloc( jpi, jpj, jpk, znitr, zfebact, zwork1 )
 
       ! Initialisation of temprary arrys
       zdepprod(:,:,:) = 1._wp
+      zdepeff (:,:,:) = 0.3_wp
       ztempbac(:,:)   = 0._wp
+      zfebact (:,:,:) = 0._wp
       zfacsib(:,:,:)  = xsilab / ( 1.0 - xsilab )
       zfacsi(:,:,:)   = xsilab
@@ -111 +115 @@
                   zdepbac (ji,jj,jk) = zdepmin**0.683 * ztempbac(ji,jj)
                   zdepprod(ji,jj,jk) = zdepmin**0.273
+                  zdepeff (ji,jj,jk) = 0.3 * zdepmin**0.3
                ENDIF
             END DO
@@ -220 +225 @@
                zbactfer = feratb *  rfact2 * prmaxp(ji,jj,jk) * xlimbacl(ji,jj,jk)             &
                   &              * biron(ji,jj,jk) / ( xkferb + biron(ji,jj,jk) )    &
-                  &              * zdepprod(ji,jj,jk) * zdepbac(ji,jj,jk)
+                  &              * zdepprod(ji,jj,jk) * zdepeff(ji,jj,jk) * zdepbac(ji,jj,jk)
 #if defined key_kriest
                tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.05
                tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.05
+               zfebact(ji,jj,jk)   = zbactfer * 0.15
 #else
-               tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.12
-               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.09
-               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.03
+               tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.39
+               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.3
+               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.09
+               zfebact(ji,jj,jk)   = zbactfer * 0.39
 #endif
             END DO
@@ -284 +291 @@
           CALL iom_put( "DENIT" , denitrc(:,:,:) * rdenit * rno3 * tmask(:,:,:) * zrfact2  )  ! Denitrification
           CALL iom_put( "NIT"   , znitr(:,:,:) * rno3 * tmask(:,:,:) * zrfact2 )  !
-          CALL iom_put( "BACT"  , zdepbac(:,:,:) * 1.E6 * tmask(:,:,:) )  ! Bacterial biomass
+          CALL iom_put( "BACT", zdepbac(:,:,:) * 1.E6 * tmask(:,:,:) )  ! Bacterial biomass
+          CALL iom_put( "FEBACT" , zfebact(:,:,:) * 1E9 * tmask(:,:,:) * zrfact2  )  ! Denitrification
       ENDIF
       !
       CALL wrk_dealloc( jpi, jpj,      ztempbac                  )
-      CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zwork1, zfacsi, zfacsib, znitr )
+      CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zdepeff, zfacsi, zfacsib )
+      CALL wrk_dealloc( jpi, jpj, jpk, zwork1, zfebact, znitr )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p5z_rem')