1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
26
27 implicit none
28 include 'med.hf'
29
30
31 integer*8 fid
32 integer cret
33 character*64 maa
34 character*200 des
35 integer nmaa, mdim ,nnoe,type,sdim
36 character*16 nomcoo(2)
37 character*16 unicoo(2)
38 character(16) :: dtunit
39 integer nstep, stype, atype,chgt,tsf
40 integer numglb(100),i
41
42
43
44 call mfiope(fid,
'test31.med',med_acc_rdonly, cret)
45 print '(I1)',cret
46 if (cret .ne. 0 ) then
47 print *,'Erreur ouverture du fichier test31.med'
48 call efexit(-1)
49 endif
50
51
52
53 call mmhmii(fid,1,maa,sdim,mdim,
type,des,dtunit,
54 & stype,nstep,atype,nomcoo,unicoo,cret)
55 print '(I1)',cret
56 if (cret .ne. 0 ) then
57 print *,'Erreur acces au premier maillage'
58 call efexit(-1)
59 endif
60
61 nnoe = 0
62 call mmhnme(fid,maa,med_no_dt,med_no_it,med_node,med_none,
63 & med_coordinate,med_no_cmode,chgt,tsf,nnoe,cret)
64 if (cret .ne. 0 ) then
65 print *,'Erreur acces au nombre de noeud du premier maillage'
66 call efexit(-1)
67 endif
68
69
70 print '(A,I1,A,A4,A,I1,A,I4)','maillage '
71 & ,0,' de nom ',maa,' et de dimension ',mdim,
72 & ' comportant le nombre de noeud ',nnoe
73
74
75
76 call mmhgnr(fid,maa,med_no_dt,med_no_it,med_node,med_none,
77 & numglb,cret)
78
79 if (cret .ne. 0 ) then
80 print *,'Erreur lecture numerotation globale '
81 call efexit(-1)
82 endif
83
84
85
86
87 do i=1,min(nnoe,100)
88 print '(A,I3,A,I4)',
89 & 'Numero global du noeud ',i,' : ',numglb(i)
90 enddo
91
92
93
95 print '(I1)',cret
96 if (cret .ne. 0 ) then
97 print *,'Erreur fermeture du fichier'
98 call efexit(-1)
99 endif
100
101 end
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
subroutine mmhnme(fid, name, numdt, numit, entype, geotype, datype, cmode, chgt, tsf, n, cret)
Cette routine permet de lire le nombre d'entités dans un maillage pour une étape de calcul donnée.
subroutine mmhmii(fid, it, name, sdim, mdim, mtype, desc, dtunit, stype, nstep, atype, aname, aunit, cret)
Cette routine permet de lire les informations relatives à un maillage dans un fichier.
subroutine mmhgnr(fid, name, numdt, numit, entype, geotype, num, cret)