diff --git a/src/single_interfaces_initial_conditions.jl b/src/single_interfaces_initial_conditions.jl
index c7ab305..7153859 100644
--- a/src/single_interfaces_initial_conditions.jl
+++ b/src/single_interfaces_initial_conditions.jl
@@ -23,7 +23,7 @@ function STSingleInterfaceInitialConditions(model, depth_of_interface, salinity,
                                             interface_smoothing = NoSmoothing,
                                             background_state = NoBackground)
 
-    eos = model.buoyancy.model.equation_of_state
+    eos = model.buoyancy.formulation.equation_of_state
 
     R_ρ = compute_R_ρ(salinity, temperature, depth_of_interface, eos)
 
diff --git a/src/staircase_initial_conditions.jl b/src/staircase_initial_conditions.jl
index ddfd001..1e9c78f 100644
--- a/src/staircase_initial_conditions.jl
+++ b/src/staircase_initial_conditions.jl
@@ -18,7 +18,7 @@ struct STStaircaseInitialConditions{T} <: AbstractStaircaseInitialConditions
 end
 function STStaircaseInitialConditions(model, number_of_interfaces, depth_of_interfaces, salinity, temperature)
 
-    eos = model.buoyancy.model.equation_of_state
+    eos = model.buoyancy.formulation.equation_of_state
 
     R_ρ = compute_R_ρ(salinity, temperature, depth_of_interfaces, eos)