pro Test_loop

rand_arr=100*randomn(10,100,100)
array_dimension=100.
THRESHOLD=200
log_file='C:\stats.txt'

;openw,1,log_file,/append

gt_thresh=where(rand_arr gt THRESHOLD,length_1)
print,size(gt_thresh)

for i=0,(length_1-1) do begin

    if ((gt_thresh(i) gt 3*array_dimension) and (gt_thresh(i) lt 97*array_dimension) $
        and ((gt_thresh(i) mod array_dimension) gt 2) and ((gt_thresh(i) mod array_dimension) lt 98)) then begin


        for radius=3,0,-1 do begin

            found_less=0

            for xy_radius=(-radius),radius do begin

                for diag_radius=(-radius+abs(xy_radius)),(radius-abs(xy_radius)) do begin

                    if ( (rand_arr(long(gt_thresh(i)+long(xy_radius*array_dimension)+diag_radius)) lt THRESHOLD) $
                        or (rand_arr(long(gt_thresh(i)+long(xy_radius*array_dimension)+diag_radius)) gt rand_arr(gt_thresh(i)))) then begin
                        found_less=1
                        break
                    endif

                endfor

          if (found_less eq 1) then begin
                    break
                endif

            endfor

        if (found_less eq 0) then begin
            act_radius=radius+1
         if act_radius gt 1 then begin
            print,'act radius',act_radius
            printf,1,log_file,FORMAT='("%d")',act_radius
         endif
            for xy_radius=(-radius),radius do begin

                for diag_radius=(-radius+abs(xy_radius)),(radius-abs(xy_radius)) do begin

              rand_arr(long(gt_thresh(i)+long(xy_radius*array_dimension)+diag_radius))=0

          endfor

         endfor
            break
       endif

    endfor

    endif

endfor

close,1

end