Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(type_input), | intent(inout) | :: | input | |||
type(type_dp_3d), | intent(inout), | pointer | :: | coordinate | ||
type(type_domain), | intent(inout) | :: | domain |
module function construct_type_thermal_3phase_2d(input, coordinate, domain) result(structure) implicit none class(abst_thermal), allocatable :: structure type(type_input), intent(inout) :: input type(type_dp_3d), intent(inout), pointer :: coordinate type(type_domain), intent(inout) :: domain integer(int32) :: i integer(int32) :: nNode integer(int32) :: ierr if (allocated(structure)) deallocate (structure) allocate (Type_Thermal_3Phase_2D :: structure) nNode = domain%get_num_nodes() call structure%KT_star_0%initialize(domain) structure%KT_l = structure%KT_star_0 structure%KT_old = structure%KT_star_0 structure%CT_l = structure%KT_star_0 structure%order = input%basic%solver_settings%bdf_order allocate (structure%CT_old(structure%order)) do i = 1, structure%order structure%CT_old(i) = structure%KT_star_0 end do call allocate_array(structure%FT, nNode) call allocate_array(structure%FT_old, nNode) call allocate_array(structure%PHIT, nNode) call allocate_array(structure%PHIT_old, nNode) call structure%T%initialize(nNode, structure%order) !--------------------------------------------------------------------------------------------------------------------------- ! 線形求解ソルバーの設定 !--------------------------------------------------------------------------------------------------------------------------- structure%solver = create_solver(input, "thermal", structure%KT_star_0, nNode) !--------------------------------------------------------------------------------------------------------------------------- end function construct_type_thermal_3phase_2d