MED fichier: f/2.3.6/test28.f

MED fichier

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

f/2.3.6/test28.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 : test28.f
20 C *
21 C * - Description : lecture des maillages structures (grille cartesienne |
22 C * grille de-structuree ) dans le fichier test27.med
23 C *
24 C *****************************************************************************
25 program test28
26 C
27 implicit none
28 include 'med.hf'
29 C
30 C
31 integer*8 fid
32 integer cret,i,j
33 C ** la dimension du maillage **
34 integer mdim,nind,nmaa,type,quoi,rep,typmaa
35 C ** nom du maillage de longueur maxi MED_TAILLE_NOM **
36 character*32 maa
37 C ** le nombre de noeuds **
38 integer nnoe
39 C ** table des coordonnees **
40 real*8 coo(8)
41 character*16 comp, comp2(2)
42 character*16 unit, unit2(2)
43 character*200 desc
44 integer strgri(2)
45 C ** grille cartesienne **
46 integer axe
47 real*8 indice(4)
48 integer tmp
49
50 C
51 C On ouvre le fichier test27.med en lecture seule
52 call efouvr(fid,'test27.med',med_lecture, cret)
53 if (cret .ne. 0 ) then
54 print *,'Erreur ouverture du fichier'
55 call efexit(-1)
56 endif
57 print *,cret
58
59 print *,'Ouverture du fichier test27.med'
60 C
61 C Combien de maillage ?
62 call efnmaa(fid,nmaa,cret)
63 print *,cret
64 if (cret .ne. 0 ) then
65 print *,'Erreur lecture du nombre de maillage'
66 call efexit(-1)
67 endif
68 C
69 C On boucle sur les maillages et on ne lit que les
70 C maillages structures
71 do 10 i=1,nmaa
72 C
73 C On repere les maillages qui nous interessent
74 C
75 call efmaai(fid,i,maa,mdim,typmaa,desc,cret)
76 print *,cret
77 if (cret .ne. 0 ) then
78 print *,'Erreur lecture maillage info'
79 call efexit(-1)
80 endif
81 print *,'Maillge de nom : ',maa
82 print *,'- Dimension : ',mdim
83 if (typmaa.eq.med_structure) then
84 print *,'- Type : MED_STRUCTURE'
85 else
86 print *,'- Type : MED_NON_STRUCTURE'
87 endif
88 C
89 C On repere le type de la grille
90 if (typmaa.eq.med_structure) then
91 call efnagl(fid,maa,type,cret)
92 print *,cret
93 if (cret .ne. 0 ) then
94 print *,'Erreur lecture nature de la grille'
95 call efexit(-1)
96 endif
97 if (type.eq.med_grille_cartesienne) then
98 print *,'- Nature de la grille :',
99 & 'MED_GRILLE_CARTESIENNE'
100 endif
101 if (type.eq.med_grille_standard) then
102 print *,'- Nature de la grille : MED_GRILLE_STANDARD'
103 endif
104 endif
105 C
106 C On regarde la structure et les coordonnees de la grille MED_GRILLE_STANDARD
107 if ((type.eq.med_grille_standard)
108 & .and. (typmaa.eq.med_structure)) then
109 C
110 call efnema(fid,maa,med_coor,med_noeud,0,0,nnoe,cret)
111 print *,cret
112 if (cret .ne. 0 ) then
113 print *,'Erreur lecture nombre de noeud'
114 call efexit(-1)
115 endif
116 print *,'- Nombre de noeuds : ',nnoe
117 C
118 call efscol(fid,maa,mdim,strgri,cret)
119 print *,cret
120 if (cret .ne. 0 ) then
121 print *,'Erreur lecture structure de la grille'
122 call efexit(-1)
123 endif
124 print *,'- Structure de la grille : ',strgri
125 C
126 call efcool(fid,maa,mdim,coo,
127 & med_full_interlace,med_all,tmp,
128 & 0,rep,comp2,unit2,cret)
129 print *,cret
130 if (cret .ne. 0 ) then
131 print *,'Erreur lecture des coordonnees des noeuds'
132 call efexit(-1)
133 endif
134 print *,'- Coordonnees :'
135 do 20 j=1,nnoe*mdim
136 print *,coo(j)
137 20 continue
138 endif
139 C
140 if ((type.eq.med_grille_cartesienne)
141 & .and. (typmaa.eq.med_structure)) then
142 C
143 do 30 axe=1,mdim
144 if (axe.eq.1) then
145 quoi = med_coor_ind1
146 endif
147 if (axe.eq.2) then
148 quoi = med_coor_ind2
149 endif
150 if (axe.eq.3) then
151 quoi = med_coor_ind3
152 endif
153 C Lecture de la taille de l'indice selon la dimension
154 C fournie par le parametre quoi
155 call efnema(fid,maa,quoi,med_noeud,0,0,nind,cret)
156 print *,cret
157 if (cret .ne. 0 ) then
158 print *,'Erreur lecture taille indice'
159 call efexit(-1)
160 endif
161 print *,'- Axe ',axe
162 print *,'- Nombre d indices : ',nind
163 C Lecture des indices des coordonnees de la grille
164 call eficol(fid,maa,mdim,indice,nind,axe,comp,unit,
165 & cret)
166 print *,cret
167 if (cret .ne. 0 ) then
168 print *,'Erreur lecture indices de coordonnées'
169 call efexit(-1)
170 endif
171 print *,'- Axe ',comp
172 print *,' unite : ',unit
173 do 40 j=1,nind
174 print *,indice(j)
175 40 continue
176 30 continue
177 C
178 endif
179 C
180 10 continue
181 C
182 C On ferme le fichier
183 call efferm (fid,cret)
184 print *,cret
185 if (cret .ne. 0 ) then
186 print *,'Erreur fermeture du fichier'
187 call efexit(-1)
188 endif
189 print *,'Fermeture du fichier'
190 C
191 end
192