initialize_type_output Subroutine

    subroutine initialize_type_output(self, input, domain, coordinate)
        implicit none
        class(type_output), intent(inout) :: self
        type(type_input), intent(in) :: input
        class(type_domain), intent(inout), optional :: domain
        type(type_dp_3d), intent(inout), pointer :: coordinate

        character(256) :: dir_path
        logical :: exists

        character(:), allocatable :: command
        integer(int32) :: i, j, idx
        integer(int32) :: total
        real(real64) :: tmp_xi, tmp_eta
        logical :: is_inside

        integer(int32) :: iObs, iElem
        integer(int32) :: nElements
        integer(int32) :: local_id, local_type
        integer(int32) :: ierr

        character(8) :: output_extentions(3) = [".dat", ".csv", ".log"]
        character(8) :: output_file_extentions(5) = [".dat", ".csv", ".vtk", ".vtu", ".log"]

        ! Path settings
        dir_path = get_project_path()

        self%dir_output = trim(adjustl(dir_path))//"Output/"
        call setup_directory(self%dir_output, output_extentions)
        self%dir_output_field = trim(adjustl(dir_path))//"Output/Files/"
        call setup_directory(self%dir_output_field, output_file_extentions)

        self%log_file_name = trim(adjustl(self%dir_output))//"run.log"

        self%is_thermal = input%basic%analysis_controls%calculate_thermal
        self%is_hydraulic = input%basic%analysis_controls%calculate_hydraulic

        if (allocated(self%observations)) deallocate (self%observations)
        allocate (self%observations(size(input%output_settings%history_output%variable_names)))
        do i = 1, size(input%output_settings%history_output%variable_names)
            call self%observations(i)%initialize(input, coordinate, domain, self%dir_output, &
                                                 input%output_settings%history_output%variable_names(i))
            call self%observations(i)%write_header(input%output_settings%history_output%output_interval_unit)
        end do

        call self%overall%initialize(input, coordinate, domain, self%dir_output_field)

    end subroutine initialize_type_output