tx_F_NP Subroutine

public subroutine tx_F_NP(model, z0, zi, np, P, beta_w, kinds_x, kind_w, X, ns, S, F, df)

Uses

  • proc~~tx_f_np~~UsesGraph proc~tx_f_np tx_F_NP iso_fortran_env iso_fortran_env proc~tx_f_np->iso_fortran_env

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

Arguments

Type IntentOptional Attributes Name
class(ArModel), intent(in) :: model
real(kind=pr), intent(in) :: z0(:)
real(kind=pr), intent(in) :: zi(:)
integer, intent(in) :: np

Number of main phases.
real(kind=pr), intent(in) :: P

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

Fraction of the reference (incipient) phase.
character(len=14), intent(in) :: kinds_x(np)

Kinds of the main phases.
character(len=14), intent(in) :: kind_w

Kind of the reference 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.

Calls

proc~~tx_f_np~~CallsGraph proc~tx_f_np tx_F_NP proc~get_z get_z proc~tx_f_np->proc~get_z proc~lnphi_pt ArModel%lnphi_pt proc~tx_f_np->proc~lnphi_pt proc~lnphi_vt ArModel%lnphi_vt proc~lnphi_pt->proc~lnphi_vt proc~volume~3 ArModel%volume proc~lnphi_pt->proc~volume~3 residual_helmholtz residual_helmholtz proc~lnphi_vt->residual_helmholtz get_v0 get_v0 proc~volume~3->get_v0 interface~newton newton proc~volume~3->interface~newton proc~newton_1d newton_1d interface~newton->proc~newton_1d

Called by

proc~~tx_f_np~~CalledByGraph proc~tx_f_np tx_F_NP proc~solve_point~7 TXEnvelMP%solve_point proc~solve_point~7->proc~tx_f_np proc~tx_envelope tx_envelope proc~tx_envelope->proc~solve_point~7

Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: K(np,size(z0))
real(kind=pr), private :: T
real(kind=pr), private :: Vl(np)
real(kind=pr), private :: Vw
real(kind=pr), private :: alpha
real(kind=pr), private :: betas(np)
real(kind=pr), private :: denom(size(z0))
real(kind=pr), private :: denomdlnK(np,size(z0),size(z0))
real(kind=pr), private, dimension(size(z0), size(z0)) :: dlnphi_dn
real(kind=pr), private, dimension(np, size(z0), size(z0)) :: dlnphi_dn_l
real(kind=pr), private, dimension(size(z0), size(z0)) :: dlnphi_dn_w
real(kind=pr), private :: dlnphi_dt(size(z0))
real(kind=pr), private, dimension(np, size(z0)) :: dlnphi_dt_l
real(kind=pr), private, dimension(size(z0)) :: dlnphi_dt_w
real(kind=pr), private :: dwda(size(z0))
real(kind=pr), private :: dwdb(np,size(z0))
real(kind=pr), private :: dwdlnK(np,size(z0))
real(kind=pr), private :: dx_l_dlnK(np,np,size(z0))
real(kind=pr), private :: dzda(size(z0))
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(z0))
real(kind=pr), private, dimension(np, size(z0)) :: lnphi_l
real(kind=pr), private, dimension(size(z0)) :: lnphi_w
real(kind=pr), private :: moles(size(z0))
integer, private :: nc
integer, private :: phase
integer, private :: ub
real(kind=pr), private, dimension(size(z0)) :: w
real(kind=pr), private, dimension(np, size(z0)) :: x_l
real(kind=pr), private :: z(size(z0))