diff --git a/src/oce_ale_pressure_bv.F90 b/src/oce_ale_pressure_bv.F90
index 4a2a19082..c6689a2eb 100644
--- a/src/oce_ale_pressure_bv.F90
+++ b/src/oce_ale_pressure_bv.F90
@@ -3176,7 +3176,6 @@ subroutine init_ref_density_advanced(tracers, partit, mesh)
         do nz=nzmin,nzmax
             call densityJM_components(ref_temp1D(nz), ref_salt1D(nz), bulk_0, bulk_pz, bulk_pz2, rhopot)
             rho = bulk_0   + auxz*bulk_pz   + auxz*bulk_pz2
-            density_ref(nzmin, node) = rho*rhopot/(rho+0.1_WP*auxz)
             auxz=Z_3d_n(nz,node)
             rho = bulk_0   + auxz*bulk_pz   + auxz*bulk_pz2
             density_ref(nz,node) = rho*rhopot/(rho+0.1_WP*auxz)
@@ -3184,4 +3183,4 @@ subroutine init_ref_density_advanced(tracers, partit, mesh)
     end do
 !$OMP END PARALLEL DO
 if(mype==0) write(*,*) ' --> compute reference density'
-end subroutine init_ref_density_advanced
\ No newline at end of file
+end subroutine init_ref_density_advanced