31 parameter(fname =
"Unittest_MEDsupportMesh_1.med")
34 parameter(sdim1=2, mdim1=2)
36 parameter(sdim2=3,mdim2=2)
37 parameter(smname1 =
"supportMesh1")
39 parameter(smname2 =
"supportMesh2")
40 character*200 description1
41 parameter(description1=
"support mesh1 description")
42 character*200 description2
43 parameter(description2=
"support mesh2 description")
44 character*16 nomcoo2D(2)
45 character*16 unicoo2D(2)
46 data nomcoo2d /
"x",
"y"/, unicoo2d /
"cm",
"cm"/
47 character*16 nomcoo3D(3)
48 character*16 unicoo3D(3)
49 data nomcoo3d /
"x",
"y",
"z"/, unicoo3d /
"cm",
"cm",
"cm"/
50 integer atype1, atype2
51 parameter(atype1=med_cartesian, atype2=med_cartesian)
54 character*16 aunit(3), aname(3)
55 character*200 description
56 integer sdim, mdim, atype
60 call mfiope(fid,fname,med_acc_rdonly,cret)
61 print *,
'Open file in RD_ONLY access mode',cret
62 if (cret .ne. 0 )
then
63 print *,
'ERROR : open file in READ_ONLY access mode'
69 call msmnan(fid,smname1,sdim,cret)
70 print *,
'Number of axis (by name) : ',sdim
71 if (cret .ne. 0 )
then
72 print *,
'ERROR : read number of axis (by name)'
75 if (sdim .ne. sdim1)
then
76 print *,
'ERROR : number of axis (by name)'
80 call msmnan(fid,smname2,sdim,cret)
81 print *,
'Number of axis (by name) : ',sdim
82 if (cret .ne. 0 )
then
83 print *,
'ERROR : read number of axis (by name)'
86 if (sdim .ne. sdim2)
then
87 print *,
'ERROR : number of axis (by name)'
93 call msmsni(fid,smname1,sdim,mdim,
96 print *,
'Support mesh information by name',cret
97 if (cret .ne. 0 )
then
98 print *,
'ERROR : read support mesh information by name'
101 if ((sdim .ne. sdim1) .or.
102 & (mdim .ne. mdim1) .or.
103 & (description .ne. description1) .or.
104 & (atype .ne. atype1) .or.
105 & (aunit(1) .ne. unicoo2d(1)) .or.
106 & (aunit(2) .ne. unicoo2d(2)) .or.
107 & (aname(1) .ne. nomcoo2d(1)) .or.
108 & (aname(2) .ne. nomcoo2d(2))
110 print *,
'ERROR : support mesh information by name'
116 call msmsni(fid,smname2,sdim,mdim,
118 & aname, aunit, cret)
119 print *,
'Support mesh information by name',cret
120 if (cret .ne. 0 )
then
121 print *,
'ERROR : read support mesh information by name'
124 if ((sdim .ne. sdim2) .or.
125 & (mdim .ne. mdim2) .or.
126 & (description .ne. description2) .or.
127 & (atype .ne. atype2) .or.
128 & (aunit(1) .ne. unicoo3d(1)) .or.
129 & (aunit(2) .ne. unicoo3d(2)) .or.
130 & (aunit(3) .ne. unicoo3d(3)) .or.
131 & (aname(1) .ne. nomcoo3d(1)) .or.
132 & (aname(2) .ne. nomcoo3d(2)) .or.
133 & (aname(3) .ne. nomcoo3d(3))
135 print *,
'ERROR : support mesh information by name'
142 print *,
'Close file',cret
143 if (cret .ne. 0 )
then
144 print *,
'ERROR : close file'
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
subroutine msmsni(fid, name, sdim, mdim, desc, atype, aname, aunit, cret)
Cette routine permet de lire les informations relatives à un maillage support en précisant son nom...
subroutine msmnan(fid, name, naxis, cret)
Cette routine permet de lire dans un maillage support le nombre d'axes du repère des coordonnées des ...
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.