Changeset 14276 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90

Timestamp:

2021-01-07T23:09:56+01:00 (3 years ago)

Author:

aumont

Message:

numerous updates to PISCES, PISCES-QUOTA and the sediment module

File:

• NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90 (modified) (13 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zsms.F90

@@ -64 +64 @@
       REAL(wp) ::  ztra
       CHARACTER (len=25) :: charout
+      REAL(wp), DIMENSION(jpi,jpj,jpk,jp_pisces) :: ztrbbio
       !!---------------------------------------------------------------------
       !
@@ -102 +103 @@
          END DO
       ENDIF
+
+      DO jn = jp_pcs0, jp_pcs1              !   Store the tracer concentrations before entering PISCES
+         ztrbbio(:,:,:,jn) = trb(:,:,:,jn)
+      END DO
       !
       IF( ll_sbc ) CALL p4z_sbc( kt )   ! external sources of nutrients 
@@ -147 +152 @@
             tra(:,:,:,jn) = 0._wp
          END DO
-         ! Euler-forward temporal scheme
-         IF( ln_top_euler ) THEN
-            DO jn = jp_pcs0, jp_pcs1
-               trn(:,:,:,jn) = trb(:,:,:,jn)
-            END DO
-         ENDIF
+         !
       END DO
-
+      !
+#endif
+      !
+      ! If ln_sediment is set to .true. then the sediment module is called
+      IF( ln_sediment ) THEN 
+         !
+         CALL sed_model( kt )     !  Main program of Sediment model
+         !
+      ENDIF
+      !
+      !
+      DO jn = jp_pcs0, jp_pcs1
+         tra(:,:,:,jn) = ( trb(:,:,:,jn) - ztrbbio(:,:,:,jn) ) * rfactr
+         trb(:,:,:,jn) = ztrbbio(:,:,:,jn)
+         ztrbbio(:,:,:,jn) = 0._wp
+      END DO
       !
       IF( l_trdtrc ) THEN
@@ -161 +176 @@
          END DO
       END IF
-#endif
-      !
-      ! If ln_sediment is set to .true. then the sediment module is called
-      IF( ln_sediment ) THEN 
-         !
-         CALL sed_model( kt )     !  Main program of Sediment model
-         ! Eulor forward temporal scheme
-         IF( ln_top_euler ) THEN
-            DO jn = jp_pcs0, jp_pcs1
-               trn(:,:,:,jn) = trb(:,:,:,jn)
-            END DO
-         ENDIF
-         !
-      ENDIF
       !
       IF( lrst_trc )  CALL p4z_rst( kt, 'WRITE' )  !* Write PISCES informations in restart file 
@@ -199 +200 @@
       INTEGER :: ios                 ! Local integer output status for namelist read
       !!
-      NAMELIST/nampisbio/ nrdttrc, wsbio, xkmort, ferat3, wsbio2, wsbio2max, wsbio2scale,    &
-         &                   ldocp, ldocz, lthet, no3rat3, po4rat3
+      NAMELIST/nampisbio/ nrdttrc, wsbio, xkmort, feratz, feratm, wsbio2, wsbio2max,    &
+         &                wsbio2scale, ldocp, ldocz, lthet, no3rat3, po4rat3
          !
       NAMELIST/nampisdmp/ ln_pisdmp, nn_pisdmp
@@ -226 +227 @@
          WRITE(numout,*) '      half saturation constant for mortality    xkmort      =', xkmort 
          IF( ln_p5z ) THEN
-            WRITE(numout,*) '      N/C in zooplankton                        no3rat3     =', no3rat3
-            WRITE(numout,*) '      P/C in zooplankton                        po4rat3     =', po4rat3
-         ENDIF
-         WRITE(numout,*) '      Fe/C in zooplankton                       ferat3      =', ferat3
+            WRITE(numout,*) '      N/C in zooplankton                     no3rat3     =', no3rat3
+            WRITE(numout,*) '      P/C in zooplankton                     po4rat3     =', po4rat3
+         ENDIF
+         WRITE(numout,*) '      Fe/C in microzooplankton                  feratz      =', feratz
+         WRITE(numout,*) '      Fe/C in microzooplankton                  feratz      =', feratm
          WRITE(numout,*) '      Big particles sinking speed               wsbio2      =', wsbio2
          WRITE(numout,*) '      Big particles maximum sinking speed       wsbio2max   =', wsbio2max
@@ -313 +315 @@
          ENDIF
 
+         ! Read the Fe3 consumption term by phytoplankton
+         IF( iom_varid( numrtr, 'Consfe3', ldstop = .FALSE. ) > 0 ) THEN
+            CALL iom_get( numrtr, jpdom_autoglo, 'Consfe3' , consfe3(:,:,:)  )
+         ELSE
+            consfe3(:,:,:) = 0._wp
+         ENDIF
+
+
          ! Read the cumulative total flux. If not in the restart file, it is set to 0          
          IF( iom_varid( numrtr, 'tcflxcum', ldstop = .FALSE. ) > 0 ) THEN  ! cumulative total flux of carbon
@@ -326 +336 @@
             sized(:,:,:) = 1.
          ENDIF
+         sized(:,:,:) = MAX( 1.0, sized(:,:,:) )
          IF( iom_varid( numrtr, 'sizen', ldstop = .FALSE. ) > 0 ) THEN
             CALL iom_get( numrtr, jpdom_autoglo, 'sizen' , sizen(:,:,:)  )
@@ -331 +342 @@
             sizen(:,:,:) = 1.
          ENDIF
+         sizen(:,:,:) = MAX( 1.0, sizen(:,:,:) )
 
          ! PISCES-QUOTA specific part
@@ -338 +350 @@
             IF( iom_varid( numrtr, 'sizep', ldstop = .FALSE. ) > 0 ) THEN
                CALL iom_get( numrtr, jpdom_autoglo, 'sizep' , sizep(:,:,:)  )
+               sizep(:,:,:) = MAX( 1.0, sizep(:,:,:) )
             ELSE
                sizep(:,:,:) = 1.
@@ -353 +366 @@
          CALL iom_rstput( kt, nitrst, numrtw, 'Silicalim', xksi(:,:) )    ! Si 1/2 saturation constant
          CALL iom_rstput( kt, nitrst, numrtw, 'Silicamax', xksimax(:,:) ) ! Si max concentration
+         CALL iom_rstput( kt, nitrst, numrtw, 'Consfe3', consfe3(:,:,:) ) ! Si max concentration
          CALL iom_rstput( kt, nitrst, numrtw, 'tcflxcum', t_oce_co2_flx_cum ) ! Cumulative CO2 flux
          CALL iom_rstput( kt, nitrst, numrtw, 'sizen', sizen(:,:,:) )  ! Size of nanophytoplankton
@@ -377 +391 @@
       !
       REAL(wp) ::  alkmean = 2426.     ! mean value of alkalinity ( Glodap ; for Goyet 2391. )
-      REAL(wp) ::  po4mean = 2.165     ! mean value of phosphate
-      REAL(wp) ::  no3mean = 30.90     ! mean value of nitrate
-      REAL(wp) ::  silmean = 91.51     ! mean value of silicate
+      REAL(wp) ::  po4mean = 2.174     ! mean value of phosphate
+      REAL(wp) ::  no3mean = 31.00     ! mean value of nitrate
+      REAL(wp) ::  silmean = 90.33     ! mean value of silicate
       !
       REAL(wp) :: zarea, zalksumn, zpo4sumn, zno3sumn, zsilsumn
@@ -534 +548 @@
          zwork(:,:,:) =   trn(:,:,:,jpfer) + trn(:,:,:,jpnfe) + trn(:,:,:,jpdfe)   &
             &         +   trn(:,:,:,jpbfe) + trn(:,:,:,jpsfe)                      &
-            &         + ( trn(:,:,:,jpzoo) + trn(:,:,:,jpmes) )  * ferat3    
+            &         + trn(:,:,:,jpzoo) * feratz + trn(:,:,:,jpmes) * feratm
          !
          ferbudget = glob_sum( 'p4zsms', zwork(:,:,:) * cvol(:,:,:)  )