1 subroutine utoraq ( somare, a1, a2, a3, a4,
3 c ______________________________________________________________________
7 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
9 c Version originale enregistree le 18 juin 1996 sous le numero 96036
10 c aupres des huissiers de justice Simart et Lavoir a Clamart
11 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
12 c aupres des huissiers de justice
13 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
15 c HOMARD est une marque deposee d'Electricite de France
21 c ______________________________________________________________________
23 c UTilitaire - ORientation des Aretes d'un Quadrangle
25 c Remarque : cela suppose que les aretes a1, a2, a3, a4 sont donnees
26 c dans l'ordre standard
27 c Remarque : programme semblable a utsoqu
28 c ______________________________________________________________________
30 c . nom . e/s . taille . description .
31 c .____________________________________________________________________.
32 c . somare . e .2*nbaret. numeros des extremites d'arete .
33 c .a1,..,a4. e . 1 . les numeros d'aretes du quadrangle .
34 c . ori . s . 1 . 1 ou -1 , selon que l'arete ari est dans .
35 c . . . . le sens du quadrangle .
36 c ______________________________________________________________________
39 c 0. declarations et dimensionnement
42 c 0.1. ==> generalites
52 integer a1, a2, a3, a4
53 integer or1, or2, or3, or4
55 c 0.4. ==> variables locales
58 c 1. recherche des orientations :
59 c l'orientation de l'arete i est positive si son second somemt est
60 c un des sommets de l'arete i+1
73 if ( somare(2,a1) .eq. somare(1,a2) .or.
74 > somare(2,a1) .eq. somare(2,a2) ) then
75 c le 2nd noeud de l'arete 1 est un sommet de a2 ;
76 c donc l'orientation est positive
86 if ( somare(2,a2) .eq. somare(1,a3) .or.
87 > somare(2,a2) .eq. somare(2,a3) ) then
88 c le 2nd noeud de l'arete 2 est un sommet de a3 ;
89 c donc l'orientation est positive
99 if ( somare(2,a3) .eq. somare(1,a4) .or.
100 > somare(2,a3) .eq. somare(2,a4) ) then
101 c le 2nd noeud de l'arete 3 est un sommet de a4 ;
102 c donc l'orientation est positive
112 if ( somare(2,a4) .eq. somare(1,a1) .or.
113 > somare(2,a4) .eq. somare(2,a1) ) then
114 c le 2nd noeud de l'arete 4 est un sommet de a1 ;
115 c donc l'orientation est positive