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)
Fermeture d'un fichier MED.
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
subroutine msmnsm(fid, n, cret)
Cette routine permet de lire le nombre de maillages support dans un fichier.
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 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 ...