diff --git a/src/OCE/stp2d.F90 b/src/OCE/stp2d.F90
index 0cea844f1e239832c563fa82418d3c6e1ee5658b..6ca69903dcca2841cf9242ab4d1f4160028b0911 100644
--- a/src/OCE/stp2d.F90
+++ b/src/OCE/stp2d.F90
@@ -19,10 +19,8 @@ MODULE stp2d
    !!----------------------------------------------------------------------
    USE step_oce       ! time stepping used modules
    USE domqco         ! quasi-eulerian coordinate      (dom_qco_r3c routine)
-   USE dynadv_cen2    ! centred flux form advection    (dyn_adv_cen2 routine)
-   USE dynadv_ubs     ! UBS flux form advection        (dyn_adv_ubs  routine)
-   USE dynkeg         ! kinetic energy gradient        (dyn_keg      routine)
    USE dynspg_ts      ! 2D mode integration
+   USE sshwzv         ! vertical speed
    USE sbc_ice , ONLY : snwice_mass, snwice_mass_b
    USE sbcapr         ! surface boundary condition: atmospheric pressure
    USE sbcwave,  ONLY : bhd_wave
@@ -111,6 +109,10 @@ CONTAINS
       !
       !                             !*  compute advection + coriolis *!
       !
+      r3t(:,:,Kaa) =  ssh(:,:,Kaa) * r1_ht_0(:,:)                           ! ratio at t-point at Kaa (n-1)
+      !
+      CALL wzv    ( kt, Kbb, Kbb, Kaa , uu(:,:,:,Kbb), vv(:,:,:,Kbb), ww )  ! ww guess at Kbb (n)
+      !
       CALL dyn_adv( kt, Kbb, Kbb      , uu, vv, Krhs)       !- vector form KEG+ZAD 
       !                                                     !- flux   form ADV
       CALL dyn_vor( kt,            Kbb, uu, vv, Krhs )      !- vector form COR+RVO