30 character*64 fname,lname1,giname1,isname1
31 character*64 giname,isname
32 parameter(fname=
"Unittest_MEDlocalization_1.med")
33 parameter(lname1 =
"Localization name")
34 parameter(giname1=med_no_interpolation)
35 parameter(isname1=med_no_mesh_support)
36 integer gtype1,sdim1,nip1
37 integer gtype,sdim,nip
38 parameter(gtype1=med_tria3)
41 real*8 ecoo1(6), ipcoo1(6), wght1(3)
42 real*8 ecoo(6), ipcoo(6), wght(3)
43 data ecoo1 / 0.0, 0.0, 1.0, 0.0, 0.0,1.0 /
44 data ipcoo1 / 0.166666, 0.166666, 0.66666, 0.166666,
45 & 0.166666, 0.666666 /
46 data wght1 / 0.166666, 0.166666, 0.166666 /
49 integer sgtype,sgtype1
50 parameter(sgtype1=med_undef_geotype)
54 call mfiope(fid,fname,med_acc_rdonly,cret)
56 if (cret .ne. 0 )
then 57 print *,
'ERROR : open file' 63 call mlclni(fid, lname1, gtype, sdim, nip,
64 & giname, isname, nsmc, sgtype, cret)
66 if (cret .ne. 0 )
then 67 print *,
'ERROR : read information' 70 if ((gtype .ne. gtype1) .or.
71 & (sdim .ne. sdim1) .or.
72 & (nip .ne. nip1) .or.
73 & (giname .ne. giname1) .or.
74 & (isname .ne. isname1) .or.
75 & (nsmc .ne. nsmc1) .or.
76 & (sgtype .ne. sgtype1) )
then 78 print *,gtype1,sdim1,nip1,
"|",giname1,
"|",
"|",
79 & isname1,
"|",nsmc1,sgtype1
80 print *,gtype,sdim,nip,
"|",giname,
"|",
"|",isname,
"|",
82 print *,
'ERROR : read information' 88 call mlclor(fid,lname1,med_full_interlace,
89 & ecoo,ipcoo,wght,cret)
91 if (cret .ne. 0 )
then 92 print *,
'ERROR : read localization' 96 if ((ecoo(1) .ne. ecoo1(1)) .or.
97 & (ecoo(2) .ne. ecoo1(2)) .or.
98 & (ecoo(3) .ne. ecoo1(3)) .or.
99 & (ecoo(4) .ne. ecoo1(4)) .or.
100 & (ecoo(5) .ne. ecoo1(5)) .or.
101 & (ecoo(6) .ne. ecoo1(6)))
then 102 print *,
'ERROR : read localization' 106 if ((ipcoo(1) .ne. ipcoo1(1)) .or.
107 & (ipcoo(2) .ne. ipcoo1(2)) .or.
108 & (ipcoo(3) .ne. ipcoo1(3)) .or.
109 & (ipcoo(4) .ne. ipcoo1(4)) .or.
110 & (ipcoo(5) .ne. ipcoo1(5)) .or.
111 & (ipcoo(6) .ne. ipcoo1(6)))
then 112 print *,
'ERROR : read localization' 116 if ((wght(1) .ne. wght1(1)) .or.
117 & (wght(2) .ne. wght1(2)) .or.
118 & (wght(3) .ne. wght1(3)))
then 119 print *,
'ERROR : read localization' 127 if (cret .ne. 0 )
then 128 print *,
'ERROR : close file' subroutine mficlo(fid, cret)
subroutine mlclni(fid, lname, gtype, sdim, nip, giname, isname, nsmc, sgtype, cret)
subroutine mlclor(fid, lname, swm, ecoo, ipcoo, wght, cret)
subroutine mfiope(fid, name, access, cret)