%GRID_BND forms a boundary file for a grid

count = 1; %cell count

row = 1;   %row count

poly = zeros(12*12*6,2);

for i = 1:12
   
   for j = 1:12
      
      I = 2*i - 1;
      
      J = 2*j - 1;
      
      poly(row,:) = [count,5];
      
      poly(row + 1,:) = [I,J];
      
      poly(row + 2,:) = [I+2,J];
      
      poly(row + 3,:) = [I+2,J+2];
      
      poly(row + 4,:) = [I,J+2];
      
      poly(row + 5,:) = [I,J];
      
      count = count + 1;
      
      row = row + 6;
      
   end
   
end