init_ustboco_ghostface.f90

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

!------------------------------------------------------------------------------!
! Procedure : init_ustboco_ghostface      Auteur : J. Gressier
!                                         Date   : Mars 2003
! Fonction                                Modif  :
!   Affectation des connectivites entre faces limites et cellules limites
!   pour le type "ghostface" (point fictif sur la face) : cela revient a
!   avoir une cellule fictive de volume nul.
!
! Defauts/Limitations/Divers :
!
!------------------------------------------------------------------------------!
subroutine init_ustboco_ghostface(ib, defboco, ust_mesh)

use TYPHMAKE
!use VARCOM
use OUTPUT
use USTMESH
use MENU_BOCO

implicit none

! -- Declaration des entrees --
integer        :: ib                     ! numero de condition aux limites
type(mnu_boco) :: defboco                ! parametres du solveur

! -- Declaration des entrees/sorties --
type(st_ustmesh) :: ust_mesh             ! maillage et connectivites

! -- Declaration des sorties --

! -- Declaration des variables internes --
integer :: if                            ! indice boco et de face
integer :: icell, iface                  ! index de cellule et de face

! -- Debut de la procedure --

! affectation de connectivite face limites -> cellules fictives
! la variable ust_mesh%ncell_lim contient le nombre courant de cellules limites affectees

! le tableau de cellules est cense pouvoir contenir le nombre de cellules fictives (test)

if ((ust_mesh%ncell_lim+ust_mesh%boco(ib)%nface)>(ust_mesh%ncell-ust_mesh%ncell_int)) then
  call erreur("Allocation","Pas assez de cellules allouees pour les cellules fictives")
endif

! -- boucle sur la liste des faces de la condition limite --

do if = 1, ust_mesh%boco(ib)%nface    
  
  ! affectation de connectivite face limites -> cellules fictives
  ust_mesh%ncell_lim = ust_mesh%ncell_lim + 1     ! nouvelle cellule limite
  icell = ust_mesh%ncell_int + ust_mesh%ncell_lim ! index de cellule limite
  iface = ust_mesh%boco(ib)%iface(if)             ! index de face

  ! definition geometrique de la cellule fictive
  ust_mesh%mesh%volume(icell,1,1) = 0._krp
  ust_mesh%mesh%centre(icell,1,1) = ust_mesh%mesh%iface(iface,1,1)%centre
  
  if (ust_mesh%facecell%fils(iface,2) == 0) then
    ust_mesh%facecell%fils(iface,2) = icell        ! affectation de la cellule fictive
  else
    call erreur("Initialisation de connectivite", &
                "Connectivite deja affectee sur face limite")
  endif

enddo

endsubroutine init_ustboco_ghostface

!------------------------------------------------------------------------------!
! Historique des modifications
!
! mars 2003 : creation de la procedure
!------------------------------------------------------------------------------!