PRO sesuds,suds,out,lim ;Search in array SUDS to eliminate UDs with d < lim (") ;and output of: no.of UDs, total area of UDs, mean and stdev d ;(array OUT, col. 0 1 (pix) 2 (") 3 ). out=fltarr(4,360) pi4=!pi/4. for i=0,359 do begin w=where((suds(0,*) eq i) and (suds(1,*) gt lim),c) out(0,i)=c out(1,i)=total(nint(pi4*(suds(1,w)/0.125)^2)) sd=stdev(suds(1,w),me) out(2,i)=me out(3,i)=sd endfor end