fill_g_Dreg_terms Subroutine

private elemental subroutine fill_g_Dreg_terms(T, ThAr)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: T
type(ThArrays_t), intent(inout) :: ThAr

Called by

proc~~fill_g_dreg_terms~~CalledByGraph proc~fill_g_dreg_terms fill_g_Dreg_terms proc~eos_pt_terms eos_pT_terms proc~eos_pt_terms->proc~fill_g_dreg_terms proc~fill_e_tpart fill_e_Tpart proc~fill_e_tpart->proc~fill_g_dreg_terms proc~r_rot r_roT proc~r_rot->proc~fill_g_dreg_terms proc~t_rop T_roP proc~t_rop->proc~fill_g_dreg_terms proc~circulator_dynamic_parameters_df circulator_dynamic_parameters_t%circulator_dynamic_parameters_df proc~circulator_dynamic_parameters_df->proc~t_rop proc~jacobian_rot jacobian_roT proc~circulator_dynamic_parameters_df->proc~jacobian_rot proc~derivative_function_of_pstar_pt derivative_function_of_pstar_PT proc~derivative_function_of_pstar_pt->proc~t_rop proc~derivative_function_of_pstar_pt->proc~jacobian_rot proc~droeint_drop droeint_droP proc~droeint_drop->proc~t_rop proc~eos_terms eos_terms proc~droeint_drop->proc~eos_terms proc~eos_e_terms eos_e_terms proc~eos_e_terms->proc~fill_e_tpart proc~ffsrclink_prelax_resolution_from_and_to_ports FFsrcLink_Prelax_resolution_from_and_to_ports proc~ffsrclink_prelax_resolution_from_and_to_ports->proc~t_rop proc~he_prop he_prop proc~ffsrclink_prelax_resolution_from_and_to_ports->proc~he_prop proc~state_rot state_roT proc~ffsrclink_prelax_resolution_from_and_to_ports->proc~state_rot proc~funct_ro_from_pt ro_from_pT_obj_t%funct_ro_from_pT proc~funct_ro_from_pt->proc~eos_pt_terms proc~funct_t_from_rop T_from_roP_obj_t%funct_T_from_roP proc~funct_t_from_rop->proc~eos_pt_terms proc~function_of_pstar_mt function_of_pstar_MT proc~function_of_pstar_mt->proc~t_rop proc~outgoing_branch_cold_circulator outgoing_branch_cold_circulator proc~outgoing_branch_cold_circulator->proc~t_rop proc~outgoing_branch_cold_circulator->proc~circulator_dynamic_parameters_df proc~state_rop state_roP proc~outgoing_branch_cold_circulator->proc~state_rop proc~circulator_dynamic_parameters_f circulator_dynamic_parameters_t%circulator_dynamic_parameters_f proc~outgoing_branch_cold_circulator->proc~circulator_dynamic_parameters_f proc~dc2_rot dc2_roT proc~outgoing_branch_cold_circulator->proc~dc2_rot proc~outgoing_branch_cold_circulator->proc~jacobian_rot proc~jacobian_star_req_rem jacobian_star_Req_rem proc~outgoing_branch_cold_circulator->proc~jacobian_star_req_rem proc~solve_boundary_mt solve_boundary_MT proc~solve_boundary_mt->proc~t_rop proc~solve_boundary_mt->proc~function_of_pstar_mt proc~solve_boundary_mt->proc~state_rop proc~solve_boundary_mt->proc~dc2_rot proc~solve_boundary_mt->proc~jacobian_rot proc~jacobian_star_req_rem~2 jacobian_star_Req_rem proc~solve_boundary_mt->proc~jacobian_star_req_rem~2 proc~solve_junction solve_junction proc~solve_junction->proc~r_rot proc~solve_junction->proc~droeint_drop proc~solve_junction->proc~state_rop proc~t_roe T_roE proc~solve_junction->proc~t_roe proc~solve_junction->proc~jacobian_rot proc~jacobian_star_req_rem_loc jacobian_star_Req_rem_loc proc~solve_junction->proc~jacobian_star_req_rem_loc proc~state_roe_withr state_roE_withR proc~state_roe_withr->proc~t_rop proc~state_roe_withr->proc~jacobian_rot proc~state_rop->proc~t_rop proc~state_rop->proc~jacobian_rot proc~state_rop_withr state_roP_withR proc~state_rop_withr->proc~t_rop proc~state_rop_withr->proc~jacobian_rot proc~t_roe->proc~fill_e_tpart proc~boundary_resolution_from_and_to_ports boundary_resolution_from_and_to_ports proc~boundary_resolution_from_and_to_ports->proc~solve_boundary_mt proc~solve_boundary_pt solve_boundary_PT proc~boundary_resolution_from_and_to_ports->proc~solve_boundary_pt proc~circulator_dynamic_parameters_f->proc~state_rop proc~circulator_resolution_from_and_to_ports circulator_resolution_from_and_to_ports proc~circulator_resolution_from_and_to_ports->proc~outgoing_branch_cold_circulator proc~incoming_branch_cold_circulator incoming_branch_cold_circulator proc~circulator_resolution_from_and_to_ports->proc~incoming_branch_cold_circulator proc~mdot_derivatives mdot_derivatives proc~circulator_resolution_from_and_to_ports->proc~mdot_derivatives proc~eos_terms->proc~eos_e_terms proc~friction_source_term friction_source_term proc~friction_source_term->proc~state_roe_withr proc~state_roe state_roE proc~friction_source_term->proc~state_roe proc~friction_source_term->proc~jacobian_rot proc~friction_source_term->proc~he_prop proc~from_sol_to_prim_channel from_sol_to_prim_channel proc~from_sol_to_prim_channel->proc~state_roe_withr proc~from_sol_to_prim_channel->proc~state_roe proc~funct_t_from_roe T_from_roE_obj_t%funct_T_from_roE proc~funct_t_from_roe->proc~eos_e_terms proc~function_of_pstar_pt function_of_pstar_PT proc~function_of_pstar_pt->proc~state_rop proc~hug_brent_f hug_brent_t%hug_brent_f proc~hug_brent_f->proc~droeint_drop proc~hugoniot_star_state hugoniot_star_state proc~hugoniot_star_state->proc~droeint_drop proc~incoming_branch_cold_circulator->proc~state_rop proc~incoming_branch_cold_circulator->proc~jacobian_rot proc~hhh_derivatives HHH_derivatives proc~incoming_branch_cold_circulator->proc~hhh_derivatives proc~incoming_branch_cold_circulator->proc~jacobian_star_req_rem proc~junction_dynamic_parameters_fcn junction_dynamic_parameters_t%junction_dynamic_parameters_fcn proc~junction_dynamic_parameters_fcn->proc~droeint_drop proc~junction_resolution_from_and_to_ports Junction_resolution_from_and_to_ports proc~junction_resolution_from_and_to_ports->proc~solve_junction proc~main_loop main_loop proc~main_loop->proc~ffsrclink_prelax_resolution_from_and_to_ports proc~main_loop->proc~boundary_resolution_from_and_to_ports proc~main_loop->proc~circulator_resolution_from_and_to_ports proc~main_loop->proc~friction_source_term proc~main_loop->proc~from_sol_to_prim_channel proc~main_loop->proc~junction_resolution_from_and_to_ports proc~source_term_definition_channel source_term_definition_channel proc~main_loop->proc~source_term_definition_channel proc~full_physics_definition_channel full_physics_definition_channel proc~main_loop->proc~full_physics_definition_channel proc~he_riemann_solver_channel He_Riemann_solver_channel proc~main_loop->proc~he_riemann_solver_channel proc~last_tasks_for_channels last_tasks_for_channels proc~main_loop->proc~last_tasks_for_channels proc~solve_boundary_pt->proc~derivative_function_of_pstar_pt proc~solve_boundary_pt->proc~state_rop proc~solve_boundary_pt->proc~state_rop_withr proc~solve_boundary_pt->proc~function_of_pstar_pt proc~solve_boundary_pt->proc~hugoniot_star_state proc~solve_boundary_pt->proc~dc2_rot proc~solve_boundary_pt->proc~jacobian_rot proc~solve_boundary_pt->proc~jacobian_star_req_rem~2 proc~solve_boundary_pt->proc~state_rot proc~sourceterms_notfriction_channel SourceTerms_NOTfriction_channel proc~sourceterms_notfriction_channel->proc~state_roe_withr proc~sourceterms_notfriction_channel->proc~state_roe proc~sourceterms_notfriction_channel->proc~jacobian_rot proc~sourceterms_notfriction_channel->proc~he_prop proc~state_roe->proc~t_roe proc~state_roe->proc~jacobian_rot proc~dc2_rot->proc~eos_terms proc~jacobian_rot->proc~eos_terms proc~source_term_definition_channel->proc~sourceterms_notfriction_channel program~reims_p reims_p program~reims_p->proc~main_loop proc~channel_init_part1 channel_init_part1 program~reims_p->proc~channel_init_part1 proc~full_physics_definition_channel->proc~dc2_rot proc~full_physics_definition_channel->proc~source_term_definition_channel proc~he_prop->proc~jacobian_rot proc~hhh_derivatives->proc~jacobian_rot proc~jacobian_channel_req_rem jacobian_channel_Req_rem proc~jacobian_channel_req_rem->proc~jacobian_rot proc~jacobian_star_req_rem->proc~jacobian_rot proc~jacobian_star_req_rem_loc->proc~jacobian_rot proc~jacobian_star_req_rem~2->proc~jacobian_rot proc~mdot_derivatives->proc~jacobian_rot proc~state_rot->proc~jacobian_rot proc~channel_init_part1->proc~state_rot proc~he_riemann_solver_channel->proc~jacobian_channel_req_rem proc~last_tasks_for_channels->proc~state_rot

