
myray=2160  & myray2=myray*2   & myraym1 = myray-1
totnum=float(myray*2.0)
fakex=findgen(myray)
acirc = !PI * 2.0 / float(myray)
xcen = 255.5  & ycen = 255.5   & pixrad = 104.44  & rnder = 0.5005
dp = findgen(myray) * acirc  +  acirc * (0.25 * (float(myray)))
mydp = findgen(myray) * acirc   * 2.0
dpx = intarr(myray)
dpy = intarr(myray)

 
B1sdp = sin(mydp)
B2sdp = sin(mydp) + sin(mydp*2.)
S1fft = fft(B1sdp,-1)
S2fft = fft(B2sdp,-1)
 plot,[0,10],[0,1.],/nodata
oplot,abs(S1fft(0:10))
oplot,abs(S2fft(0:10)),linestyle=2., thick = 10

plot, B2sdp
newfft=S1fft*S2fft*2.
plot,abs(newfft(0:10))
newB=fft(newfft,1)
plot,newB

