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' subroutine mfiope(fid, name, access, cret)
subroutine msmnax(fid, it, naxis, cret)
subroutine msmsmi(fid, it, name, sdim, mdim, desc, atype, aname, aunit, cret)
subroutine msmnsm(fid, n, cret)
subroutine mficlo(fid, cret)