Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/oce_sed.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/oce_sed.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/oce_sed.F90	(revision 13540)
@@ -13,5 +13,8 @@
    USE dom_oce , ONLY :   glamt     =>   glamt          !: longitude of t-point (degre)
    USE dom_oce , ONLY :   gphit     =>   gphit          !: latitude  of t-point (degre)
+!!st 
+#if ! defined key_qco
    USE dom_oce , ONLY :   e3t       =>   e3t            !: latitude  of t-point (degre)
+#endif
    USE dom_oce , ONLY :   e3t_1d    =>   e3t_1d         !: reference depth of t-points (m)
    USE dom_oce , ONLY :   gdepw_0   =>   gdepw_0        !: reference depth of t-points (m)
@@ -53,4 +56,2 @@
 
 END MODULE oce_sed
-
-
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedchem.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedchem.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedchem.F90	(revision 13540)
@@ -138,5 +138,5 @@
          CALL sed_chem_cst
       ELSE
-         DO_2D_11_11
+         DO_2D( 1, 1, 1, 1 )
             ikt = mbkt(ji,jj) 
             IF ( tmask(ji,jj,ikt) == 1 ) THEN
@@ -577,5 +577,5 @@
          saltprac(:) = salt(:) * 35.0 / 35.16504
       ELSE
-         saltprac(:) = temp(:)
+         saltprac(:) = salt(:)
       ENDIF
 
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/seddta.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/seddta.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/seddta.F90	(revision 13540)
@@ -24,4 +24,5 @@
    !! * Substitutions
 #  include "do_loop_substitute.h90"
+#  include "domzgr_substitute.h90"
    !! $Id$
 CONTAINS
@@ -95,5 +96,5 @@
       !    -----------------------------------------------------------
       IF (ln_sediment_offline) THEN
-         DO_2D_11_11
+         DO_2D( 1, 1, 1, 1 )
             ikt = mbkt(ji,jj)
             zwsbio4(ji,jj) = wsbio2 / rday
@@ -101,5 +102,5 @@
          END_2D
       ELSE
-         DO_2D_11_11
+         DO_2D( 1, 1, 1, 1 )
             ikt = mbkt(ji,jj)
             zdep = e3t(ji,jj,ikt,Kmm) / rDt_trc
@@ -110,5 +111,5 @@
 
       trc_data(:,:,:) = 0.
-      DO_2D_11_11
+      DO_2D( 1, 1, 1, 1 )
          ikt = mbkt(ji,jj)
          IF ( tmask(ji,jj,ikt) == 1 ) THEN
@@ -164,5 +165,5 @@
       CALL pack_arr ( jpoce,  rainrm_dta(1:jpoce,jscal), trc_data(1:jpi,1:jpj,14), iarroce(1:jpoce) )
       rainrm_dta(1:jpoce,jscal) = rainrm_dta(1:jpoce,jscal) * 1e-4
-      ! vector temperature [°C] and salinity 
+      ! vector temperature [�C] and salinity 
       CALL pack_arr ( jpoce,  temp(1:jpoce), trc_data(1:jpi,1:jpj,15), iarroce(1:jpoce) )
       CALL pack_arr ( jpoce,  salt(1:jpoce), trc_data(1:jpi,1:jpj,16), iarroce(1:jpoce) )
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedini.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedini.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedini.F90	(revision 13540)
@@ -135,5 +135,5 @@
       ! Determination of sediments number of points and allocate global variables
       epkbot(:,:) = 0.
-      DO_2D_11_11
+      DO_2D( 1, 1, 1, 1 )
          ikt = mbkt(ji,jj) 
          IF( tmask(ji,jj,ikt) == 1 ) epkbot(ji,jj) = e3t_1d(ikt)
@@ -247,5 +247,5 @@
       ! Computation of 1D array of sediments points
       indoce = 0
-      DO_2D_11_11
+      DO_2D( 1, 1, 1, 1 )
          IF (  epkbot(ji,jj) > 0. ) THEN
             indoce          = indoce + 1
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedinorg.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedinorg.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedinorg.F90	(revision 13540)
@@ -89,4 +89,5 @@
             zsolcpcl = zsolcpcl + solcp(ji,jk,jsclay) * dz(jk)
          END DO
