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