diff --git a/src/TOP/PISCES/P4Z/p4zsed.F90 b/src/TOP/PISCES/P4Z/p4zsed.F90
index 80e1e3fb042fcc822fa0df2bf2b40132d069f562..629df28d037f7cc3db49f1a13bc9fcdeb1943e94 100644
--- a/src/TOP/PISCES/P4Z/p4zsed.F90
+++ b/src/TOP/PISCES/P4Z/p4zsed.F90
@@ -90,16 +90,16 @@ CONTAINS
zsedcal (:,:) = 0.e0
zsedc (:,:) = 0.e0
- IF( .NOT.lk_sed ) THEN
- ! OA: Warning, the following part is necessary to avoid CFL problems above the sediments
- ! --------------------------------------------------------------------
- DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
- ikt = mbkt(ji,jj)
- zdep = e3t(ji,jj,ikt,Kmm) / xstep
- zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) )
- zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) )
- END_2D
+ ! OA: Warning, the following part is necessary to avoid CFL problems above the sediments
+ ! --------------------------------------------------------------------
+ DO_2D( nn_hls, nn_hls, nn_hls, nn_hls )
+ ikt = mbkt(ji,jj)
+ zdep = e3t(ji,jj,ikt,Kmm) / xstep
+ zwsbio4(ji,jj) = MIN( 0.99 * zdep, wsbio4(ji,jj,ikt) )
+ zwsbio3(ji,jj) = MIN( 0.99 * zdep, wsbio3(ji,jj,ikt) )
+ END_2D
+ IF( .NOT.lk_sed ) THEN
! Computation of the sediment denitrification proportion: The metamodel from midlleburg (2006) is being used
! Computation of the fraction of organic matter that is permanently buried from Dunne's model
! -------------------------------------------------------
@@ -226,7 +226,7 @@ CONTAINS
DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1)
! ! Potential nitrogen fixation dependant on temperature and iron
ztemp = ts(ji,jj,jk,jp_tem,Kmm)
- zmudia = MAX( 0.,-0.001096*ztemp**2 + 0.057*ztemp -0.637 ) * 7.625
+ zmudia = MAX( 0.,-0.001096*ztemp**2 + 0.057*ztemp -0.637 ) / rno3
! Potential nitrogen fixation dependant on temperature and iron
xdianh4 = tr(ji,jj,jk,jpnh4,Kbb) / ( concnnh4 + tr(ji,jj,jk,jpnh4,Kbb) )
xdiano3 = tr(ji,jj,jk,jpno3,Kbb) / ( concnno3 + tr(ji,jj,jk,jpno3,Kbb) ) * (1. - xdianh4)
diff --git a/src/TOP/PISCES/SED/sedsfc.F90 b/src/TOP/PISCES/SED/sedsfc.F90
index 3c86dc99e9af35718bd2e61a1502963139c4ba3c..460d760d94900d1ae159b9ec11b90db666977c74 100644
--- a/src/TOP/PISCES/SED/sedsfc.F90
+++ b/src/TOP/PISCES/SED/sedsfc.F90
@@ -5,6 +5,7 @@ MODULE sedsfc
!!=====================================================================
!! * Modules used
USE sed ! sediment global variable
+ USE sedini
USE sedarr
USE seddta
@@ -53,11 +54,11 @@ CONTAINS
IF ( tmask(ji,jj,ikt) == 1 ) THEN
tr(ji,jj,ikt,jptal,Kbb) = trc_data(ji,jj,1)
tr(ji,jj,ikt,jpdic,Kbb) = trc_data(ji,jj,2)
- tr(ji,jj,ikt,jpno3,Kbb) = trc_data(ji,jj,3) * 7.625
- tr(ji,jj,ikt,jppo4,Kbb) = trc_data(ji,jj,4) * 122.
+ tr(ji,jj,ikt,jpno3,Kbb) = trc_data(ji,jj,3) * redC / redNo3
+ tr(ji,jj,ikt,jppo4,Kbb) = trc_data(ji,jj,4) * redC
tr(ji,jj,ikt,jpoxy,Kbb) = trc_data(ji,jj,5)
tr(ji,jj,ikt,jpsil,Kbb) = trc_data(ji,jj,6)
- tr(ji,jj,ikt,jpnh4,Kbb) = trc_data(ji,jj,7) * 7.625
+ tr(ji,jj,ikt,jpnh4,Kbb) = trc_data(ji,jj,7) * redC / redNo3
tr(ji,jj,ikt,jpfer,Kbb) = trc_data(ji,jj,8)
ENDIF
END_2D