ftdss.F90 Source File

module Main_FTDSS
    use, intrinsic :: iso_fortran_env
    use :: stdlib_logger
    use :: module_core
    use :: Inout_Input
    use :: module_control, only:type_time, type_iteration
    use :: module_output, only:type_output
    use :: module_domain, only:type_domain
    use :: module_properties, only:type_proereties_manager
    use :: module_boundary, only:type_bc
    use :: module_initial, only:type_ic

    use :: Main_Thermal
    implicit none

    type :: type_ftdss
        ! type(Type_Input) :: Input

        type(type_dp_3d), pointer :: coordinate
        type(type_domain) :: domain
        ! type(Belonging), allocatable :: NodeBelonging(:)
        class(abst_thermal), allocatable :: Thermal

        type(type_proereties_manager) :: property
        type(type_bc) :: bc
        type(type_ic) :: ic

        type(type_variable) :: phi

        type(type_time) :: time
        type(type_iteration) :: Iteration
        type(Type_output) :: output

    contains
        procedure, pass(self) :: initialize => initialize_type_ftdss
    end type type_ftdss

contains
    subroutine initialize_type_ftdss(self)
        implicit none
        class(type_ftdss), intent(inout) :: self

        type(type_input) :: input

        integer(int32) :: nsize
        integer(int32) :: iN

        integer(int32) :: ierr

        character(len=10), allocatable :: profiler_labels(:)

        ! ★ 計測したいセクション名を定義
        profiler_labels = [character(len=10) :: "IO", "Setup", "Assemble", "Solve", "Total"]
        call self%time%initialize(profiler_sections=profiler_labels)
        call self%time%Record("Start")
        call self%time%Profile_Start("Total")
        call self%time%Profile_Start("IO")

        ! Initialize the FDTSS module
        ! This is where you would set up any necessary parameters or configurations
        call input%initialize()
        call self%time%initialize(input=input)

        if (input%output_settings%standard_output%print_progress) then
            call global_logger%configure(level=information_level, &
                                         time_stamp=.true., &
                                         max_width=0)
        else
            call global_logger%configure(level=warning_level, &
                                         time_stamp=.false., &
                                         max_width=0)
        end if
        call setup_handler()

        !---------------------------------------------------------------------------------------------------------------------------
        !
        !---------------------------------------------------------------------------------------------------------------------------
        nsize = input%geometry%vtk%num_points

        ! Initialize the Structure
        allocate (self%coordinate)
        call self%coordinate%initialize(nsize)
        self%coordinate = input%geometry%vtk%POINTS

        call self%domain%initialize(input, self%coordinate, ierr)
        if (ierr /= 0) then
            print *, "Error initializing domain in Type_Thermal_3Phase_2D_Construct"
            return
        end if

        call self%bc%initialize(input, self%domain)
        call self%ic%initialize(input)

        call global_logger%log_information(message="Boundary and Initial Conditions set up.")

        self%Thermal = Type_Thermal_3Phase_2D(input, self%coordinate, self%domain)

        call self%property%initialize(input, ierr)

        call self%output%initialize(input, self%domain, self%coordinate)

        call self%phi%initialize(nsize, input%basic%solver_settings%bdf_order)
        call self%ic%apply('porosity', self%domain, self%phi)

        call self%output%output_coloring(self%domain)

        call self%time%Profile_Stop("IO")
        call global_logger%log_information(message="FTDSS module initialized successfully.")
    end subroutine initialize_type_ftdss

end module Main_FTDSS