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 msmnsm(fid,nsmesh,cret)
70 print *,
'Read number of support mesh : ',nsmesh
72 if (cret .ne. 0 )
then 73 print *,
'ERROR : read number of support mesh' 76 if (nsmesh .ne. 2)
then 77 print *,
'ERROR : number of support mesh' 85 call msmnax(fid,i,sdim,cret)
86 if (cret .ne. 0 )
then 87 print *,
'ERROR : read number of axis ' 90 print *,
'Number of axis : ',sdim
93 if (sdim .ne. sdim1)
then 94 print *,
'ERROR : support mesh information' 99 if (sdim .ne. sdim2)
then 100 print *,
'ERROR : support mesh information' 105 call msmsmi(fid,i,smname,sdim,mdim,
107 & atype,aname,aunit,cret)
108 print *,
'Support mesh information',cret
109 if (cret .ne. 0 )
then 110 print *,
'ERROR : read support mesh information' 115 if ((sdim .ne. sdim1) .or.
116 & (mdim .ne. mdim1) .or.
117 & (description .ne. description1) .or.
118 & (atype .ne. atype1) .or.
119 & (aunit(1) .ne. unicoo2d(1)) .or.
120 & (aunit(2) .ne. unicoo2d(2)) .or.
121 & (aname(1) .ne. nomcoo2d(1)) .or.
122 & (aname(2) .ne. nomcoo2d(2))
124 print *,
'ERROR : support mesh information by name' 130 if ((sdim .ne. sdim2) .or.
131 & (mdim .ne. mdim2) .or.
132 & (description .ne. description2) .or.
133 & (atype .ne. atype2) .or.
134 & (aunit(1) .ne. unicoo3d(1)) .or.
135 & (aunit(2) .ne. unicoo3d(2)) .or.
136 & (aunit(3) .ne. unicoo3d(3)) .or.
137 & (aname(1) .ne. nomcoo3d(1)) .or.
138 & (aname(2) .ne. nomcoo3d(2)) .or.
139 & (aname(3) .ne. nomcoo3d(3))
141 print *,
'ERROR : support mesh information by name' 151 print *,
'Close file',cret
152 if (cret .ne. 0 )
then 153 print *,
'ERROR : close file' subroutine mficlo(fid, cret)
subroutine msmnsm(fid, n, cret)
subroutine msmsmi(fid, it, name, sdim, mdim, desc, atype, aname, aunit, cret)
subroutine msmnax(fid, it, naxis, cret)
subroutine mfiope(fid, name, access, cret)