--- /dev/null
+ subroutine cmcp1b ( nulofa, lepent,
+ > aretri, nivtri,
+ > filqua,
+ > facpen, cofape,
+ > niveau,
+ > trifad, cotrvo, areqtr,
+ > ulsort, langue, codret )
+c
+c ______________________________________________________________________
+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 Creation du Maillage - Conformite - decoupage des Pentaedres
+c - - - -
+c - cas 1, phase B
+c - -
+c Reperage des aretes et des triangles sur les faces externes
+c ______________________________________________________________________
+c . . . . .
+c . nom . e/s . taille . description .
+c .____________________________________________________________________.
+c . nulofa . e . 4 . numero local des faces a traiter .
+c . lepent . e . 1 . pentaedre a decouper .
+c . aretri . e .nouvtr*3. numeros des 3 aretes des triangles .
+c . nivtri . e . nouvtr . niveau des triangles .
+c . filqua . e . nouvqu . premier fils des quadrangles .
+c . facpen . e .nouvpf*5. numeros des faces des pentaedres .
+c . cofape . e .nouvpf*5. codes des faces des pentaedres .
+c . niveau . s . 1 . niveau des faces issus du decoupage .
+c . trifad . s .(2,0:2) . triangles traces sur les faces decoupees .
+c . cotrvo . s .(2,0:2) . code des triangles dans les volumes .
+c . areqtr . s . (2,2) . aretes tri tracees sur les faces decoupees .
+c . ulsort . e . 1 . unite logique de la sortie generale .
+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 . . . . 1 : aucune arete ne correspond .
+c ______________________________________________________________________
+c
+c====
+c 0. declarations et dimensionnement
+c====
+c
+c 0.1. ==> generalites
+c
+ implicit none
+ save
+c
+ character*6 nompro
+ parameter ( nompro = 'CMCP1B' )
+c
+#include "nblang.h"
+c
+c 0.2. ==> communs
+c
+#include "envex1.h"
+c
+#include "nouvnb.h"
+c
+c 0.3. ==> arguments
+c
+ integer lepent, nulofa(2)
+ integer aretri(nouvtr,3), nivtri(nouvtr)
+ integer filqua(nouvqu)
+ integer facpen(nouvpf,5), cofape(nouvpf,5)
+ integer niveau
+ integer trifad(2,0:2), cotrvo(2,0:2), areqtr(2,2)
+c
+ integer ulsort, langue, codret
+c
+c 0.4. ==> variables locales
+c
+ integer iaux, jaux
+c
+ integer nbmess
+ parameter ( nbmess = 10 )
+ character*80 texte(nblang,nbmess)
+c
+c 0.5. ==> initialisations
+c ______________________________________________________________________
+c
+c====
+c 1. initialisations
+c====
+c
+#include "impr01.h"
+c
+#ifdef _DEBUG_HOMARD_
+ write (ulsort,texte(langue,1)) 'Entree', nompro
+ call dmflsh (iaux)
+#endif
+c
+#include "impr03.h"
+#include "impr04.h"
+c
+ codret = 0
+c
+c====
+c 2. Triangles et aretes tracees sur le quadrangle dont le
+c triangle central sera la face 2 de la pyramide
+c trifad(1,0) : triangle central de ce decoupage
+c trifad(1,1) : triangle du cote de la face F1 du pentaedre
+c trifad(1,2) : triangle du cote de la face F2 du pentaedre
+c cotrvo(1,0/1/2) : futur code du triangle trifad(p,0/1/2) dans la
+c description des fils
+c areqtr(1,1) : arete interne au quadrangle de bord et bordant le
+c triangle trifad(p,1)
+c areqtr(1,2) : arete interne au quadrangle de bord et bordant le
+c triangle trifad(p,2)
+c====
+c
+ iaux = facpen(lepent,nulofa(1))
+ jaux = cofape(lepent,nulofa(1))
+ trifad(1,0) = -filqua(iaux)
+ if ( jaux.lt.5 ) then
+ cotrvo(1,0) = 2
+ trifad(1,1) = trifad(1,0) + 2
+ cotrvo(1,1) = 1
+ trifad(1,2) = trifad(1,0) + 1
+ cotrvo(1,2) = 4
+ areqtr(1,1) = aretri(trifad(1,0),3)
+ areqtr(1,2) = aretri(trifad(1,0),1)
+ else
+ cotrvo(1,0) = 4
+ trifad(1,1) = trifad(1,0) + 1
+ cotrvo(1,1) = 4
+ trifad(1,2) = trifad(1,0) + 2
+ cotrvo(1,2) = 1
+ areqtr(1,1) = aretri(trifad(1,0),1)
+ areqtr(1,2) = aretri(trifad(1,0),3)
+ endif
+#ifdef _DEBUG_HOMARD_
+ write(ulsort,90006) 'Quadrangle = ', iaux,', code =', jaux
+ do 2222 , iaux = 0, 2
+ write (ulsort,90015) 'trifad(1,0/1/2) =', trifad(1,iaux),
+ > ', aretes', (aretri(trifad(1,iaux),jaux),jaux=1,3)
+ 2222 continue
+ write(ulsort,90006) 'cotrvo(1,0) = ', cotrvo(1,0),
+ > 'cotrvo(1,1) = ', cotrvo(1,1),
+ > 'cotrvo(1,2) = ', cotrvo(1,2)
+ write(ulsort,90002) 'areqtr(1,1) = ', areqtr(1,1)
+ write(ulsort,90002) 'areqtr(1,2) = ', areqtr(1,2)
+#endif
+c
+c====
+c 3. Triangles et aretes tracees sur le quadrangle dont le
+c triangle central sera la face 4 de la pyramide
+c trifad(2,0) : triangle central de ce decoupage
+c trifad(2,1) : triangle du cote de la face F1 du pentaedre
+c trifad(2,2) : triangle du cote de la face F2 du pentaedre
+c cotrvo(2,0/1/2) : futur code du triangle trifad(p,0/1/2) dans la
+c description des fils
+c areqtr(2,1) : arete interne au quadrangle de bord et bordant le
+c triangle trifad(p,1)
+c areqtr(2,2) : arete interne au quadrangle de bord et bordant le
+c triangle trifad(p,2)
+c====
+c
+ iaux = facpen(lepent,nulofa(2))
+ jaux = cofape(lepent,nulofa(2))
+ trifad(2,0) = -filqua(iaux)
+ if ( jaux.lt.5 ) then
+ cotrvo(2,0) = 1
+ trifad(2,1) = trifad(2,0) + 1
+ cotrvo(2,1) = 4
+ trifad(2,2) = trifad(2,0) + 2
+ cotrvo(2,2) = 1
+ areqtr(2,1) = aretri(trifad(2,0),1)
+ areqtr(2,2) = aretri(trifad(2,0),3)
+ else
+ cotrvo(2,0) = 5
+ trifad(2,1) = trifad(2,0) + 2
+ cotrvo(2,1) = 1
+ trifad(2,2) = trifad(2,0) + 1
+ cotrvo(2,2) = 4
+ areqtr(2,1) = aretri(trifad(2,0),3)
+ areqtr(2,2) = aretri(trifad(2,0),1)
+ endif
+#ifdef _DEBUG_HOMARD_
+ write(ulsort,90006) 'Quadrangle = ', iaux,', code =', jaux
+ do 3333 , iaux = 0, 2
+ write (ulsort,90015) 'trifad(2,0/1/2) =', trifad(2,iaux),
+ > ', aretes', (aretri(trifad(2,iaux),jaux),jaux=1,3)
+ 3333 continue
+ write(ulsort,90006) 'cotrvo(1,0) = ', cotrvo(1,0),
+ > 'cotrvo(1,1) = ', cotrvo(1,1),
+ > 'cotrvo(1,2) = ', cotrvo(1,2)
+ write(ulsort,90002) 'areqtr(2,1) = ', areqtr(1,1)
+ write(ulsort,90002) 'areqtr(2,2) = ', areqtr(1,2)
+#endif
+c
+c====
+c 4. niveau des triangles des conformites des faces
+c====
+c
+ niveau = nivtri(trifad(1,0))
+#ifdef _DEBUG_HOMARD_
+ write(ulsort,90002) 'niveau', niveau
+#endif
+c
+c====
+c 5. 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