Source Code

elemental subroutine fill_g_Dreg_terms(T, ThAr)
    real(dp), intent(in)  :: T
    type(ThArrays_t), intent(inout) :: ThAr

    ThAr%g(1) =nn(1) *T+nn(2)*sqrt(T)+nn(3)+nn(4)/T+nn(5)/(T**2)
    ThAr%g(2) =nn(6) *T+nn(7)+nn(8)/T+nn(9)/(T**2)
    ThAr%g(3) =nn(10)*T+nn(11)+nn(12)/T
    ThAr%g(4) =nn(13)
    ThAr%g(5) =nn(14)/T+nn(15)/(T**2)
    ThAr%g(6) =nn(16)/T
    ThAr%g(7) =nn(17)/T+nn(18)/(T**2)
    ThAr%g(8) =nn(19)/(T**2)
    ThAr%g(9) =nn(20)/(T**2)+nn(21)/(T**3)
    ThAr%g(10)=nn(22)/(T**2)+nn(23)/(T**4)
    ThAr%g(11)=nn(24)/(T**2)+nn(25)/(T**3)
    ThAr%g(12)=nn(26)/(T**2)+nn(27)/(T**4)
    ThAr%g(13)=nn(28)/(T**2)+nn(29)/(T**3)
    ThAr%g(14)=nn(30)/(T**2)+nn(31)/(T**3)+nn(32)/(T**4)

    ThAr%Derg(1) =nn(1)+nn(2)/(2.0_dp*sqrt(T))-nn(4)/(T**2)-2.0_dp*nn(5)/(T**3)
    ThAr%Derg(2) =nn(6)-nn(8)/(T**2)-2.0_dp*nn(9)/(T**3)
    ThAr%Derg(3) =nn(10)-nn(12)/(T**2)
    ThAr%Derg(4) =0.0_dp
    ThAr%Derg(5) =-nn(14)/(T**2)-2.0_dp*nn(15)/(T**3)
    ThAr%Derg(6) =-nn(16)/(T**2)
    ThAr%Derg(7) =-nn(17)/(T**2)-2.0_dp*nn(18)/(T**3)
    ThAr%Derg(8) =-2.0_dp*nn(19)/(T**3)
    ThAr%Derg(9) =-2.0_dp*nn(20)/(T**3)-3.0_dp*nn(21)/(T**4)
    ThAr%Derg(10)=-2.0_dp*nn(22)/(T**3)-4.0_dp*nn(23)/(T**5)
    ThAr%Derg(11)=-2.0_dp*nn(24)/(T**3)-3.0_dp*nn(25)/(T**4)
    ThAr%Derg(12)=-2.0_dp*nn(26)/(T**3)-4.0_dp*nn(27)/(T**5)
    ThAr%Derg(13)=-2.0_dp*nn(28)/(T**3)-3.0_dp*nn(29)/(T**4)
    ThAr%Derg(14)=-2.0_dp*nn(30)/(T**3)-3.0_dp*nn(31)/(T**4)-4.0_dp*nn(32)/(T**5)
end subroutine fill_g_Dreg_terms