1 subroutine pcsi00 ( nbfop2, profho, vap2ho,
4 c ______________________________________________________________________
8 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
10 c Version originale enregistree le 18 juin 1996 sous le numero 96036
11 c aupres des huissiers de justice Simart et Lavoir a Clamart
12 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
13 c aupres des huissiers de justice
14 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
16 c HOMARD est une marque deposee d'Electricite de France
22 c ______________________________________________________________________
24 c aPres adaptation - Conversion de Solution -
26 c interpolation iso/p2 sur les noeuds - phase 00
28 c Moyenne des valeurs aux extremites d'un ensmble d'aretes
29 c ______________________________________________________________________
31 c . nom . e/s . taille . description .
32 c .____________________________________________________________________.
33 c . nbfop2 . e . 1 . nombre de fonctions P2 .
34 c . profho . es . * . pour chaque entite en numerotation homard :.
35 c . . . . 0 : l'entite est absente du profil .
36 c . . . . 1 : l'entite est presente dans le profil .
37 c . vap2ho . es . nbfop2*. variables p2 numerotation homard .
39 c . somare . e .2*nbarto. numeros des extremites d'arete .
40 c . np2are . e . nbarto . numero des noeuds p2 milieux d'aretes .
41 c . nbarco . e . 1 . nombre d'aretes concernees .
42 c . nuaret . e . nbarco . numero des aretes a traiter .
43 c ______________________________________________________________________
46 c 0. declarations et dimensionnement
49 c 0.1. ==> generalites
64 integer somare(2,nbarto), np2are(nbarto)
65 integer nbarco, nuaret(nbarco)
67 double precision vap2ho(nbfop2,*)
69 c 0.4. ==> variables locales
73 integer s1, s2, sm, nuv
74 c ______________________________________________________________________
78 cgn write(1,90002) 'Aretes', (nuaret(iaux),iaux=1,nbarco)
80 do 10 , iaux = 1 , nbarco
83 c 1. L'arete concernee
86 cgn write(1,90002) 'larete',larete
95 cgn write(1,90002) 'sm =',sm
99 do 21 , nuv = 1, nbfop2
101 cgn write(1,90002) 'sommets',s1, s2)
102 vap2ho(nuv,sm) = unsde * ( vap2ho(nuv,s1) + vap2ho(nuv,s2) )
103 cgn write(1,*) 'vap2ho(nuv,',sm,') =',vap2ho(nuv,sm)