1 subroutine pcsiar ( nbfop2, profho, vap2ho,
2 > hetare, somare, np2are, filare )
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 aPres adaptation - Conversion de Solution -
25 c interpolation iso-p2 sur les noeuds lors du decoupage des ARetes
27 c remarque : pcs2ar et pcsiar sont des clones
28 c ______________________________________________________________________
30 c . nom . e/s . taille . description .
31 c .____________________________________________________________________.
32 c . nbfop2 . e . 1 . nombre de fonctions P2 .
33 c . profho . es . * . pour chaque entite en numerotation homard :.
34 c . . . . 0 : l'entite est absente du profil .
35 c . . . . 1 : l'entite est presente dans le profil .
36 c . vap2ho . es . nbfop2*. variables iso-p2 numerotation homard .
38 c . somare . e .2*nbarto. numeros des extremites d'arete .
39 c . np2are . e . nbarto . numero des noeuds p2 milieux d'aretes .
40 c . filare . e . nbarto . premiere fille des aretes .
41 c . hetare . e . nbarto . historique de l'etat des aretes .
42 c ______________________________________________________________________
45 c 0. declarations et dimensionnement
48 c 0.1. ==> generalites
63 integer profho(nbnoto)
64 integer hetare(nbarto), somare(2,nbarto), np2are(nbarto)
65 integer filare(nbarto)
67 double precision vap2ho(nbfop2,*)
69 c 0.4. ==> variables locales
71 integer larete, a1, a2, sm, s1, s2, m1, m2, nuv
72 c ______________________________________________________________________
75 c 1. interpolation pour les aretes qui viennent d'etre decoupees
78 do 10 , larete = 1, nbarto
80 if ( hetare(larete).eq.2 ) then
82 c recuperation des aretes filles
87 c recuperation des sommets de l'arete
92 c recuperation du nouveau noeud sommet
96 if ( profho(s1).eq.1 .and. profho(s2).eq.1 .and.
97 > profho(sm).eq.1 ) then
99 c recuperation des nouveaux noeuds milieux
106 c interpolation p1 : interpolee (ui,i=1,2) = 1/2 (u1+u2)
108 do 11, nuv = 1, nbfop2
111 > unsde * ( vap2ho(nuv,s1) + vap2ho(nuv,sm) )
113 > unsde * ( vap2ho(nuv,s2) + vap2ho(nuv,sm) )