Unittest_MEDinterp_3.f
Aller à la documentation de ce fichier.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022 program MEDinterp2
00023
00024 implicit none
00025 include 'med.hf'
00026
00027
00028 integer cret
00029 integer fid
00030 character*64 fname
00031 parameter (fname = "Unittest_MEDinterp_1.med")
00032 integer n,ni
00033 parameter (ni=1)
00034 integer it
00035 character *64 name1,name
00036 parameter (name1="Interpolation family name")
00037 integer gtype1,gtype
00038 parameter (gtype1=MED_TRIA3)
00039 integer cnode1,cnode
00040 parameter (cnode1=MED_FALSE)
00041 integer nvar1,maxd1,nmaxc1
00042 integer nvar,maxd,nmaxc
00043 parameter (nvar1=2,maxd1=1,nmaxc1=3)
00044 integer nbf,nbf1
00045 parameter (nbf1=3)
00046
00047
00048
00049 call mfiope(fid,fname,MED_ACC_RDONLY,cret)
00050 print *,'Open file',cret
00051 if (cret .ne. 0 ) then
00052 print *,'ERROR : file creation'
00053 call efexit(-1)
00054 endif
00055
00056
00057
00058 call mipnip(fid,n,cret)
00059 print *,'Number of interpolation',cret
00060 if (cret .ne. 0 ) then
00061 print *,'ERROR : number of interpolation'
00062 call efexit(-1)
00063 endif
00064 if (n .ne. ni) then
00065 print *,'ERROR : number of interpolation'
00066 call efexit(-1)
00067 endif
00068
00069
00070
00071 do it=1,n
00072 call mipipi(fid,it,name,gtype,cnode,
00073 & nbf,nvar,maxd,nmaxc,cret)
00074 print *,'interpolation information',cret
00075 if (cret .ne. 0 ) then
00076 print *,'ERROR : interpolation information'
00077 call efexit(-1)
00078 endif
00079
00080 if (it .eq. 1) then
00081 if ( (gtype .ne. gtype1) .or.
00082 & (cnode .ne. cnode1) .or.
00083 & (nbf .ne. nbf1) .or.
00084 & (nvar .ne. nvar1) .or.
00085 & (maxd .ne. maxd1) .or.
00086 & (nmaxc .ne. nmaxc1) ) then
00087 print *,'ERROR : interpolation information'
00088 call efexit(-1)
00089 endif
00090 endif
00091
00092 enddo
00093
00094
00095
00096
00097 call mficlo(fid,cret)
00098 print *,'Close file',cret
00099 if (cret .ne. 0 ) then
00100 print *,'ERROR : close file'
00101 call efexit(-1)
00102 endif
00103
00104
00105
00106 end
00107