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'
subroutine mfiope(fid, name, access, cret)
Ouverture d'un fichier MED.
subroutine mficlo(fid, cret)
Fermeture d'un fichier MED.
subroutine mlclni(fid, lname, gtype, sdim, nip, giname, isname, nsmc, sgtype, cret)
Cette routine permet d'obtenir la description d'une localisation de points d'intégration nommée local...
subroutine mlclor(fid, lname, swm, ecoo, ipcoo, wght, cret)