fluxlimite.f90

国外大名顶顶的“台风”并行计算流体力学CFD软件的早期版本的源代码

!------------------------------------------------------------------------------!
! Procedure : flux_limite                 Auteur : J. Gressier/E. Radenac
!                                         Date   : Juin 2004
! Fonction                                Modif  : (cf Historique)
!  Flux aux faces limites quand necessaire
!
! Defauts/Limitations/Divers :
!
!------------------------------------------------------------------------------!
subroutine fluxlimite(defsolver, domaine, flux)

use TYPHMAKE
use OUTPUT
use VARCOM
use MENU_SOLVER
use MENU_BOCO
use USTMESH
use DEFFIELD 

implicit none

! -- Declaration des entrees --
type(mnu_solver) :: defsolver        ! type d'equation a resoudre
type(st_ustmesh) :: domaine          ! domaine non structure a integrer

! -- Declaration des sorties --
type(st_genericfield)   :: flux        ! flux physiques

! -- Declaration des variables internes --
integer          :: ifb, if, ib, idef     ! index de liste, index de face limite et parametres

! -- Debut de la procedure --

do ib = 1, domaine%nboco
  idef = domaine%boco(ib)%idefboco
  if ( (defsolver%boco(idef)%typ_boco == bc_wall_adiab).or. &
       (defsolver%boco(idef)%typ_boco == bc_wall_flux) .or. &
       (defsolver%boco(idef)%typ_boco == bc_wall_hconv) ) then
    do ifb = 1, domaine%boco(ib)%nface
      if = domaine%boco(ib)%iface(ifb)
      flux%tabscal(1)%scal(if) = domaine%boco(ib)%bocofield%tabscal(1)%scal(ifb)
    enddo
  endif
enddo

endsubroutine fluxlimite

!------------------------------------------------------------------------------!
! Historique des modifications
!
! juin 2004 : creation de la procedure
!------------------------------------------------------------------------------!