fid = fopen('MontereyBayTopoUTM.asc');
fgetl(fid)
fgetl(fid)
fgetl(fid)
fgetl(fid)
fgetl(fid)
fgetl(fid)
i=1:10
depth(i:ncols) = 