function  read_block, image, i, j, SIDE=side

;   side of a square block
if not(keyword_set(side)) then  side = 32

;   reading of the block
block = image[side*i:side*(i+1)-1, side*j:side*(j+1)-1]


return,  block
end