heat_mat_ass_boundary.f90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
! Copyright (c) 2019 FrontISTR Commons
! This software is released under the MIT License, see LICENSE.txt
!-------------------------------------------------------------------------------
module m_heat_mat_ass_boundary
contains
 
  subroutine heat_mat_ass_boundary(hecMESH, hecMAT, hecMESHmpc, hecMATmpc, &
     & fstrSOLID, fstrHEAT, next_time, delta_time)
    use m_fstr
    use m_heat_mat_ass_bc_cflux
    use m_heat_mat_ass_bc_dflux
    use m_heat_mat_ass_bc_fixt
    use m_heat_mat_ass_bc_film
    use m_heat_mat_ass_bc_radiate
    use hecmw_ebc_defer
    implicit none
    type(fstr_solid) :: fstrSOLID
    type(fstr_heat) :: fstrHEAT
    type(hecmwst_matrix) :: hecMAT
    type(hecmwst_matrix), pointer :: hecMATmpc
    type(hecmwst_local_mesh) :: hecMESH
    type(hecmwst_local_mesh), pointer :: hecMESHmpc
    type(hecmwst_ebc) :: hecEBC
    real(kind=kreal) :: next_time, delta_time, beta
 
    beta = fstrheat%beta
 
    call heat_mat_ass_bc_cflux(hecmat, fstrheat, next_time, delta_time, beta)
 
    call heat_mat_ass_bc_dflux(hecmesh, hecmat, fstrsolid, fstrheat, next_time, delta_time, beta)
 
    call heat_mat_ass_bc_film(hecmesh, hecmat, fstrsolid, fstrheat, next_time, delta_time, beta)
 
    call heat_mat_ass_bc_radiate(hecmesh, hecmat, fstrheat, next_time, delta_time, beta)
 
    call hecmw_ebc_init(hecmat, hecebc)
    call heat_mat_ass_bc_fixt(hecmat, hecebc, fstrheat, next_time, delta_time, beta)
 
    call hecmw_mpc_mat_ass(hecmesh, hecmat, hecmeshmpc, hecmatmpc)
    call hecmw_mpc_trans_rhs(hecmesh, hecmat, hecmatmpc)
    call hecmw_ebc_apply(hecmeshmpc, hecmatmpc, hecebc)
    call hecmw_ebc_finalize(hecebc)
 
  end subroutine heat_mat_ass_boundary
end module m_heat_mat_ass_boundary