Subroutine MFSPF
C   =====================
C
C This is an interface between /mfcts/ and the local common blocks.
C The arrays in the common blocks cannot be completely filled, although
C all the relevant information is there.
C
#include "partap.inc"
#include "mfcts.inc"
#include "fwhmu.inc"
#include "mfgeox.inc"
#include "mfsigm.inc"
C
      integer i
      integer qua(8) /1,1,2,2,3,3,4,4/
      real dphi,r
      save qua
      integer nspline
c
      real cosphi,sinphi
c
      nspline=coutab(mfcts)
      ntmf=min(nspline,mxtm)
      if(nspline.lt.1) return
c
      do i=1,ntmf
       call fettab(mfcts,ID,i)
       r=mfcts_par(4)*zlt(qua(mfcts_oct),1)+mfcts_par(5)
       cosphi=cos(mfcts_avphi)
       sinphi=sin(mfcts_avphi)
       xtmf(1,i)=r*cosphi
       xtmf(2,i)=r*sinphi
       xtmf(3,i)=zlt(qua(mfcts_oct),1)
       xtmf(4,i)=atan(mfcts_par(4))
       xtmf(5,i)=mfcts_avphi
       xtmf(7,i)=mfcts_oct
       xtmf(6,i)=mfcts_p
       xtmf(8,i)=mfcts_id
       xtmf(9,i)=0.
       etmf(1,i)=cosphi**2*sig2r+xtmf(2,i)**2*sig2p
       etmf(2,i)=sinphi**2*sig2r+xtmf(1,i)**2*sig2p
       etmf(3,i)=sig2r*xtmf(3,i)**2/(xtmf(3,i)**2+r**2)**2
       etmf(4,i)=sig2p
       etmf(5,i)=0.25*xtmf(6,i)**2
      enddo
 9900 return
      end
*