MED fichier: f/2.3.6/test25.f

MED fichier

- Documentation MED
- Guides d'utilisation
- Guides de référence

f/2.3.6/test25.f
1 C* This file is part of MED.
2 C*
3 C* COPYRIGHT (C) 1999 - 2019 EDF R&D, CEA/DEN
4 C* MED is free software: you can redistribute it and/or modify
5 C* it under the terms of the GNU Lesser General Public License as published by
6 C* the Free Software Foundation, either version 3 of the License, or
7 C* (at your option) any later version.
8 C*
9 C* MED is distributed in the hope that it will be useful,
10 C* but WITHOUT ANY WARRANTY; without even the implied warranty of
11 C* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 C* GNU Lesser General Public License for more details.
13 C*
14 C* You should have received a copy of the GNU Lesser General Public License
15 C* along with MED. If not, see <http://www.gnu.org/licenses/>.
16 C*
17
18 C *******************************************************************************
19 C * - Nom du fichier : test25.f
20 C *
21 C * - Description : ecriture de mailles MED_POLYEDRE dans un maillage MED
22 C *
23 C ******************************************************************************
24 program test25
25 C
26 implicit none
27 include 'med.hf'
28 C
29 integer*8 fid
30 integer cret,mdim
31 parameter(mdim = 3)
32 character*32 maa
33 integer n
34 parameter(n=2)
35 C Connectivite nodale
36 integer np,nf
37 parameter(nf=9,np=3)
38 integer indexp(np),indexf(nf)
39 integer conn(24)
40 C Connectivite descendante
41 integer np2,nf2
42 parameter(nf2=8,np2=3)
43 integer indexp2(np2),indexf2(nf2)
44 integer conn2(nf2)
45 character*16 nom(n)
46 integer num(n),fam(n)
47 C
48 data indexp / 1,5,9 /
49 data indexf / 1,4,7,10,13,16,19,22,25 /
50 data conn / 1,2,3,4,5,6,7,8,9,10,11,12,13,14,
51 & 15,16,17,18,19,20,21,22,23,24 /
52 data indexp2 / 1,5,9 /
53 data indexf2 / med_tria3,med_tria3,med_tria3,med_tria3,
54 & med_tria3,med_tria3,med_tria3,med_tria3 /
55 data conn2 / 1,2,3,4,5,6,7,8 /
56 data nom / "poly1", "poly2"/
57 data num / 1,2 /, fam / 0,-1 /
58 data maa /"maa1"/
59
60 C ** Creation du fichier test25.med **
61 call efouvr(fid,'test25.med',med_lecture_ecriture, cret)
62 print *,cret
63 if (cret .ne. 0 ) then
64 print *,'Erreur creation du fichier'
65 call efexit(-1)
66 endif
67 print *,'Creation du fichier test25.med'
68
69 C ** Creation du maillage **
70 call efmaac(fid,maa,mdim,med_non_structure,
71 & 'un maillage pour test25',cret)
72 if (cret .ne. 0 ) then
73 print *,'Erreur creation du maillage'
74 call efexit(-1)
75 endif
76 print *,cret
77 print *,'Creation du maillage'
78
79 C ** Ecriture des connectivites des mailles polyedres en mode nodal **
80 call efpece(fid,maa,indexp,np,indexf,nf,conn,med_nod,cret)
81 print *,cret
82 if (cret .ne. 0 ) then
83 print *,'Erreur ecriture connectivite des polyedres'
84 call efexit(-1)
85 endif
86 print *,'Ecriture des connectivites des mailles
87 & de type MED_POLYEDRE'
88 print *,'Description nodale'
89
90 C ** Ecriture des connectivites des mailles polyedres en mode descendant **
91 call efpece(fid,maa,indexp2,np2,indexf2,nf2,conn2,med_desc,cret)
92 print *,cret
93 if (cret .ne. 0 ) then
94 print *,'Erreur ecriture connectivite des polyedres'
95 call efexit(-1)
96 endif
97 print *,'Ecriture des connectivites des mailles
98 & de type MED_POLYEDRE'
99 print *,'Description descendante'
100
101 C ** Ecriture des noms des mailles polyedres **
102 call efnome(fid,maa,nom,n,med_maille,med_polyedre,
103 & cret)
104 print *,cret
105 if (cret .ne. 0 ) then
106 print *,'Erreur ecriture noms des polyedres'
107 call efexit(-1)
108 endif
109 print *,'Ecriture des noms des polyedress'
110
111 C ** Ecriture des numeros des mailles polyedres **
112 call efnume(fid,maa,num,n,med_maille,med_polyedre,
113 & cret)
114 print *,cret
115 if (cret .ne. 0 ) then
116 print *,'Erreur ecriture numeros des polyedres'
117 call efexit(-1)
118 endif
119 print *,'Ecriture des numeros des polyedres'
120
121 C ** Ecriture des numeros des familles des segments **
122 call effame(fid,maa,fam,n,
123 & med_maille,med_polyedre,cret)
124 print *,cret
125 if (cret .ne. 0 ) then
126 print *,'Erreur ecriture numeros de familles polyedres'
127 call efexit(-1)
128 endif
129 print *,'Ecriture des numeros de familles des polyedres'
130
131 C ** Fermeture du fichier **
132 call efferm (fid,cret)
133 print *,cret
134 if (cret .ne. 0 ) then
135 print *,'Erreur fermeture du fichier'
136 call efexit(-1)
137 endif
138 print *,'Fermeture du fichier'
139 C
140 end