FSlink_init_part2 Subroutine

public subroutine FSlink_init_part2(me, krn, cfg)

Arguments

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

Calls

proc~~fslink_init_part2~~CallsGraph proc~fslink_init_part2 FSlink_init_part2 interface~to_ptr_arr to_ptr_arr proc~fslink_init_part2->interface~to_ptr_arr proc~array_conv array_conv proc~fslink_init_part2->proc~array_conv proc~check_if_conversion_required check_if_conversion_required proc~fslink_init_part2->proc~check_if_conversion_required proc~input_dbl input_t%input_dbl proc~fslink_init_part2->proc~input_dbl proc~input_dbl1d input_t%input_dbl1d proc~fslink_init_part2->proc~input_dbl1d proc~input_dict1d input_t%input_dict1d proc~fslink_init_part2->proc~input_dict1d proc~input_has_key input_t%input_has_key proc~fslink_init_part2->proc~input_has_key proc~input_int1d input_t%input_int1d proc~fslink_init_part2->proc~input_int1d proc~input_str input_t%input_str proc~fslink_init_part2->proc~input_str proc~input_str1d input_t%input_str1d proc~fslink_init_part2->proc~input_str1d proc~krn_coo_add_link krn_t%krn_coo_add_link proc~fslink_init_part2->proc~krn_coo_add_link proc~nusselt_init nusselt_init proc~fslink_init_part2->proc~nusselt_init proc~port_list_find port_list_t%port_list_find proc~fslink_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~count_chars count_chars proc~check_if_conversion_required->proc~count_chars proc~input_dbl->proc~input_has_key proc~continue_if continue_if proc~input_dbl->proc~continue_if proc~get0d_hdf5 get0d_hdf5 proc~input_dbl->proc~get0d_hdf5 proc~input_node input_t%input_node proc~input_dbl->proc~input_node to_real to_real proc~input_dbl->to_real proc~input_dbl1d->proc~input_has_key interface~to_str to_str proc~input_dbl1d->interface~to_str proc~input_dbl1d->proc~continue_if proc~get1d_hdf5 get1d_hdf5 proc~input_dbl1d->proc~get1d_hdf5 proc~input_dbl1d->proc~input_node proc~select_scalar select_scalar proc~input_dbl1d->proc~select_scalar proc~input_dbl1d->to_real proc~input_dict1d->proc~input_has_key get_string get_string proc~input_dict1d->get_string 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 get get proc~input_has_key->get proc~input_has_key->proc~select_dict proc~input_int1d->proc~input_has_key proc~input_int1d->interface~to_str proc~input_int1d->proc~continue_if proc~get1d_range get1d_range proc~input_int1d->proc~get1d_range proc~input_int1d->proc~input_node proc~input_int1d->proc~select_scalar to_integer to_integer proc~input_int1d->to_integer proc~input_str->proc~input_has_key proc~input_str->proc~input_node proc~input_str->proc~select_scalar proc~input_str1d->proc~input_has_key proc~input_str1d->proc~input_node proc~input_str1d->proc~select_scalar proc~nusselt_init->proc~input_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_has_key proc~get0d_hdf5->proc~input_str 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~get1d_hdf5->proc~input_has_key proc~get1d_hdf5->proc~input_str proc~get1d_hdf5->h5dclose_f proc~get1d_hdf5->h5dopen_f proc~get1d_hdf5->h5dread_f proc~get1d_hdf5->h5fclose_f proc~get1d_hdf5->h5fopen_f proc~get1d_hdf5->h5ltget_dataset_info_f proc~get1d_range->proc~input_has_key proc~input_node->get proc~input_node->proc~select_dict 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~~fslink_init_part2~~CalledByGraph proc~fslink_init_part2 FSlink_init_part2 program~reims_p reims_p program~reims_p->proc~fslink_init_part2

