32 parameter(fname =
"Unittest_MEDinterp_1.med")
34 parameter(name1=
"Interpolation family name")
36 parameter(gtype1=med_tria3)
38 parameter(cnode1=med_false)
39 integer nvar1,maxd1,nmaxc1
40 parameter(nvar1=2,maxd1=1,nmaxc1=3)
41 integer ncoef1,ncoef2,ncoef3
42 parameter(ncoef1=3,ncoef2=1,ncoef3=1)
43 integer power1(6),power2(2),power3(2)
44 data power1 / 0,0, 1,0, 0,1 /
47 real*8 coef1(3), coef2(1), coef3(1)
48 data coef1 / 1., -1., -1. /
54 call mfiope(fid,fname,med_acc_creat,cret)
55 print *,
'Open file',cret
56 if (cret .ne. 0 )
then
57 print *,
'ERROR : file creation'
63 call mipcre(fid,name1,gtype1,cnode1,nvar1,
65 print *,
'interpolation creation',cret
66 if (cret .ne. 0 )
then
67 print *,
'ERROR : interpolation creation'
73 call mipbfw(fid,name1,1,ncoef1,power1,coef1,cret)
74 print *,
'function creation',cret
75 if (cret .ne. 0 )
then
76 print *,
'ERROR : function creation'
80 call mipbfw(fid,name1,2,ncoef2,power2,coef2,cret)
81 print *,
'function creation',cret
82 if (cret .ne. 0 )
then
83 print *,
'ERROR : function creation'
87 call mipbfw(fid,name1,3,ncoef3,power3,coef3,cret)
88 print *,
'function creation',cret
89 if (cret .ne. 0 )
then
90 print *,
'ERROR : function creation'
97 print *,
'Close file',cret
98 if (cret .ne. 0 )
then
99 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 mipcre(fid, name, gtype, cnode, nvar, maxd, nmaxc, cret)
Cette routine permet de créer une nouvelle fonction d'interpolation polynômiale nommée interpname.
subroutine mipbfw(fid, name, it, nc, pw, co, cret)