%   mfile initplot
%   set up graphics
%   Assign constants
%   Used to refer to all data
       ALL = (1:d);
%   Used to refer to ``real'' data, i.e. without boundaries
       OCEAN = (2:d-1);
%   ZZ stores locations of each cell in ocean as complex number
       II = sqrt(-1);
       [XX,YY]=meshgrid(ALL,ALL);
       ZZ=XX+II*YY;
       ZZocean=ZZ(OCEAN,OCEAN);