ro_pT Function

public function ro_pT(p, T) result(ro)

Solve p(ro,T) = p for ro.

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: p
real(kind=dp), intent(in) :: T

Return Value real(kind=dp)


Calls

proc~~ro_pt~~CallsGraph proc~ro_pt ro_pT proc~brent brent proc~ro_pt->proc~brent f f proc~brent->f proc~set_error set_error proc~brent->proc~set_error proc~zero zero proc~brent->proc~zero proc~zero->f

Called by

proc~~ro_pt~~CalledByGraph proc~ro_pt ro_pT proc~channel_init_part1 channel_init_part1 proc~channel_init_part1->proc~ro_pt proc~solve_boundary_pt solve_boundary_PT proc~solve_boundary_pt->proc~ro_pt proc~boundary_resolution_from_and_to_ports boundary_resolution_from_and_to_ports proc~boundary_resolution_from_and_to_ports->proc~solve_boundary_pt program~reims_p reims_p program~reims_p->proc~channel_init_part1 proc~main_loop main_loop program~reims_p->proc~main_loop proc~main_loop->proc~boundary_resolution_from_and_to_ports

Source Code

function ro_pT(p, T) result(ro)
    !! Solve p(ro,T) = p for ro.
    real(dp), intent(in) :: p, T
    real(dp) :: ro
    type(ro_from_pT_obj_t) :: obj

    obj%T = T
    obj%p = p*1.0e-6_dp
    ro = brent(obj, ro_from_pT_min_Arp, ro_from_pT_max_Arp, 'ro_pT')
    ro = ro*MolMass*1000.0_dp
end function ro_pT