pt_F_NP Subroutine

public subroutine pt_F_NP(model, z, np, beta_w, X, ns, S, F, df)

Function to solve at each point of a multi-phase envelope.

Arguments

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

Model to use.

real(kind=pr), intent(in) :: z(:)

Mixture global composition.

integer, intent(in) :: np

Number of main phases.

real(kind=pr), intent(in) :: beta_w

Fraction of the reference (incipient) phase.

real(kind=pr), intent(in) :: X(:)

Vector of variables.

integer, intent(in) :: ns

Number of specification.

real(kind=pr), intent(in) :: S

Specification value.

real(kind=pr), intent(out) :: F(size(X))

Vector of functions valuated.

real(kind=pr), intent(out) :: df(size(X),size(X))

Jacobian matrix.


Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: K(np,size(z))
real(kind=pr), private :: P
real(kind=pr), private :: T
real(kind=pr), private :: Vl(np)
real(kind=pr), private :: Vw
real(kind=pr), private :: betas(np)
real(kind=pr), private :: denom(size(z))
real(kind=pr), private :: denomdlnK(np,size(z),size(z))
real(kind=pr), private, dimension(size(z), size(z)) :: dlnphi_dn
real(kind=pr), private, dimension(np, size(z), size(z)) :: dlnphi_dn_l
real(kind=pr), private, dimension(size(z), size(z)) :: dlnphi_dn_w
real(kind=pr), private :: dlnphi_dp(size(z))
real(kind=pr), private, dimension(np, size(z)) :: dlnphi_dp_l
real(kind=pr), private, dimension(size(z)) :: dlnphi_dp_w
real(kind=pr), private :: dlnphi_dt(size(z))
real(kind=pr), private, dimension(np, size(z)) :: dlnphi_dt_l
real(kind=pr), private, dimension(size(z)) :: dlnphi_dt_w
real(kind=pr), private :: dwdb(np,size(z))
real(kind=pr), private :: dwdlnK(np,size(z))
real(kind=pr), private :: dx_l_dlnK(np,np,size(z))
integer, private :: i
integer, private :: idx_1
integer, private :: idx_2
integer, private :: j
integer, private :: l
integer, private :: lb
real(kind=pr), private :: lnphi(size(z))
real(kind=pr), private, dimension(np, size(z)) :: lnphi_l
real(kind=pr), private, dimension(size(z)) :: lnphi_w
real(kind=pr), private :: moles(size(z))
integer, private :: nc
integer, private :: phase
integer, private :: ub
real(kind=pr), private, dimension(size(z)) :: w
real(kind=pr), private, dimension(np, size(z)) :: x_l