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'