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)
subroutine msmsni(fid, name, sdim, mdim, desc, atype, aname, aunit, cret)
subroutine msmnan(fid, name, naxis, cret)
subroutine mfiope(fid, name, access, cret)