--- /dev/null
+ subroutine sfcofa ( bilan, nbbasc, libasc,
+ > lenoeu, larete,
+ > nufade, nufafi, nbvoto,
+ > coonoe,
+ > somare, filare, np2are,
+ > facare,
+ > hettri, aretri,
+ > voltri,
+ > hetqua, arequa, filqua,
+ > volqua,
+ > ulsort, langue, codret)
+c ______________________________________________________________________
+c H O M A R D
+c
+c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
+c
+c Version originale enregistree le 18 juin 1996 sous le numero 96036
+c aupres des huissiers de justice Simart et Lavoir a Clamart
+c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
+c aupres des huissiers de justice
+c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
+c
+c HOMARD est une marque deposee d'Electricite de France
+c
+c Copyright EDF 1996
+c Copyright EDF 1998
+c Copyright EDF 2002
+c Copyright EDF 2020
+c ______________________________________________________________________
+c
+c Suivi de Frontiere - COntroles des FAces
+c - - -- --
+c ______________________________________________________________________
+c . . . . .
+c . nom . e/s . taille . description .
+c .____________________________________________________________________.
+c . bilan . s . 1 . bilan du controle de l'arete .
+c . . . . 0 : pas de probleme .
+c . . . . 1 : probleme .
+c . nbbasc . s . 1 . nombre de bascule a faire .
+c . libasc . s . * . liste des aretes a basculer .
+c . lenoeu . e . 1 . noeud qui bouge .
+c . larete . e . 1 . arete a controler .
+c . nufade . e . 1 . numero face depart des voisines de larete .
+c . nufafi . e . 1 . numero face fin des voisines de larete .
+c . nbvoto . e . 1 . nombre de volumes total .
+c . coonoe . e . nbnoto . coordonnees des noeuds .
+c . . . *sdim . .
+c . somare . e .2*nbarto. numeros des extremites d'arete .
+c . filare . e . nbarto . premiere fille des aretes .
+c . np2are . e . nbarto . noeud milieux des aretes .
+c . facare . es . nbfaar . liste des faces contenant une arete .
+c . hettri . e . nbtrto . historique de l'etat des triangles .
+c . aretri . e .nbtrto*3. numeros des 3 aretes des triangles .
+c . voltri . e .2*nbtrto. numeros des 2 volumes par triangle .
+c . . . . voltri(i,k) definit le i-eme voisin de k .
+c . . . . 0 : pas de voisin .
+c . . . . j>0 : tetraedre j .
+c . . . . j<0 : pyramide/pentaedre dans pypetr(1/2,j).
+c . hetqua . e . nbquto . historique de l'etat des quadrangles .
+c . arequa . e .nbquto*4. numeros des 4 aretes des quadrangles .
+c . filqua . e . nbquto . premier fils des quadrangles .
+c . volqua . e .2*nbquto. numeros des 2 volumes par quadrangle .
+c . . . . volqua(i,k) definit le i-eme voisin de k .
+c . . . . 0 : pas de voisin .
+c . . . . j>0 : hexaedre j .
+c . . . . j<0 : pyramide/pentaedre dans pypequ(1/2,j).
+c . ulsort . e . 1 . numero d'unite logique de la liste standard.
+c . langue . e . 1 . langue des messages .
+c . . . . 1 : francais, 2 : anglais .
+c . codret . es . 1 . code de retour des modules .
+c . . . . 0 : pas de probleme .
+c . . . . x : probleme .
+c ______________________________________________________________________
+c
+c====
+c 0. declarations et dimensionnement
+c====
+c
+c 0.1. ==> generalites
+c
+ implicit none
+ save
+c
+ character*6 nompro
+ parameter ( nompro = 'SFCOFA' )
+c
+#include "nblang.h"
+c
+c 0.2. ==> communs
+c
+#include "envex1.h"
+c
+#include "envca1.h"
+#include "nombno.h"
+#include "nombar.h"
+#include "nombqu.h"
+#include "nombtr.h"
+#include "impr02.h"
+c
+c 0.3. ==> arguments
+c
+ integer bilan, nbbasc, libasc(*)
+ integer lenoeu, larete
+ integer nufade, nufafi, nbvoto
+c
+ integer somare(2,nbarto), filare(nbarto), np2are(nbarto)
+ integer facare(nbfaar)
+ integer hettri(nbtrto), aretri(nbtrto,3)
+ integer voltri(2,nbtrto)
+ integer hetqua(nbquto), arequa(nbquto,4), filqua(nbquto)
+ integer volqua(2,nbquto)
+c
+ double precision coonoe(nbnoto,sdim)
+c
+ integer ulsort, langue, codret
+c
+c 0.4. ==> variables locales
+c
+ integer iaux
+ integer laface
+c
+ logical bascul
+c
+ integer nbmess
+ parameter ( nbmess = 10 )
+ character*80 texte(nblang,nbmess)
+c
+c 0.5. ==> initialisations
+c ______________________________________________________________________
+c
+c====
+c 1. messages
+c====
+c
+#include "impr01.h"
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,1)) 'Entree', nompro
+ call dmflsh (iaux)
+#endif
+c
+ texte(1,4) = '(/,''.. Examen du '',a,i10)'
+ texte(1,5) = '(''.. Probleme.'')'
+ texte(1,6) = '(''.. Bascule a faire.'')'
+c
+ texte(2,4) = '(/,''.. Examination of '',a,'' # '',i10)'
+ texte(2,5) = '(''. Problem.'')'
+ texte(2,6) = '(''.. Swapping.'')'
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,4)) mess14(langue,1,1), larete
+ write (ulsort,texte(langue,4)) mess14(langue,1,-1), lenoeu
+#endif
+c
+ codret = 0
+c
+ bilan = 0
+ nbbasc = 0
+c
+c====
+c 2. boucle sur les faces s'appuyant sur l'arete
+c On ne s'interesse qu'aux aretes qui viennent d'etre decoupees et
+c qui font partie d'une frontiere reconnue
+c On ne s'interesse qu'aux faces qui ne bordent aucun volume
+c====
+c
+ do 21 , iaux = nufade, nufafi
+c
+ if ( codret.eq.0 ) then
+c
+ bascul = .false.
+c
+ laface = facare(iaux)
+c
+c 2.1. ==> si la face voisine est un triangle
+c
+ if ( laface.gt.0 ) then
+cgn write (ulsort,*)'.. Face voisine : triangle ', laface
+c
+ if ( nbvoto.eq.0 .or. voltri(1,laface).eq.0 ) then
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,3)) 'SFTQTR', nompro
+#endif
+ call sftqtr ( bilan, bascul,
+ > lenoeu, larete, laface,
+ > coonoe,
+ > somare, filare, np2are,
+ > hettri, aretri,
+ > ulsort, langue, codret)
+c
+ endif
+c
+c 2.2. ==> si la face voisine est un quadrangle
+c
+ else
+cgn write(ulsort,*)'.. Face voisine : quadrangle ',-laface
+c
+ if ( nbvoto.eq.0 .or. volqua(1,-laface).eq.0 ) then
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,3)) 'SFTQQU', nompro
+#endif
+ call sftqqu ( bilan,
+ > lenoeu, larete, -laface,
+ > coonoe,
+ > somare, filare, np2are,
+ > hetqua, arequa, filqua,
+ > ulsort, langue, codret)
+c
+ endif
+c
+ endif
+c
+c 2.3. ==> Memorisations
+c
+ if ( codret.eq.0 ) then
+c
+ if ( bilan.ne.0 ) then
+ goto 30
+ endif
+c
+ if ( bascul ) then
+ nbbasc = nbbasc + 1
+ libasc(nbbasc) = laface
+ endif
+c
+ endif
+c
+ endif
+c
+ 21 continue
+c
+c====
+c 3. Bilan
+c====
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,*) '3. Bilan ; codret = ', codret
+#endif
+c
+ 30 continue
+c
+#ifdef _DEBUG_HOMARD_
+ if ( codret.eq.0 ) then
+ if ( bilan.ne.0 ) then
+ write (ulsort,texte(langue,4)) mess14(langue,1,1), larete
+ write (ulsort,texte(langue,5))
+ endif
+ if ( bascul ) then
+ write (ulsort,texte(langue,4)) mess14(langue,1,1), larete
+ write (ulsort,texte(langue,6))
+ endif
+ endif
+#endif
+c
+c====
+c 4. La fin
+c====
+c
+ if ( codret.ne.0 ) then
+c
+#include "envex2.h"
+c
+ write (ulsort,texte(langue,1)) 'Sortie', nompro
+ write (ulsort,texte(langue,2)) codret
+c
+ endif
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,1)) 'Sortie', nompro
+ call dmflsh (iaux)
+#endif
+c
+ end