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 msmnsm(fid,nsmesh,cret)
71 print *,
'Read number of support mesh : ',nsmesh
73 if (cret .ne. 0 )
then
74 print *,
'ERROR : read number of support mesh'
77 if (nsmesh .ne. 2)
then
78 print *,
'ERROR : number of support mesh'
86 call msmnax(fid,i,sdim,cret)
87 if (cret .ne. 0 )
then
88 print *,
'ERROR : read number of axis '
91 print *,
'Number of axis : ',sdim
94 if (sdim .ne. sdim1)
then
95 print *,
'ERROR : support mesh information'
100 if (sdim .ne. sdim2)
then
101 print *,
'ERROR : support mesh information'
106 call msmsmi(fid,i,smname,sdim,mdim,
108 & atype,aname,aunit,cret)
109 print *,
'Support mesh information',cret
110 if (cret .ne. 0 )
then
111 print *,
'ERROR : read support mesh information'
116 if ((sdim .ne. sdim1) .or.
117 & (mdim .ne. mdim1) .or.
118 & (description .ne. description1) .or.
119 & (atype .ne. atype1) .or.
120 & (aunit(1) .ne. unicoo2d(1)) .or.
121 & (aunit(2) .ne. unicoo2d(2)) .or.
122 & (aname(1) .ne. nomcoo2d(1)) .or.
123 & (aname(2) .ne. nomcoo2d(2))
125 print *,
'ERROR : support mesh information by name'
131 if ((sdim .ne. sdim2) .or.
132 & (mdim .ne. mdim2) .or.
133 & (description .ne. description2) .or.
134 & (atype .ne. atype2) .or.
135 & (aunit(1) .ne. unicoo3d(1)) .or.
136 & (aunit(2) .ne. unicoo3d(2)) .or.
137 & (aunit(3) .ne. unicoo3d(3)) .or.
138 & (aname(1) .ne. nomcoo3d(1)) .or.
139 & (aname(2) .ne. nomcoo3d(2)) .or.
140 & (aname(3) .ne. nomcoo3d(3))
142 print *,
'ERROR : support mesh information by name'
152 print *,
'Close file',cret
153 if (cret .ne. 0 )
then
154 print *,
'ERROR : close file'