31 integer cret,mdim,nmaa,npoly,i,j,k,l
32 integer nfaces, nnoeuds
38 integer np,nf,np2,nf2,taille,tmp
39 parameter(np=3,nf=9,np2=3,nf2=8)
40 integer indexp(np),indexf(nf)
42 integer indexp2(np2),indexf2(nf2)
49 call efouvr(fid,
'test25.med',med_lecture, cret)
51 if (cret .ne. 0 )
then
52 print *,
'Erreur ouverture du fichier'
55 print *,
'Ouverture du fichier test25.med'
58 call efnmaa(fid,nmaa,cret)
60 if (cret .ne. 0 )
then
61 print *,
'Erreur lecture du nombre de maillage'
64 print *,
'Nombre de maillages : ',nmaa
71 call efmaai(fid,i,maa,mdim,
type,desc,cret)
73 if (cret .ne. 0 )
then
74 print *,
'Erreur infos maillage'
77 print *,
'Maillage : ',maa
78 print *,
'Dimension : ',mdim
81 call efnema(fid,maa,med_conn,med_maille,med_polyedre,
84 if (cret .ne. 0 )
then
85 print *,
'Erreur lecture nombre de polyedre'
88 print *,
'Nombre de mailles MED_POLYEDRE : ',npoly
91 call efpyei(fid,maa,med_nod,tmp,taille,cret)
93 if (cret .ne. 0 )
then
94 print *,
'Erreur infos sur les polyedres'
97 print *,
'Taille de la connectivite : ',taille
98 print *,
'Taille du tableau indexf : ',tmp
101 call efpecl(fid,maa,indexp,npoly+1,indexf,tmp,conn,
104 if (cret .ne. 0 )
then
105 print *,
'Erreur lecture connectivites polyedres'
108 print *,
'Lecture de la connectivite des polyedres'
109 print *,
'Connectivite nodale'
112 call efpecl(fid,maa,indexp2,npoly+1,indexf2,tmp,conn2,
115 if (cret .ne. 0 )
then
116 print *,
'Erreur lecture connectivite des polyedres'
119 print *,
'Lecture de la connectivite des polyedres'
120 print *,
'Connectivite descendante'
123 call efnoml(fid,maa,nom,npoly,med_maille,med_polyedre,
126 if (cret .ne. 0 )
then
127 print *,
'Erreur lecture noms des polyedres'
130 print *,
'Lecture des noms'
133 call efnuml(fid,maa,num,npoly,med_maille,med_polyedre,
136 if (cret .ne. 0 )
then
137 print *,
'Erreur lecture des numeros des polyedres'
140 print *,
'Lecture des numeros'
143 call effaml(fid,maa,fam,npoly,med_maille,med_polyedre,
146 if (cret .ne. 0 )
then
147 print *,
'Erreur lecture numeros de famille polyedres'
150 print *,
'Lecture des numeros de famille'
153 print *,
'Affichage des resultats'
156 print *,
'>> Maille polygone ',j
157 print *,
'---- Connectivite nodale ---- : '
158 nfaces = indexp(j+1) - indexp(j)
164 ind2 = indexf(ind1+k-1)
165 nnoeuds = indexf(ind1+k) - indexf(ind1+k-1)
168 print *,
' ',conn(ind2+l-1)
171 print *,
'---- Connectivite descendante ---- : '
172 nfaces = indexp2(j+1) - indexp2(j)
177 print *,
' => Numero : ',conn2(ind1+k-1)
178 print *,
' => Type : ',indexf2(ind1+k-1)
180 print *,
'---- Nom ---- : ',nom(j)
181 print *,
'---- Numero ----: ',num(j)
182 print *,
'---- Numero de famille ---- : ',fam(j)
189 call efferm (fid,cret)
191 if (cret .ne. 0 )
then
192 print *,
'Erreur fermeture du fichier'
195 print *,
'Fermeture du fichier'