32 parameter(fname =
"Unittest_MEDsupportMesh_1.med")
35 parameter(sdim1=2, mdim1=2)
37 parameter(sdim2=3,mdim2=2)
38 parameter(smname1 =
"supportMesh1")
40 parameter(smname2 =
"supportMesh2")
41 character*200 description1
42 parameter(description1=
"support mesh1 description")
43 character*200 description2
44 parameter(description2=
"support mesh2 description")
45 character*16 nomcoo2d(2)
46 character*16 unicoo2d(2)
47 data nomcoo2d /
"x",
"y"/, unicoo2d /
"cm",
"cm"/
48 character*16 nomcoo3d(3)
49 character*16 unicoo3d(3)
50 data nomcoo3d /
"x",
"y",
"z"/, unicoo3d /
"cm",
"cm",
"cm"/
51 integer atype1, atype2
52 parameter(atype1=med_cartesian, atype2=med_cartesian)
55 character*16 aunit(3), aname(3)
56 character*200 description
57 integer sdim, mdim, atype
61 call mfiope(fid,fname,med_acc_rdonly,cret)
62 print *,
'Open file in RD_ONLY access mode',cret
63 if (cret .ne. 0 )
then
64 print *,
'ERROR : open file in READ_ONLY access mode'
70 call msmnan(fid,smname1,sdim,cret)
71 print *,
'Number of axis (by name) : ',sdim
72 if (cret .ne. 0 )
then
73 print *,
'ERROR : read number of axis (by name)'
76 if (sdim .ne. sdim1)
then
77 print *,
'ERROR : number of axis (by name)'
81 call msmnan(fid,smname2,sdim,cret)
82 print *,
'Number of axis (by name) : ',sdim
83 if (cret .ne. 0 )
then
84 print *,
'ERROR : read number of axis (by name)'
87 if (sdim .ne. sdim2)
then
88 print *,
'ERROR : number of axis (by name)'
94 call msmsni(fid,smname1,sdim,mdim,
97 print *,
'Support mesh information by name',cret
98 if (cret .ne. 0 )
then
99 print *,
'ERROR : read support mesh information by name'
102 if ((sdim .ne. sdim1) .or.
103 & (mdim .ne. mdim1) .or.
104 & (description .ne. description1) .or.
105 & (atype .ne. atype1) .or.
106 & (aunit(1) .ne. unicoo2d(1)) .or.
107 & (aunit(2) .ne. unicoo2d(2)) .or.
108 & (aname(1) .ne. nomcoo2d(1)) .or.
109 & (aname(2) .ne. nomcoo2d(2))
111 print *,
'ERROR : support mesh information by name'
117 call msmsni(fid,smname2,sdim,mdim,
119 & aname, aunit, cret)
120 print *,
'Support mesh information by name',cret
121 if (cret .ne. 0 )
then
122 print *,
'ERROR : read support mesh information by name'
125 if ((sdim .ne. sdim2) .or.
126 & (mdim .ne. mdim2) .or.
127 & (description .ne. description2) .or.
128 & (atype .ne. atype2) .or.
129 & (aunit(1) .ne. unicoo3d(1)) .or.
130 & (aunit(2) .ne. unicoo3d(2)) .or.
131 & (aunit(3) .ne. unicoo3d(3)) .or.
132 & (aname(1) .ne. nomcoo3d(1)) .or.
133 & (aname(2) .ne. nomcoo3d(2)) .or.
134 & (aname(3) .ne. nomcoo3d(3))
136 print *,
'ERROR : support mesh information by name'
143 print *,
'Close file',cret
144 if (cret .ne. 0 )
then
145 print *,
'ERROR : close file'