1 subroutine deiarr ( nivmax,
8 > ulsort, langue, codret)
9 c ______________________________________________________________________
13 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
15 c Version originale enregistree le 18 juin 1996 sous le numero 96036
16 c aupres des huissiers de justice Simart et Lavoir a Clamart
17 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
18 c aupres des huissiers de justice
19 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
21 c HOMARD est une marque deposee d'Electricite de France
27 c ______________________________________________________________________
29 c traitement des DEcisions - Initialisation de l'indicateur entier
31 c - cas des ARetes - Raffinement
33 c ______________________________________________________________________
35 c . nom . e/s . taille . description .
36 c .____________________________________________________________________.
37 c . nivmax . e . 1 . niveau max a ne pas depasser en raffinement.
38 c . decare . s .0:nbarto. decisions des aretes .
39 c . hetare . e . nbarto . historique de l'etat des aretes .
40 c . posifa . e . nbarto . pointeur sur tableau facare .
41 c . facare . e . nbfaar . liste des faces contenant une arete .
42 c . nivtri . e . nbtrto . niveau des triangles .
43 c . nivqua . e . nbquto . niveau des quadrangles .
44 c . arsupp . e . nbarto . support pour les aretes .
45 c . arindi . e . nbarto . valeurs entieres pour les aretes .
46 c . ulsort . e . 1 . numero d'unite logique de la liste standard.
47 c . langue . e . 1 . langue des messages .
48 c . . . . 1 : francais, 2 : anglais .
49 c . codret . es . 1 . code de retour des modules .
50 c . . . . 0 : pas de probleme .
51 c . . . . 2 : probleme dans le traitement .
52 c ______________________________________________________________________
55 c 0. declarations et dimensionnement
58 c 0.1. ==> generalites
64 parameter ( nompro = 'DEIARR' )
79 integer decare(0:nbarto)
80 integer hetare(nbarto)
81 integer posifa(0:nbarto), facare(nbfaar)
82 integer nivtri(nbtrto)
83 integer nivqua(nbquto)
84 integer arsupp(nbarto), arindi(nbarto)
86 integer ulsort, langue, codret
88 c 0.4. ==> variables locales
92 integer iaux, jaux, kaux, ideb, ifin
95 parameter (nbmess = 10 )
96 character*80 texte(nblang,nbmess)
97 c ______________________________________________________________________
103 c 1.1. ==> Les messages
107 #ifdef _DEBUG_HOMARD_
108 write (ulsort,texte(langue,1)) 'Entree', nompro
117 c 2. traitement des indicateurs portant sur les aretes
118 c pour le filtrage sur les niveaux, on tient compte du fait que
119 c le niveau d'une arete est identifie a celui de l'une quelconque de
120 c ses faces voisines quand elle en a. Sinon, on ne filtre pas.
123 #ifdef _DEBUG_HOMARD_
124 write(ulsort,texte(langue,4)) mess14(langue,3,1)
127 #ifdef _DEBUG_HOMARD_
128 write(ulsort,texte(langue,5))
133 do 21 , larete = 1, nbarto
135 if ( arsupp(larete).ne.0 ) then
136 etat = mod(hetare(larete),10)
137 if ( etat.eq.0 ) then
138 if ( arindi(larete).eq.1 ) then
139 ideb = posifa(larete-1)+1
140 ifin = posifa(larete)
142 if ( ifin.ge.ideb .and. nivmax.ge.0 ) then
143 if ( facare(ideb).gt.0 ) then
144 kaux = nivtri(facare(ideb))
146 kaux = nivqua(-facare(ideb))
148 if ( kaux.ge.nivmax ) then
152 if ( jaux.eq.0 ) then
154 cgn write(ulsort,*) 'mise a 2 de decare pour arete', larete
164 if ( iaux.ne.0 ) then
165 write(ulsort,texte(langue,10))
166 write(ulsort,texte(langue,4)) mess14(langue,3,1)
167 write(ulsort,texte(langue,7)) nivmax
168 write(ulsort,texte(langue,9)) iaux
175 if ( codret.ne.0 ) then
179 write (ulsort,texte(langue,1)) 'Sortie', nompro
180 write (ulsort,texte(langue,2)) codret
184 #ifdef _DEBUG_HOMARD_
185 write (ulsort,texte(langue,1)) 'Sortie', nompro