look_for_cep Subroutine

private subroutine look_for_cep(model, z0, zi, Pc, Vc, Tc, a, u, found, CEP)

Uses

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in) :: model

Equation of state model
real(kind=pr), intent(in) :: z0(:)

Molar fractions of the first fluid
real(kind=pr), intent(in) :: zi(:)

Molar fractions of the second fluid
real(kind=pr), intent(in) :: Pc

Pressure [bar]
real(kind=pr), intent(in) :: Vc

Volume [L/mol]
real(kind=pr), intent(in) :: Tc

Temperature [K]
real(kind=pr), intent(in) :: a

Molar fraction of the second fluid
real(kind=pr), intent(inout) :: u(:)

Eigen-vector
logical, intent(out) :: found

Found a Critical End Point
type(EquilibriumState), intent(out) :: CEP

Critical End Point

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: Fcep(size(z0)+4)
real(kind=pr), private :: V_cep
real(kind=pr), private :: Xcep(size(z0)+4)
real(kind=pr), private :: dFcep(size(z0)+4,size(z0)+4)
real(kind=pr), private :: dXcep(size(z0)+4)
integer, private :: its
real(kind=pr), private :: y_cep(size(z0))