tx_envelope Function

public function tx_envelope(model, z0, zi, np, P, x_l0, w0, betas0, T0, alpha0, ns0, dS0, beta_w, points)

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
real(kind=pr), intent(in) :: P
real(kind=pr), intent(in) :: x_l0(np,size(z0))
real(kind=pr), intent(in) :: w0(size(z0))
real(kind=pr), intent(in) :: betas0(np)
real(kind=pr), intent(in) :: T0
real(kind=pr), intent(in) :: alpha0
integer, intent(in) :: ns0
real(kind=pr), intent(in) :: dS0
real(kind=pr), intent(in) :: beta_w

Fraction of the reference (incipient) phase.

integer, intent(in), optional :: points

Return Value type(TXEnvelMP)


Variables

Type Visibility Attributes Name Initial
real(kind=pr), private :: F(size(z0)*np+np+2)
real(kind=pr), private :: S

Specified value

real(kind=pr), private :: T
real(kind=pr), private :: X(size(z0)*np+np+2)
real(kind=pr), private :: X0(size(X))

Initial guess for the point

real(kind=pr), private :: alpha
real(kind=pr), private, allocatable :: alphas(:)
real(kind=pr), private :: betas(np)
real(kind=pr), private :: dF(size(z0)*np+np+2,size(z0)*np+np+2)
real(kind=pr), private :: dS

Step size of the specification for the next point

real(kind=pr), private :: dX(size(z0)*np+np+2)
real(kind=pr), private :: dXdS(size(z0)*np+np+2)
type(MPPoint), private, allocatable :: env_points(:)
integer, private :: i

Point calculation index

integer, private :: inner

Number of times a failed point is retried to converge

integer, private :: its
integer, private :: lb

Lower bound, index of the first component of a phase

integer, private :: max_iterations = 100
integer, private :: nc
integer, private :: ns

Number of the specified variable

integer, private :: number_of_points
type(MPPoint), private :: point
integer, private :: ub

Upper bound, index of the last component of a phase

real(kind=pr), private :: w(size(z0))
real(kind=pr), private :: x_l(np,size(z0))