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)
Ouverture d'un fichier MED.
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
subroutine msmsmi(fid, it, name, sdim, mdim, desc, atype, aname, aunit, cret)
Cette routine permet de lire les informations relatives à un maillage support dans un fichier.
subroutine msmnsm(fid, n, cret)
Cette routine permet de lire le nombre de maillages support dans un fichier.
subroutine msmnax(fid, it, naxis, cret)
Cette routine permet de lire dans un maillage support le nombre d'axes du repère des coordonnées des ...