all_junction_allocation Subroutine

public subroutine all_junction_allocation(me)

Arguments

Type IntentOptional Attributes Name
type(junction_t), intent(inout) :: me

Called by

proc~~all_junction_allocation~~CalledByGraph proc~all_junction_allocation all_junction_allocation proc~junction_init_part2 junction_init_part2 proc~junction_init_part2->proc~all_junction_allocation program~reims_p reims_p program~reims_p->proc~junction_init_part2

Source Code

subroutine all_junction_allocation(me)
    type(junction_t), intent(inout) :: me
    integer :: j

    allocate(me%br(me%NbTotBr))
    allocate(me%FlJ(me%NbTotBr))
    allocate(me%Matrix(Nb_VarC,Nb_VarC,me%NbTotBr,me%NbTotBr))
    me%Matrix = 0.0_dp
    allocate(me%DerFlx_dBr(Nb_VarC,Nb_VarC,me%NbTotBr,me%NbTotBr))
    me%DerFlx_dBr = 0.0_dp
    allocate(me%idxKappaRef(me%NbTotBr))
    me%idxKappaRef = .false.; me%idxKappaRef(1) = .true.
    allocate(me%DerVel_dBr(Nb_VarC,me%NbTotBr,me%NbTotBr))
    me%DerVel_dBr=0.0_dp
    
    allocate(me%dynamic%br(me%NbTotBr))
    allocate(me%dynamic%p_star_prev(me%NbTotBr))
    me%dynamic%p_star_prev = 0.0_dp
    allocate(me%dynamic%alpha_prev(me%NbTotBr))
    me%dynamic%alpha_prev = 1.0e-1_dp

    do j=1,me%NbTotBr
        allocate(me%FlJ(j)%DerFlx(Nb_VarC,Nb_VarC,me%NbTotBr))
        me%FlJ(j)%DerFlx = 0.0_dp
        allocate(me%FlJ(j)%DerVit(Nb_VarC,me%NbTotBr))
        me%FlJ(j)%DerVit = 0.0_dp
    enddo
end subroutine all_junction_allocation