interpolatestate.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	interpolatestate (ne, ipkon, kon, lakon, ne0, mi, xstate, pslavsurf, nstate_, xstateini, islavsurf, islavsurfold, pslavsurfold, tieset, ntie, itiefac)

Function/Subroutine Documentation

interpolatestate()

subroutine interpolatestate	(	integer	ne,
		integer, dimension(*)	ipkon,
		integer, dimension(*)	kon,
		character*8, dimension(*)	lakon,
		integer	ne0,
		integer, dimension(*)	mi,
		real*8, dimension(nstate_,mi(1),*)	xstate,
		real*8, dimension(3,*)	pslavsurf,
		integer	nstate_,
		real*8, dimension(nstate_,mi(1),*)	xstateini,
		integer, dimension(2,*)	islavsurf,
		integer, dimension(2,*)	islavsurfold,
		real*8, dimension(3,*)	pslavsurfold,
		character*81, dimension(3,*)	tieset,
		integer	ntie,
		integer, dimension(2,*)	itiefac
	)

!
      implicit none
!
      character*8 lakon(*),lakonl
      character*81 tieset(3,*)
!
      integer ipkon(*),kon(*),ne,i,n,mi(*),indexc,ne0,indexcj,
     &  nstate_,kk,nopespring,iface,ifacej,ielemslave,ll,
     &  numpts,islavsurf(2,*),islavsurfold(2,*),ntie,itiefac(2,*)
!
      real*8 xstate(nstate_,mi(1),*),pslavsurf(3,*),pslavsurfold(3,*),
     &  xstateini(nstate_,mi(1),*)
!
      do i=1,ntie
         if(tieset(1,i)(81:81).ne.'C') cycle
         do kk=itiefac(1,i),itiefac(2,i)
            numpts=islavsurfold(2,kk+1)-islavsurfold(2,kk)
            if(numpts.gt.2) then
               call interpolateinface(kk,xstate,xstateini,numpts,
     &              nstate_,mi,islavsurf,pslavsurf,
     &              ne0,islavsurfold,pslavsurfold)
            endif
         enddo
      enddo
!     
      return