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'
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
subroutine msmnan(fid, name, naxis, cret)
Cette routine permet de lire dans un maillage support le nombre d'axes du repère des coordonnées des ...
subroutine msmsni(fid, name, sdim, mdim, desc, atype, aname, aunit, cret)
Cette routine permet de lire les informations relatives à un maillage support en précisant son nom.