29 integer cret, fid,mdim
37 integer indexp(np),indexf(nf)
41 parameter(nf2=8,np2=3)
42 integer indexp2(np2),indexf2(nf2)
48 data indexf / 1,4,7,10,13,16,19,22,25 /
49 data conn / 1,2,3,4,5,6,7,8,9,10,11,12,13,14,
50 & 15,16,17,18,19,20,21,22,23,24 /
51 data indexp2 / 1,5,9 /
52 data indexf2 / med_tria3,med_tria3,med_tria3,med_tria3,
53 & med_tria3,med_tria3,med_tria3,med_tria3 /
54 data conn2 / 1,2,3,4,5,6,7,8 /
55 data nom /
"poly1",
"poly2"/
56 data num / 1,2 /, fam / 0,-1 /
60 call efouvr(fid,
'test25.med',med_lecture_ecriture, cret)
62 if (cret .ne. 0 )
then
63 print *,
'Erreur creation du fichier'
66 print *,
'Creation du fichier test25.med'
69 call efmaac(fid,maa,mdim,med_non_structure,
70 &
'un maillage pour test25',cret)
71 if (cret .ne. 0 )
then
72 print *,
'Erreur creation du maillage'
76 print *,
'Creation du maillage'
79 call efpece(fid,maa,indexp,np,indexf,nf,conn,med_nod,cret)
81 if (cret .ne. 0 )
then
82 print *,
'Erreur ecriture connectivite des polyedres'
85 print *,
'Ecriture des connectivites des mailles
86 & de type MED_POLYEDRE'
87 print *,
'Description nodale'
90 call efpece(fid,maa,indexp2,np2,indexf2,nf2,conn2,med_desc,cret)
92 if (cret .ne. 0 )
then
93 print *,
'Erreur ecriture connectivite des polyedres'
96 print *,
'Ecriture des connectivites des mailles
97 & de type MED_POLYEDRE'
98 print *,
'Description descendante'
101 call efnome(fid,maa,nom,n,med_maille,med_polyedre,
104 if (cret .ne. 0 )
then
105 print *,
'Erreur ecriture noms des polyedres'
108 print *,
'Ecriture des noms des polyedress'
111 call efnume(fid,maa,num,n,med_maille,med_polyedre,
114 if (cret .ne. 0 )
then
115 print *,
'Erreur ecriture numeros des polyedres'
118 print *,
'Ecriture des numeros des polyedres'
121 call effame(fid,maa,fam,n,
122 & med_maille,med_polyedre,cret)
124 if (cret .ne. 0 )
then
125 print *,
'Erreur ecriture numeros de familles polyedres'
128 print *,
'Ecriture des numeros de familles des polyedres'
131 call efferm (fid,cret)
133 if (cret .ne. 0 )
then
134 print *,
'Erreur fermeture du fichier'
137 print *,
'Fermeture du fichier'