31 parameter(fname =
"Unittest_MEDstructElement_1.med")
32 character*64 mname1, mname2, mname3
33 parameter(mname1 =
"model name 1")
34 parameter(mname2 =
"model name 2")
35 parameter(mname3 =
"model name 3")
36 integer dim1, dim2, dim3
41 parameter(smname1=med_no_name)
43 parameter(smname2=
"support mesh name")
45 parameter(setype1=med_none)
47 parameter(setype2=med_node)
49 parameter(setype3=med_cell)
51 parameter(sgtype1=med_no_geotype)
53 parameter(sgtype2=med_no_geotype)
55 parameter(sgtype3=med_seg2)
56 integer mtype1,mtype2,mtype3
67 integer ncatt1,profile1,nvatt1
75 integer mgtype,mdim,setype,snnode,sncell
76 integer sgtype,ncatt,nvatt,profile
77 character*64 smname,mname
81 call mfiope(fid,fname,med_acc_rdonly,cret)
82 print *,
'Open file',cret
83 if (cret .ne. 0 )
then 84 print *,
'ERROR : file creation' 91 print *,
'Read number of struct model',nsmr,cret
92 if (cret .ne. 0 )
then 93 print *,
'ERROR : number of struct model' 96 if (nsmr .ne. nsm)
then 97 print *,
'ERROR : number of struct model' 105 call msesei(fid,it,mname,mgtype,mdim,smname,
106 & setype,snnode,sncell,sgtype,
107 & ncatt,profile,nvatt,cret)
108 print *,
'Read information about struct element',cret
109 if (cret .ne. 0 )
then 110 print *,
'ERROR : information about struct element' 115 if ( (mname .ne. mname1) .or.
116 & (mgtype .ne. mtype1) .or.
117 & (mdim .ne. dim1) .or.
118 & (smname .ne. smname1) .or.
119 & (setype .ne. setype1) .or.
120 & (snnode .ne. nnode1) .or.
121 & (sncell .ne. ncell1) .or.
122 & (sgtype .ne. sgtype1) .or.
123 & (ncatt .ne. ncatt1) .or.
124 & (profile .ne. profile1) .or.
125 & (nvatt .ne. nvatt1)
127 print *,
'ERROR : information about struct element' 133 if ( (mname .ne. mname2) .or.
134 & (mgtype .ne. mtype2) .or.
135 & (mdim .ne. dim2) .or.
136 & (smname .ne. smname2) .or.
137 & (setype .ne. setype2) .or.
138 & (snnode .ne. nnode2) .or.
139 & (sncell .ne. ncell1) .or.
140 & (sgtype .ne. sgtype2) .or.
141 & (ncatt .ne. ncatt1) .or.
142 & (profile .ne. profile1) .or.
143 & (nvatt .ne. nvatt1)
145 print *,
'ERROR : information about struct element ' 151 if ( (mname .ne. mname3) .or.
152 & (mgtype .ne. mtype3) .or.
153 & (mdim .ne. dim3) .or.
154 & (smname .ne. smname2) .or.
155 & (setype .ne. setype3) .or.
156 & (snnode .ne. nnode2) .or.
157 & (sncell .ne. ncell2) .or.
158 & (sgtype .ne. sgtype3) .or.
159 & (ncatt .ne. ncatt1) .or.
160 & (profile .ne. profile1) .or.
161 & (nvatt .ne. nvatt1)
163 print *,
'ERROR : information about struct element' 172 call msesen(fid,mtype1,mname,cret)
173 print *,
'Read struct element name from the type',cret
174 if (cret .ne. 0 )
then 175 print *,
'ERROR : struct element name from the type' 178 if (mname .ne. mname1)
then 179 print *,
'ERROR : struct element name from the type' 183 call msesen(fid,mtype2,mname,cret)
184 print *,
'Read struct element name from the type',cret
185 if (cret .ne. 0 )
then 186 print *,
'ERROR : struct element name from the type' 189 if (mname .ne. mname2)
then 190 print *,
'ERROR : struct element name from the type' 194 call msesen(fid,mtype3,mname,cret)
195 print *,
'Read struct element name from the type',cret
196 if (cret .ne. 0 )
then 197 print *,
'ERROR : struct element name from the type' 200 if (mname .ne. mname3)
then 201 print *,
'ERROR : struct element name from the type' 208 print *,
'Close file',cret
209 if (cret .ne. 0 )
then 210 print *,
'ERROR : close file' subroutine mficlo(fid, cret)
program medstructelement3
subroutine msense(fid, n, cret)
subroutine msesei(fid, it, mname, mgtype, mdim, smname, setype, snnode, sncell, sgtype, ncatt, ap, nvatt, cret)
subroutine mfiope(fid, name, access, cret)
subroutine msesen(fid, mgtype, mname, cret)