Changeset 15075 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/SED/sedstp.F90

Timestamp:

2021-07-02T17:15:57+02:00 (3 years ago)

Author:

aumont

Message:

major update of the sediment module

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/SED/sedstp.F90 (modified) (5 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/SED/sedstp.F90

@@ -9 +9 @@
    USE sedco3   ! carbonate in sediment pore water
    USE sedsol   ! Organic reactions and diffusion
-   USE sedbtb   ! bioturbation
    USE sedadv   ! vertical advection
    USE sedsfc   ! sediment surface data
    USE sedrst   ! restart
    USE sedwri   ! outputs
+   USE sedini
    USE trcdmp_sed
    USE lib_mpp         ! distribued memory computing library
@@ -23 +23 @@
    !! * Routine accessibility
    PUBLIC sed_stp  ! called by step.F90
+
+   !! * Substitutions
+#  include "do_loop_substitute.h90"
+#  include "domzgr_substitute.h90"
 
    !! $Id$
@@ -42 +46 @@
       !!        !  06-04 (C. Ethe)  Re-organization
       !!----------------------------------------------------------------------
-      INTEGER, INTENT(in) ::   kt                ! number of iteration
-      INTEGER, INTENT(in) ::   Kbb, Kmm, Krhs    ! time level indices
+      INTEGER, INTENT(in) ::   kt       ! number of iteration
+      INTEGER, INTENT(in) ::   Kbb, Kmm, Krhs  ! time level indices
+
+      INTEGER :: ji,jk,js,jn,jw,jkmax,jsmax
       !!----------------------------------------------------------------------
       IF( ln_timing )           CALL timing_start('sed_stp')
@@ -53 +59 @@
 
       dtsed  = rDt_trc
-!      dtsed2 = dtsed
       IF (kt /= nitsed000) THEN
-         CALL sed_dta( kt, Kbb, Kmm )       ! Load  Data for bot. wat. Chem and fluxes
+         CALL sed_dta( kt, Kbb, Kmm )    ! Load  Data for bot. wat. Chem and fluxes
       ENDIF
 
@@ -63 +68 @@
          ENDIF
 
-         CALL sed_btb( kt )        ! 1st pass of bioturbation at t+1/2
          CALL sed_sol( kt )        ! Solute diffusion and reactions 
-         CALL sed_btb( kt )        ! 2nd pass of bioturbation at t+1
          CALL sed_adv( kt )         ! advection
          CALL sed_co3( kt )         ! pH actualization for saving
-         IF (ln_sed_2way) CALL sed_sfc( kt, Kbb )         ! Give back new bottom wat chem to tracer model
+
+         IF (ln_sed_2way) CALL sed_sfc( kt, Kbb )   ! Give back new bottom wat chem to tracer model
       ENDIF
       CALL sed_wri( kt )         ! outputs
       IF( kt == nitsed000 ) THEN
           CALL iom_close( numrsr )       ! close input tracer restart file
-          IF(lrxios) CALL iom_context_finalize(      cr_sedrst_cxt  )
-!         IF(lwm) CALL FLUSH( numont )   ! flush namelist output
+!          IF(lwm) CALL FLUSH( numont )   ! flush namelist output
       ENDIF
       IF( lrst_sed )            CALL sed_rst_wri( kt )   ! restart file output
 
-      IF( kt == nitsedend )  CLOSE( numsed )
+      IF( kt == nitsedend )     CLOSE( numsed )
 
-      IF( ln_timing )   CALL timing_stop('sed_stp')
+      IF( ln_timing )           CALL timing_stop('sed_stp')
 
    END SUBROUTINE sed_stp