MED fichier: f/2.3.6/test24.f

MED fichier

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

f/2.3.6/test24.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 : test24.f
20 C *
21 C * - Description : lecture de mailles MED_POLYGONE dans le maillage MED
22 C * du fichier test23.med
23 C *
24 C ******************************************************************************
25 program test23
26 C
27 implicit none
28 include 'med.hf'
29 C
30 integer*8 fid
31 integer cret,mdim,nmaa,npoly,i,j,k,taille
32 character*32 maa
33 character*200 desc
34 integer ni, n
35 parameter(ni=4, n=3)
36 integer index(ni),ind1,ind2
37 character*16 nom(n)
38 integer num(n),fam(n)
39 integer con(16)
40 integer type
41 C
42 C Ouverture du fichier test23.med en lecture seule
43 call efouvr(fid,'test23.med',med_lecture, cret)
44 print *,cret
45 if (cret .ne. 0 ) then
46 print *,'Erreur ouverture du fichier'
47 call efexit(-1)
48 endif
49 print *,'Ouverture du fichier test23.med'
50 C
51 C Lecture du nombre de maillages
52 call efnmaa(fid,nmaa,cret)
53 print *,cret
54 if (cret .ne. 0 ) then
55 print *,'Erreur lecture nombre de maillage'
56 call efexit(-1)
57 endif
58 print *,'Nombre de maillages : ',nmaa
59 C
60 C Lecture de toutes les mailles MED_POLYGONE
61 C dans chaque maillage
62 do 10 i=1,nmaa
63 C
64 C Info sur chaque maillage
65 call efmaai(fid,i,maa,mdim,type,desc,cret)
66 if (cret .ne. 0 ) then
67 print *,'Erreur lecture infos maillage'
68 call efexit(-1)
69 endif
70 print *,cret
71 print *,'Maillage : ',maa
72 print *,'Dimension : ',mdim
73 C
74 C Combien de mailles polygones
75 call efnema(fid,maa,med_conn,med_maille,med_polygone,
76 & med_nod,npoly,cret)
77 print *,cret
78 if (cret .ne. 0 ) then
79 print *,'Erreur lecture du nombre de polygone'
80 call efexit(-1)
81 endif
82 print *,'Nombre de mailles MED_POLYGONE : ',npoly
83 C
84 C Taille des connectivites
85 call efpygi(fid,maa,med_maille,med_nod,taille,cret)
86 print *,cret
87 if (cret .ne. 0 ) then
88 print *,'Erreur lecture infos polygones'
89 call efexit(-1)
90 endif
91 print *,'Taille de la connectivite : ',taille
92 C
93 C Lecture de la connectivite
94 call efpgcl(fid,maa,index,npoly+1,con,med_maille,
95 & med_nod,cret)
96 print *,cret
97 if (cret .ne. 0 ) then
98 print *,'Erreur lecture des connectivites polygones'
99 call efexit(-1)
100 endif
101 print *,'Lecture de la connectivite des polygones'
102 C
103 C Lecture des noms
104 call efnoml(fid,maa,nom,npoly,med_maille,med_polygone,
105 & cret)
106 print *,cret
107 if (cret .ne. 0 ) then
108 print *,'Erreur lecture des noms des polygones'
109 call efexit(-1)
110 endif
111 print *,'Lecture des noms'
112 C
113 C Lecture des numeros
114 call efnuml(fid,maa,num,npoly,med_maille,med_polygone,
115 & cret)
116 print *,cret
117 if (cret .ne. 0 ) then
118 print *,'Erreur lecture des numeros des polygones'
119 call efexit(-1)
120 endif
121 print *,'Lecture des numeros'
122 C
123 C Lecture des numeros de familles
124 call effaml(fid,maa,fam,npoly,med_maille,med_polygone,
125 & cret)
126 print *,cret
127 if (cret .ne. 0 ) then
128 print *,'Erreur lecture des numeros de famille des
129 & polygones'
130 call efexit(-1)
131 endif
132 print *,'Lecture des numeros de famille'
133 C
134 C Affichage des resultats
135 print *,'Affichage des resultats'
136 do 20 j=1,npoly
137 C
138 print *,'>> Maille polygone ',j
139 print *,'---- Connectivite ---- : '
140 ind1 = index(j)
141 ind2 = index(j+1)
142 do 30 k=ind1,ind2-1
143 print *,con(k)
144 30 continue
145 print *,'---- Nom ---- : ',nom(j)
146 print *,'---- Numero ----: ',num(j)
147 print *,'---- Numero de famille ---- : ',fam(j)
148 C
149 20 continue
150 C
151 10 continue
152 C
153 C Fermeture du fichier
154 call efferm (fid,cret)
155 print *,cret
156 if (cret .ne. 0 ) then
157 print *,'Erreur fermeture du fichier'
158 call efexit(-1)
159 endif
160 print *,'Fermeture du fichier'
161 C
162 end