MED fichier
f/2.3.6/test25.f
1 C* This file is part of MED.
2 C*
3 C* COPYRIGHT (C) 1999 - 2017 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 cret, fid,mdim
30  parameter(mdim = 3)
31  character*32 maa
32  integer n
33  parameter(n=2)
34 C Connectivite nodale
35  integer np,nf
36  parameter(nf=9,np=3)
37  integer indexp(np),indexf(nf)
38  integer conn(24)
39 C Connectivite descendante
40  integer np2,nf2
41  parameter(nf2=8,np2=3)
42  integer indexp2(np2),indexf2(nf2)
43  integer conn2(nf2)
44  character*16 nom(n)
45  integer num(n),fam(n)
46 C
47  data indexp / 1,5,9 /
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 /
57  data maa /"maa1"/
58 
59 C ** Creation du fichier test25.med **
60  call efouvr(fid,'test25.med',med_lecture_ecriture, cret)
61  print *,cret
62  if (cret .ne. 0 ) then
63  print *,'Erreur creation du fichier'
64  call efexit(-1)
65  endif
66  print *,'Creation du fichier test25.med'
67 
68 C ** Creation du maillage **
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'
73  call efexit(-1)
74  endif
75  print *,cret
76  print *,'Creation du maillage'
77 
78 C ** Ecriture des connectivites des mailles polyedres en mode nodal **
79  call efpece(fid,maa,indexp,np,indexf,nf,conn,med_nod,cret)
80  print *,cret
81  if (cret .ne. 0 ) then
82  print *,'Erreur ecriture connectivite des polyedres'
83  call efexit(-1)
84  endif
85  print *,'Ecriture des connectivites des mailles
86  & de type MED_POLYEDRE'
87  print *,'Description nodale'
88 
89 C ** Ecriture des connectivites des mailles polyedres en mode descendant **
90  call efpece(fid,maa,indexp2,np2,indexf2,nf2,conn2,med_desc,cret)
91  print *,cret
92  if (cret .ne. 0 ) then
93  print *,'Erreur ecriture connectivite des polyedres'
94  call efexit(-1)
95  endif
96  print *,'Ecriture des connectivites des mailles
97  & de type MED_POLYEDRE'
98  print *,'Description descendante'
99 
100 C ** Ecriture des noms des mailles polyedres **
101  call efnome(fid,maa,nom,n,med_maille,med_polyedre,
102  & cret)
103  print *,cret
104  if (cret .ne. 0 ) then
105  print *,'Erreur ecriture noms des polyedres'
106  call efexit(-1)
107  endif
108  print *,'Ecriture des noms des polyedress'
109 
110 C ** Ecriture des numeros des mailles polyedres **
111  call efnume(fid,maa,num,n,med_maille,med_polyedre,
112  & cret)
113  print *,cret
114  if (cret .ne. 0 ) then
115  print *,'Erreur ecriture numeros des polyedres'
116  call efexit(-1)
117  endif
118  print *,'Ecriture des numeros des polyedres'
119 
120 C ** Ecriture des numeros des familles des segments **
121  call effame(fid,maa,fam,n,
122  & med_maille,med_polyedre,cret)
123  print *,cret
124  if (cret .ne. 0 ) then
125  print *,'Erreur ecriture numeros de familles polyedres'
126  call efexit(-1)
127  endif
128  print *,'Ecriture des numeros de familles des polyedres'
129 
130 C ** Fermeture du fichier **
131  call efferm (fid,cret)
132  print *,cret
133  if (cret .ne. 0 ) then
134  print *,'Erreur fermeture du fichier'
135  call efexit(-1)
136  endif
137  print *,'Fermeture du fichier'
138 C
139  end