close all
fprintf(1,'Reading Arcgrid ...\n');
%[Z,R] = arcgridread('montbays15asc.asc');

x1 = R(3,1);
y1 = R(3,2);

[nRows, NCols] = size(Z);
mbs15.x = x1:15:x1 + ((nCols-1)*15);
mbs15.y = y1:-15:y1 + ((nRows-1)*-15);

 