;return_entirepowerspectrum
;
;PURPOSE:
;       This script takes a PhaseDim x PhaseDim array of floating point 
;	ApDimbers in a phase screen and computes the power spectrum.
;RETURN:
;	  EllipVals(0,*)=High Passed Ellipticities
;         EllipVals(1,*)=High Passed Pointing Angles
;         EllipVals(2,*)=Low Passed Ellipticities
;         EllipVals(3,*)-Low Passed Pointing Angles

function return_soarpowerspectrum,XSlopes,YSlopes,$
	FittedOnly,ApDim

KMax=floor(float(ApDim)/2+.5)
KMax=ApDim/2
YSlopesFft=complexarr(ApDim,ApDim)
YSlopesFft(*,*)=fft(YSlopes(*,*))
XSlopesFft=complexarr(ApDim,ApDim)
XSlopesFft(*,*)=fft(XSlopes(*,*))
XSlopesFftRe=complexarr(ApDim,ApDim)
YSlopesFftRe=complexarr(ApDim,ApDim)

;ReArrange
XSlopesFftRe=shift(XSlopesFft,KMax,KMax)
YSlopesFftRe=shift(YSlopesFft,KMax,Kmax)

PowerSpectra=fltarr(ApDim,ApDim,3)
PowerSpectra(*,*,0)=(Abs(XSlopesFftRe))^2
PowerSpectra(*,*,1)=(Abs(YSlopesFftRe))^2

return,PowerSpectra

end