MED fichier
f/2.3.6/test27.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 : test27.f
20 C *
21 C * - Description : creation de maillages structures (grille cartesienne |
22 C * grille standard ) dans le fichier test27.med
23 C *
24 C *****************************************************************************
25  program test27
26 C
27  implicit none
28  include 'med.hf'
29 C
30 C
31  integer cret, fid
32 C ** la dimension du maillage **
33  integer mdim
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,nind
46  real*8 indice(4)
47 
48 C
49 C
50  data coo /0.0,0.0,1.0,0.0,0.0,1.0,1.0,1.0/
51  data comp2 /"x","y"/, unit2 /"cm","cm"/
52 C
53 C Creation du fichier test27.med
54  call efouvr(fid,'test27.med',med_lecture_ecriture, cret)
55  print *,cret
56  if (cret .ne. 0 ) then
57  print *,'Erreur creation du fichier'
58  call efexit(-1)
59  endif
60  print *,'Creation du fichier test27.med'
61 C
62 C Creation d'un maillage MED_NON_STRUCTURE
63  mdim = 3
64  maa = 'maillage vide'
65  desc = 'un maillage vide'
66  call efmaac(fid,maa,mdim,med_non_structure,desc,cret)
67  print *,cret
68  if (cret .ne. 0 ) then
69  print *,'Erreur creation du maillage'
70  call efexit(-1)
71  endif
72 C
73 C Creation d'une grille cartesienne
74  mdim = 2
75  maa = 'grille cartesienne'
76  desc = 'un exemple de grille cartesienne'
77  call efmaac(fid,maa,mdim,med_structure,desc,cret)
78  print *,cret
79  if (cret .ne. 0 ) then
80  print *,'Erreur creation du maillage'
81  call efexit(-1)
82  endif
83  print *,'Creation d un maillage MED_STRUCTURE'
84 
85 C
86 C On specifie la nature du maillage structure
87  call efnage(fid,maa,med_grille_cartesienne,cret)
88  print *,cret
89  print *,'On definit la nature de la grille :
90  & MED_GRILLE_CARTESIENNE'
91  if (cret .ne. 0 ) then
92  print *,'Erreur ecriture de la nature de la grille'
93  call efexit(-1)
94  endif
95 C
96 C On definit les indices de la grille selon chaque dimension
97  indice(1) = 1.1d0
98  indice(2) = 1.2d0
99  indice(3) = 1.3d0
100  indice(4) = 1.4d0
101  nind = 4
102  axe = 1
103  comp = 'X'
104  unit = 'cm'
105  call eficoe(fid,maa,mdim,indice,nind,axe,comp,unit,cret)
106  print *,cret
107  if (cret .ne. 0 ) then
108  print *,'Erreur ecriture des indices'
109  call efexit(-1)
110  endif
111  print *,'Ecriture des indices des coordonnees selon axe X'
112 C
113  indice(1) = 2.1d0
114  indice(2) = 2.2d0
115  indice(3) = 2.3d0
116  indice(4) = 2.4d0
117  nind = 4
118  axe = 2
119  comp = 'Y'
120  unit = 'cm'
121  call eficoe(fid,maa,mdim,indice,nind,axe,comp,unit,cret)
122  print *,cret
123  if (cret .ne. 0 ) then
124  print *,'Erreur ecriture des indices'
125  call efexit(-1)
126  endif
127  print *,'Ecriture des indices des coordonnees selon axe Y'
128 C
129 C Creation d'une grille MED_GRILLE_STANDARD de dimension 2
130  maa = 'grille standard'
131  mdim = 2
132  desc = 'un exemple de grille standard'
133  call efmaac(fid,maa,mdim,med_structure,desc,cret)
134  print *,cret
135  if (cret .ne. 0 ) then
136  print *,'Erreur creation de maillage'
137  call efexit(-1)
138  endif
139  print *,'Nouveau maillage MED_STRUCTURE'
140 C
141  call efnage(fid,maa,med_grille_standard,cret)
142  print *,cret
143  if (cret .ne. 0 ) then
144  print *,'Erreur ecriture de la nature de la grille'
145  call efexit(-1)
146  endif
147  print *,'On definit la nature du maillage : MED_GRILLE_STANDARD'
148 C
149 C On ecrit les coordonnes de la grille
150  nnoe = 4
151  call efcooe(fid,maa,mdim,coo,med_full_interlace,nnoe,med_cart,
152  & comp2,unit2,cret)
153  print *,cret
154  if (cret .ne. 0 ) then
155  print *,'Erreur ecriture des coordonnees des noeuds'
156  call efexit(-1)
157  endif
158  print *,'Ecriture des coordonnees de la grille'
159 C
160 C On definit la structure des coordonnees de la grille
161  strgri(1) = 2
162  strgri(2) = 2
163  call efscoe(fid,maa,mdim,strgri,cret)
164  print *,cret
165  if (cret .ne. 0 ) then
166  print *,'Erreur ecriture de la structure'
167  call efexit(-1)
168  endif
169  print *,'Ecriture de la structure de la grille : / 2,2 /'
170 C
171 C On ferme le fichier
172  call efferm (fid,cret)
173  print *,cret
174  if (cret .ne. 0 ) then
175  print *,'Erreur fermeture du fichier'
176  call efexit(-1)
177  endif
178  print *,'Fermeture du fichier'
179 C
180  end
181 
182 
183 
184 
185 
186