det_index=[5]
N=1024  & roll_angle=findgen(N)*2*!pi/N 
mapHwidth=32. & R=600. & AZ=0 
R0=256+80
psf=fltarr(N,mapHwidth*2,mapHwidth*2)
print,'R starts at ',R+R0
for j=0,63 do begin
   R1=R+R0+j
   psf1=hsi_polar_psf(R1,AZ,roll_angle,R0,det_index,mapHwidth,/plot)
   psf(*,*,j)=psf1
   print,R1,format='(I4,$)'
endfor
;save,psf,r,r0,mapHwidth,N,det_index,'psf5.sav'

end