diff --git a/src/drivers/0d/fabm0d.F90 b/src/drivers/0d/fabm0d.F90
index 55062748..30b02115 100644
--- a/src/drivers/0d/fabm0d.F90
+++ b/src/drivers/0d/fabm0d.F90
@@ -14,7 +14,7 @@ module fabm_0d
 ! !USES:
    use time
    use input
-   use eqstate,only:rho_feistel
+   use density,only:get_rho
 
    use fabm
    use fabm_driver
@@ -616,7 +616,7 @@ subroutine update_environment(n)
       if (model_type==1) column_depth = mixed_layer_depth%value
 
       ! Compute density from temperature and salinity, if required by biogeochemistry.
-      if (compute_density) dens = rho_feistel(salt%value,temp%value,0.5_rk*column_depth,.true.)
+      if (compute_density) dens = get_rho(salt%value,temp%value,0.5_rk*column_depth)
    end subroutine update_environment
 
    subroutine update_light()