Source Code

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

    type(input_t), allocatable :: links(:)
    real(dp), allocatable :: distances(:), surfaces(:)
    integer, allocatable :: nodes1(:), nodes2(:), conv(:)
    type(str_ptr), allocatable :: labels(:)
    integer :: i,offset1,offset2,max_value,j
    character(:), allocatable :: name2
    logical :: conversion_required

    links = cfg%dict1d('link')
    nodes1=links(1)%int1d('node')
    offset1 = krn%FS_list%find(links(1)%str('id'), 1) - 1

    if(links(2)%has_key('node')) then
        allocate(me%lk(2,me%sizeNodes))
        allocate(me%DerSrcP_SedT(me%sizeNodes),me%DerSrcP_SrdT(me%sizeNodes),me%DerSrcSdUi(Nb_VarC,me%sizeNodes))
        allocate(me%MatrixPS(2,me%sizeNodes),me%MatrixSP(Nb_VarC,me%sizeNodes),me%SurfCont(me%sizeNodes))
        me%MatrixPS=0.0_dp; me%MatrixSP=0.0_dp; me%SurfCont=0.0_dp; me%DerSrcP_SedT=0.0_dp; me%DerSrcP_SrdT=0.0_dp
        me%DerSrcSdUi=0.0_dp

        do i = 1, me%sizeNodes
            me%lk(1,i)%p => krn%FS_ports(offset1 + nodes1(i))
        enddo

        nodes2=links(2)%int1d('node')
        name2=links(2)%str('id')
        offset2 = krn%FS_list%find(name2, 1) - 1
    
        call check_if_conversion_required(krn%PortIdxForFSLink_List,name2,conversion_required,max_value)
        if(conversion_required) then
            allocate(conv(max_value))
            call array_conv(krn%PortIdxForFSLink_List,name2,conv)
            do i = 1, me%sizeNodes
                me%lk(2,i)%p => krn%FS_ports(offset2 + conv(nodes2(i)))
            enddo
        else
            do i = 1, me%sizeNodes
                me%lk(2,i)%p => krn%FS_ports(offset2 + nodes2(i))
            enddo
        endif
        
        if(me%lk(2,1)%p%typeS=="solid") then
            distances=links(2)%dbl1d('distance', n=me%sizeNodes)
            surfaces=cfg%dbl1d('contact_surface', n=me%sizeNodes)
            do i = 1, me%sizeNodes
              me%lk(2,i)%p%Dist4Grad=distances(i)
              me%SurfCont(i)=surfaces(i)
            enddo
        else if(me%lk(2,1)%p%typeS=="strand") then
            me%WetPer=cfg%dbl('wetted_perimeter')
        endif
        
        call nusselt_init(me%nuss,cfg) ! nusselt required in the link  
        if(.NOT.associated(me%nuss)) error stop 'Nusselt does not exist in thermalink'       
              
        do i = 1, me%sizeNodes
            call krn%coo_add_link(me%lk(2,i)%p, me%lk(1,i)%p, idx_1x4_col, idx_1x4_row, to_ptr_arr(me%MatrixSP(:,i)))
            call krn%coo_add_link(me%lk(1,i)%p, me%lk(2,i)%p, idx_2x1_col, idx_2x1_row, to_ptr_arr(me%MatrixPS(:,i)))
        enddo

    else if(links(2)%has_key('label')) then

        allocate(me%lk(maxval(me%nb_2D_Ports(:))+1,me%sizeNodes))
        allocate(me%DerSrcP_SedT2D(maxval(me%nb_2D_Ports(:)),me%sizeNodes))
        allocate(me%DerSrcP_SrdT2D(maxval(me%nb_2D_Ports(:)),me%sizeNodes))
        allocate(me%derFlx2DdUi(Nb_VarC,maxval(me%nb_2D_Ports(:)),me%sizeNodes))
        allocate(me%MatrixPS_2D(2,maxval(me%nb_2D_Ports(:)),me%sizeNodes))
        allocate(me%MatrixSP_2D(Nb_VarC,maxval(me%nb_2D_Ports(:)),me%sizeNodes))

        labels=links(2)%str1d('label')
        do i = 1, me%sizeNodes
            me%lk(1,i)%p => krn%FS_ports(offset1 + nodes1(i))
            do j=1,me%nb_2D_Ports(i)
              me%lk(1+j,i)%p => krn%FS_ports(krn%FS_list%find(links(2)%str('id')//"+"//trim(labels(i)%p),j))
            enddo
        enddo

        me%Ht=cfg%dbl('heat_coefficient')    

        do i = 1, me%sizeNodes
            do j=1,me%nb_2D_Ports(i)
                call krn%coo_add_link(me%lk(1+j,i)%p, me%lk(1,i)%p, idx_1x4_col, idx_1x4_row, to_ptr_arr(me%MatrixSP_2D(:,j,i)))
                call krn%coo_add_link(me%lk(1,i)%p, me%lk(1+j,i)%p, idx_2x1_col, idx_2x1_row, to_ptr_arr(me%MatrixPS_2D(:,j,i)))
            enddo
        enddo        

    endif

end subroutine FSlink_init_part2