assemble_type_thermal_3phase_2d Module Subroutine

    module subroutine assemble_type_thermal_3phase_2d(self, domain, property, porosity, dt, step, iter)
        implicit none
        class(type_thermal_3phase_2d), intent(inout) :: self
        type(type_domain), intent(inout) :: domain
        type(type_proereties_manager), intent(inout) :: property
        real(real64), intent(in) :: porosity(:)
        real(real64), intent(in) :: dt
        integer(int32), intent(in) :: step
        integer(int32), intent(in) :: iter

        ! if (step >= 2) then
        !     self%CT_old(2)%val(:) = self%CT_old(1)%val(:)
        !     self%CT_old(1)%val(:) = self%CT_l%val(:)
        ! end if

        self%CT_l%val(:) = 0.0d0
        self%KT_l%val(:) = 0.0d0
        self%KT_star_0%val(:) = 0.0d0
        ! ! if (step == 1) then

        ! ! end if
        !---------------------------------------------------------------------------------------------------------------------------
        ! 各剛性行列の組み立て
        !---------------------------------------------------------------------------------------------------------------------------
        !!!-----------------------------------------------------------------------
        ! (A, domain, Temperature, porosity, Propeties)
        call Assemble_Mass_Heat_1_Parallel(self%CT_l, domain, self%T%new, porosity, property)
        ! call Assemble_Mass_Heat_1(self%CT_l, domain, self%T%new, porosity, property)
        !
        call Assemble_Diffusion_Heat_1_Parallel(self%KT_l, domain, self%T%new, porosity, property)
        ! call Assemble_Diffusion_Heat_1(self%KT_l, domain, self%T%new, porosity, property)
        !---------------------------------------------------------------------------------------------------------------------------
        ! if (step == 1) then
        self%KT_star_0 = dt * self%KT_l + self%CT_l

        !     if (iter == 1) then
        ! self%CT_old(1)%val(:) = self%CT_l%val(:)
        !         self%KT_old%val(:) = self%KT_l%val(:)
        !         self%PHIT(:) = 0.0d0
        ! self%PHIT_old(:) = -self%CT_old(1) * self%T%old(:, 1)
        !     end if
        ! print *, size(self%T%old(:, 1))
        self%PHIT(:) = self%CT_l * self%T%old(:, 1)
        ! stop
        ! self%PHIT(:) = -self%CT_old(1) * self%T%old(:, 1)
        ! self%PHIT(:) = dt * (self%KT_l * self%T%pre(:)) + self%CT_l * self%T%pre(:) + self%PHIT_old(:)
        ! else
        !     self%KT_star_0%val(:) = 2.0d0 * dt * self%KT_l%val(:) + 3.0d0 * self%CT_l%val(:)
        !     if (iter == 1) then
        !         self%PHIT_old(:) = -4.0d0 * (self%CT_old(1) * self%T%old(:, 1)) + self%CT_old(2) * self%T%old(:, 2)
        !     end if
        !     self%PHIT(:) = 2.0d0 * dt * (self%KT_l * self%T%pre(:)) + 3.0d0 * (self%CT_l * self%T%pre(:)) + self%PHIT_old(:)
        ! end if
        ! self%PHIT(:) = self%PHIT_old(:)
    end subroutine assemble_type_thermal_3phase_2d