

function ffs=polte(ky1,dky,nky,te,rho,ff)

[m,M]=size(ff);

i_l=8;
i_r=8192;
xintp=intp_coef(i_l,i_r);

  kyy=rho(:)*sin(te);
  kyy=kyy.';
  ff=ff.';

  dky0=dky/i_r;

%%INTERPOLATION

    ffs=zeros(nky,m);

    for j=1:M;
	j	
	fky=kyy(j,:)-ky1;
	iky=floor(fky./(dky*ones(1,m)));
	dif=fky-(dky*iky);
	iky=iky+1;
	iky=ones(i_l,1)*iky+[-floor(i_l/2)+1:floor(i_l/2)]'*ones(1,m);
	iky=max(iky,ones(i_l,m));
	iky=min(iky,nky*ones(i_l,m));
	iky=iky+nky*ones(i_l,1)*[0:m-1];
	iky0=round(dif/dky0)+1;

	ffs(iky)=ffs(iky)+ ((ones(i_l,1)*ff(j,:)).*xintp(:,iky0));
    end;

  ffs=ffs.';
  ff=ff.';
