Changeset 9450 for branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

Timestamp:

2018-03-30T18:08:29+02:00 (6 years ago)

Author:

aumont

Message:

debug PISCES code

File:

• branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90 (modified) (6 diffs)

branches/CNRS/dev_r8832_PISCO/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90

@@ -16 +16 @@
    USE sms_pisces      !  PISCES Source Minus Sink variables
    USE p4zche          !  chemical model
-   USE p4zprod         !  Growth rate of the 2 phyto groups
    USE p4zlim
    USE prtctl_trc      !  print control for debugging
@@ -70 +69 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) :: ztempbac
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zdepbac, zolimi, zdepprod, zfacsi, zw3d, zfacsib
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: zfebact, zdepeff
       !!---------------------------------------------------------------------
       !
@@ -76 +76 @@
       ! Allocate temporary workspace
       CALL wrk_alloc( jpi, jpj,      ztempbac                  )
-      CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib )
+      CALL wrk_alloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib, zfebact, zdepeff )
 
       ! Initialisation of temprary arrys
       zdepprod(:,:,:) = 1._wp
+      zdepeff (:,:,:) = 0.3_wp
       ztempbac(:,:)   = 0._wp
+      zfebact (:,:,:) = 0._wp
       zfacsib(:,:,:)  = xsilab / ( 1.0 - xsilab )
       zfacsi(:,:,:)   = xsilab
@@ -100 +102 @@
                   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 +223 @@
                zbactfer = feratb *  rfact2 * prmax(ji,jj,jk) * xlimbacl(ji,jj,jk)             &
                   &              * trb(ji,jj,jk,jpfer) / ( xkferb + trb(ji,jj,jk,jpfer) )    &
-                  &              * zdepprod(ji,jj,jk) * zdepbac(ji,jj,jk)
-               tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.16
-               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.12
-               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.04
+                  &              * zdepprod(ji,jj,jk) * zdepeff(ji,jj,jk) * zdepbac(ji,jj,jk)
+               tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zbactfer*0.33
+               tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + zbactfer*0.25
+               tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + zbactfer*0.08
+               zfebact(ji,jj,jk)   = zbactfer * 0.33
             END DO
          END DO
@@ -283 +287 @@
           ENDIF
           !
+          IF( iom_use( "BACT" ) )  THEN
+               zw3d(:,:,:) = zdepbac(:,:,:) * 1.E6 * tmask(:,:,:)  ! Bacterial biomass
+               CALL iom_put( "BACT", zw3d )
+          ENDIF
+          IF( iom_use( "FEBACT" ) )  THEN
+               zw3d(:,:,:) = zfebact(:,:,:) * 1E9 * tmask(:,:,:) * zrfact2   ! Bacterial iron consumption
+               CALL iom_put( "FEBACT" , zw3d )  
+          ENDIF
           CALL wrk_dealloc( jpi, jpj, jpk, zw3d )
        ENDIF
       !
       CALL wrk_dealloc( jpi, jpj,      ztempbac                  )
-      CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib )
+      CALL wrk_dealloc( jpi, jpj, jpk, zdepbac, zdepprod, zolimi, zfacsi, zfacsib, zfebact, zdepeff )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_rem')