31 character*64 fname,lname1,giname1,isname1
32 character*64 giname,isname
33 parameter(fname=
"Unittest_MEDlocalization_1.med")
34 parameter(lname1 =
"Localization name")
35 parameter(giname1=med_no_interpolation)
36 parameter(isname1=med_no_mesh_support)
37 integer gtype1,sdim1,nip1
38 integer gtype,sdim,nip
39 parameter(gtype1=med_tria3)
42 real*8 ecoo1(6), ipcoo1(6), wght1(3)
43 real*8 ecoo(6), ipcoo(6), wght(3)
44 data ecoo1 / 0.0, 0.0, 1.0, 0.0, 0.0,1.0 /
45 data ipcoo1 / 0.166666, 0.166666, 0.66666, 0.166666,
46 & 0.166666, 0.666666 /
47 data wght1 / 0.166666, 0.166666, 0.166666 /
50 integer sgtype,sgtype1
51 parameter(sgtype1=med_undef_geotype)
55 call mfiope(fid,fname,med_acc_rdonly,cret)
57 if (cret .ne. 0 )
then
58 print *,
'ERROR : open file'
64 call mlclni(fid, lname1, gtype, sdim, nip,
65 & giname, isname, nsmc, sgtype, cret)
67 if (cret .ne. 0 )
then
68 print *,
'ERROR : read information'
71 if ((gtype .ne. gtype1) .or.
72 & (sdim .ne. sdim1) .or.
73 & (nip .ne. nip1) .or.
74 & (giname .ne. giname1) .or.
75 & (isname .ne. isname1) .or.
76 & (nsmc .ne. nsmc1) .or.
77 & (sgtype .ne. sgtype1) )
then
79 print *,gtype1,sdim1,nip1,
"|",giname1,
"|",
"|",
80 & isname1,
"|",nsmc1,sgtype1
81 print *,gtype,sdim,nip,
"|",giname,
"|",
"|",isname,
"|",
83 print *,
'ERROR : read information'
89 call mlclor(fid,lname1,med_full_interlace,
90 & ecoo,ipcoo,wght,cret)
92 if (cret .ne. 0 )
then
93 print *,
'ERROR : read localization'
97 if ((ecoo(1) .ne. ecoo1(1)) .or.
98 & (ecoo(2) .ne. ecoo1(2)) .or.
99 & (ecoo(3) .ne. ecoo1(3)) .or.
100 & (ecoo(4) .ne. ecoo1(4)) .or.
101 & (ecoo(5) .ne. ecoo1(5)) .or.
102 & (ecoo(6) .ne. ecoo1(6)))
then
103 print *,
'ERROR : read localization'
107 if ((ipcoo(1) .ne. ipcoo1(1)) .or.
108 & (ipcoo(2) .ne. ipcoo1(2)) .or.
109 & (ipcoo(3) .ne. ipcoo1(3)) .or.
110 & (ipcoo(4) .ne. ipcoo1(4)) .or.
111 & (ipcoo(5) .ne. ipcoo1(5)) .or.
112 & (ipcoo(6) .ne. ipcoo1(6)))
then
113 print *,
'ERROR : read localization'
117 if ((wght(1) .ne. wght1(1)) .or.
118 & (wght(2) .ne. wght1(2)) .or.
119 & (wght(3) .ne. wght1(3)))
then
120 print *,
'ERROR : read localization'
128 if (cret .ne. 0 )
then
129 print *,
'ERROR : close file'