+         zsolcpsi = MAX( zsolcpsi, rtrn )
          zsieq(ji) = sieqs(ji) * MAX(0.25, 1.0 - (0.045 * zsolcpcl / zsolcpsi )**0.58 )
          zsieq(ji) = MAX( rtrn, sieqs(ji) )
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedrst.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedrst.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedrst.F90	(revision 13540)
@@ -80,5 +80,5 @@
          IF(lwp) WRITE(numsed,*) &
              '             open sed restart.output NetCDF file: ',TRIM(clpath)//clname
-         CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed )
+         CALL iom_open( TRIM(clpath)//TRIM(clname), numrsw, ldwrt = .TRUE., kdlev = jpksed, cdcomp = 'SED' )
          lrst_sed = .TRUE.
       ENDIF
@@ -123,5 +123,5 @@
          cltra = TRIM(sedtrcd(jn))
          IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN
-            CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta(:,:,:,jn) )
+            CALL iom_get( numrsr, jpdom_auto, TRIM(cltra), zdta(:,:,:,jn) )
          ELSE
             zdta(:,:,:,jn) = 0.0
@@ -142,5 +142,5 @@
          cltra = TRIM(seddia3d(jn))
          IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN
-            CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta1(:,:,:,jn) )
+            CALL iom_get( numrsr, jpdom_auto, TRIM(cltra), zdta1(:,:,:,jn) )
          ELSE
             zdta1(:,:,:,jn) = 0.0
@@ -169,5 +169,5 @@
       cltra = "dbioturb"
       IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN
-         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:) )
+         CALL iom_get( numrsr, jpdom_auto, TRIM(cltra), zdta2(:,:,:) )
       ELSE
          zdta2(:,:,:) = 0.0
@@ -179,5 +179,5 @@
       cltra = "irrig"
       IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN
-         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:) )
+         CALL iom_get( numrsr, jpdom_auto, TRIM(cltra), zdta2(:,:,:) )
       ELSE
          zdta2(:,:,:) = 0.0
@@ -189,5 +189,5 @@
       cltra = "sedligand"
       IF( iom_varid( numrsr, TRIM(cltra) , ldstop = .FALSE. ) > 0 ) THEN
-         CALL iom_get( numrsr, jpdom_autoglo, TRIM(cltra), zdta2(:,:,:) )
+         CALL iom_get( numrsr, jpdom_auto, TRIM(cltra), zdta2(:,:,:) )
       ELSE
          zdta2(:,:,:) = 0.0
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedsfc.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedsfc.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/sedsfc.F90	(revision 13540)
@@ -48,5 +48,5 @@
 
 
-      DO_2D_11_11
+      DO_2D( 1, 1, 1, 1 )
          ikt = mbkt(ji,jj)
          IF ( tmask(ji,jj,ikt) == 1 ) THEN
Index: NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/trcdmp_sed.F90
===================================================================
--- NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/trcdmp_sed.F90	(revision 12511)
+++ NEMO/branches/2020/r12377_ticket2386/src/TOP/PISCES/SED/trcdmp_sed.F90	(revision 13540)
@@ -21,5 +21,5 @@
    USE trc             ! ocean passive tracers variables
    USE trcdta
-   USE prtctl_trc      ! Print control for debbuging
+   USE prtctl          ! Print control for debbuging
    USE iom
 
@@ -93,5 +93,5 @@
                CALL trc_dta( kt, Kmm, sf_trcdta(jl), rf_trfac(jl), ztrcdta )   ! read tracer data at nit000
                !
-               DO_2D_11_11
+               DO_2D( 1, 1, 1, 1 )
                   ikt = mbkt(ji,jj)
                   tr(ji,jj,ikt,jn,Kbb) = ztrcdta(ji,jj,ikt) + ( tr(ji,jj,ikt,jn,Kbb) -  ztrcdta(ji,jj,ikt) )     &
@@ -107,6 +107,6 @@
       IF( sn_cfctl%l_prttrc ) THEN
          WRITE(charout, FMT="('dmp ')")
-         CALL prt_ctl_trc_info(charout)
-         CALL prt_ctl_trc( tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm, clinfo2='trd' )
+         CALL prt_ctl_info( charout, cdcomp = 'top' )
+         CALL prt_ctl( tab4d_1=tr(:,:,:,:,Krhs), mask1=tmask, clinfo=ctrcnm, clinfo3='trd' )
       ENDIF
       !
