junction_init_part2 Subroutine

public subroutine junction_init_part2(me, krn, cfg)

Arguments

Type IntentOptional Attributes Name
type(junction_t), intent(inout) :: me
type(krn_t), intent(inout), target :: krn
class(input_t), intent(in), pointer :: cfg

Calls

proc~~junction_init_part2~~CallsGraph proc~junction_init_part2 junction_init_part2 interface~to_ptr_arr to_ptr_arr proc~junction_init_part2->interface~to_ptr_arr proc~all_junction_allocation all_junction_allocation proc~junction_init_part2->proc~all_junction_allocation proc~convert2int_port convert2int_port proc~junction_init_part2->proc~convert2int_port proc~input_dbl input_t%input_dbl proc~junction_init_part2->proc~input_dbl proc~input_dict1d input_t%input_dict1d proc~junction_init_part2->proc~input_dict1d proc~input_str input_t%input_str proc~junction_init_part2->proc~input_str proc~krn_coo_add_link krn_t%krn_coo_add_link proc~junction_init_part2->proc~krn_coo_add_link proc~port_list_find port_list_t%port_list_find proc~junction_init_part2->proc~port_list_find proc~to_ptr_arr0d to_ptr_arr0d interface~to_ptr_arr->proc~to_ptr_arr0d proc~to_ptr_arr1d to_ptr_arr1d interface~to_ptr_arr->proc~to_ptr_arr1d proc~to_ptr_arr2d to_ptr_arr2d interface~to_ptr_arr->proc~to_ptr_arr2d proc~continue_if continue_if proc~input_dbl->proc~continue_if proc~get0d_hdf5 get0d_hdf5 proc~input_dbl->proc~get0d_hdf5 proc~input_has_key input_t%input_has_key proc~input_dbl->proc~input_has_key proc~input_node input_t%input_node proc~input_dbl->proc~input_node to_real to_real proc~input_dbl->to_real get_string get_string proc~input_dict1d->get_string proc~input_dict1d->proc~input_has_key proc~input_dict1d->proc~input_node proc~select_dict select_dict proc~input_dict1d->proc~select_dict proc~select_list select_list proc~input_dict1d->proc~select_list proc~input_str->proc~input_has_key proc~input_str->proc~input_node proc~select_scalar select_scalar proc~input_str->proc~select_scalar interface~to_str to_str proc~port_list_find->interface~to_str proc~to_int to_int proc~port_list_find->proc~to_int proc~str_from_int str_from_int interface~to_str->proc~str_from_int proc~str_from_real str_from_real interface~to_str->proc~str_from_real proc~get0d_hdf5->proc~input_str proc~get0d_hdf5->proc~input_has_key h5close_f h5close_f proc~get0d_hdf5->h5close_f h5dclose_f h5dclose_f proc~get0d_hdf5->h5dclose_f h5dopen_f h5dopen_f proc~get0d_hdf5->h5dopen_f h5dread_f h5dread_f proc~get0d_hdf5->h5dread_f h5fclose_f h5fclose_f proc~get0d_hdf5->h5fclose_f h5fopen_f h5fopen_f proc~get0d_hdf5->h5fopen_f h5ltget_dataset_info_f h5ltget_dataset_info_f proc~get0d_hdf5->h5ltget_dataset_info_f h5open_f h5open_f proc~get0d_hdf5->h5open_f proc~input_has_key->proc~select_dict get get proc~input_has_key->get proc~input_node->proc~select_dict proc~input_node->get proc~dictionary_traverse dictionary_traverse proc~input_node->proc~dictionary_traverse proc~dictionary_traverse->proc~dictionary_traverse get_dictionary get_dictionary proc~dictionary_traverse->get_dictionary proc~abort_if abort_if proc~dictionary_traverse->proc~abort_if

Called by

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

Source Code

subroutine junction_init_part2(me,krn,cfg)
    type(junction_t), intent(inout) :: me
    type(krn_t), target, intent(inout) :: krn
    class(input_t), pointer, intent(in) :: cfg

    type(input_t), allocatable :: branches(:)
    integer :: i,j

    call all_junction_allocation(me)

    branches = cfg%dict1d('link')
    do i = 1, me%NbTotBr
        me%br(i)%p => krn%FF_flux_ports(krn%FF_flux_list%find(branches(i)%str('id'),&
                                                              convert2int_port(branches(i)%str('node'))))
        if(me%br(i)%p%connected) then
            print*,'Port ',branches(i)%str('node'),' belonging to ',branches(i)%str('id'),' already connected to a link'
            stop
        endif
        me%br(i)%p%connected=.true.
        me%br(i)%angle = branches(i)%dbl('angle',0.0_dp)
    enddo

    me%kappa=branches(1)%dbl('kappa',1.0_dp)

    do j=1,me%NbTotBr
      do i=1,me%NbTotBr
        if(j/=i) then
            call krn%coo_add_link(me%br(j)%p, me%br(i)%p, idx_4x4_col, &
                    idx_4x4_row, to_ptr_arr(me%Matrix(:,:,i,j)))
        endif
      enddo
    enddo
end subroutine junction_init_part2