function s=p2g2sa(lm,km,g)

  [jm,im]=size(g);
  w=fft2(g)/(im*jm);

  k=1:km;

  s(1+lm,1+km)=w(1,1);
  for l=1:lm
    s(1+lm+l,1+km)=real(w(1+l,1));
    s(1+lm-l,1+km)=imag(w(1+l,1));

	s(1+lm+l,1+km+k)=real(w(1+l,   1+k));
	s(1+lm-l,1+km-k)=imag(w(1+l,   1+k));
	s(1+lm-l,1+km+k)=real(w(1+jm-l,1+k));
	s(1+lm+l,1+km-k)=imag(w(1+jm-l,1+k));
  end
  s(1+lm,1+km+k)=real(w(1,1+k));
  s(1+lm,1+km-k)=imag(w(1,1+k));

