diff --git a/src/methods/property_solvers/multicomponent/bubble_point.jl b/src/methods/property_solvers/multicomponent/bubble_point.jl index b83840a31..035a39484 100755 --- a/src/methods/property_solvers/multicomponent/bubble_point.jl +++ b/src/methods/property_solvers/multicomponent/bubble_point.jl @@ -294,7 +294,7 @@ end ###Bubble Temperature function __x0_bubble_temperature(model::EoSModel,p,x,Tx0 = nothing,volatiles = FillArrays.Fill(true,length(model)),pure = split_model(model),crit = nothing) - sat = extended_saturation_temperature.(pure,p,volatiles,crit_retry = false) + sat = extended_saturation_temperature.(pure,p,crit,volatiles,crit_retry = false) if crit === nothing _crit = __crit_pure.(sat,pure,volatiles) else diff --git a/src/methods/property_solvers/multicomponent/dew_point.jl b/src/methods/property_solvers/multicomponent/dew_point.jl index b57cba5a2..0930f65fa 100755 --- a/src/methods/property_solvers/multicomponent/dew_point.jl +++ b/src/methods/property_solvers/multicomponent/dew_point.jl @@ -133,7 +133,7 @@ end function __x0_dew_temperature(model::EoSModel,p,y,Tx0 = nothing,condensables = FillArrays.Fill(true,length(model)),pure = split_model(model),crit = nothing) - sat = extended_saturation_temperature.(pure,p,condensables) + sat = extended_saturation_temperature.(pure,p,crit,condensables) if crit === nothing _crit = __crit_pure.(sat,pure,condensables) else