30 integer cret,fid,mdim,nmaa,npoly,i,j,k,l
31 integer nfaces, nnoeuds
37 integer np,nf,np2,nf2,taille,tmp
38 parameter(np=3,nf=9,np2=3,nf2=8)
39 integer indexp(np),indexf(nf)
41 integer indexp2(np2),indexf2(nf2)
48 call efouvr(fid,
'test25.med',med_lecture, cret)
50 if (cret .ne. 0 )
then
51 print *,
'Erreur ouverture du fichier'
54 print *,
'Ouverture du fichier test25.med'
57 call efnmaa(fid,nmaa,cret)
59 if (cret .ne. 0 )
then
60 print *,
'Erreur lecture du nombre de maillage'
63 print *,
'Nombre de maillages : ',nmaa
70 call efmaai(fid,i,maa,mdim,
type,desc,cret)
72 if (cret .ne. 0 )
then
73 print *,
'Erreur infos maillage'
76 print *,
'Maillage : ',maa
77 print *,
'Dimension : ',mdim
80 call efnema(fid,maa,med_conn,med_maille,med_polyedre,
83 if (cret .ne. 0 )
then
84 print *,
'Erreur lecture nombre de polyedre'
87 print *,
'Nombre de mailles MED_POLYEDRE : ',npoly
90 call efpyei(fid,maa,med_nod,tmp,taille,cret)
92 if (cret .ne. 0 )
then
93 print *,
'Erreur infos sur les polyedres'
96 print *,
'Taille de la connectivite : ',taille
97 print *,
'Taille du tableau indexf : ',tmp
100 call efpecl(fid,maa,indexp,npoly+1,indexf,tmp,conn,
103 if (cret .ne. 0 )
then
104 print *,
'Erreur lecture connectivites polyedres'
107 print *,
'Lecture de la connectivite des polyedres'
108 print *,
'Connectivite nodale'
111 call efpecl(fid,maa,indexp2,npoly+1,indexf2,tmp,conn2,
114 if (cret .ne. 0 )
then
115 print *,
'Erreur lecture connectivite des polyedres'
118 print *,
'Lecture de la connectivite des polyedres'
119 print *,
'Connectivite descendante'
122 call efnoml(fid,maa,nom,npoly,med_maille,med_polyedre,
125 if (cret .ne. 0 )
then
126 print *,
'Erreur lecture noms des polyedres'
129 print *,
'Lecture des noms'
132 call efnuml(fid,maa,num,npoly,med_maille,med_polyedre,
135 if (cret .ne. 0 )
then
136 print *,
'Erreur lecture des numeros des polyedres'
139 print *,
'Lecture des numeros'
142 call effaml(fid,maa,fam,npoly,med_maille,med_polyedre,
145 if (cret .ne. 0 )
then
146 print *,
'Erreur lecture numeros de famille polyedres'
149 print *,
'Lecture des numeros de famille'
152 print *,
'Affichage des resultats'
155 print *,
'>> Maille polygone ',j
156 print *,
'---- Connectivite nodale ---- : '
157 nfaces = indexp(j+1) - indexp(j)
163 ind2 = indexf(ind1+k-1)
164 nnoeuds = indexf(ind1+k) - indexf(ind1+k-1)
167 print *,
' ',conn(ind2+l-1)
170 print *,
'---- Connectivite descendante ---- : '
171 nfaces = indexp2(j+1) - indexp2(j)
176 print *,
' => Numero : ',conn2(ind1+k-1)
177 print *,
' => Type : ',indexf2(ind1+k-1)
179 print *,
'---- Nom ---- : ',nom(j)
180 print *,
'---- Numero ----: ',num(j)
181 print *,
'---- Numero de famille ---- : ',fam(j)
188 call efferm (fid,cret)
190 if (cret .ne. 0 )
then
191 print *,
'Erreur fermeture du fichier'
194 print *,
'Fermeture du fichier'