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