function bat=makesabbat clear datafile=dir('/scratch2/bathy/*.asc'); datafile={datafile(:).name}; lo1=[-81;-82;-82;-82;-82;-82;-80;-79;-77]; lo2=[-79;-79;-79;-79;-79;-77;-76;-75;-73]; for i=1:length(datafile) ZZ=load(datafile{i});ib=find(ZZ>10000);ZZ(ib)=NaN; la=(27+i-1):6/3600:(28+i-1);la=fliplr(la); lo=lo1(i):6/3600:lo2(i); hm=contour(lo,la,ZZ);hold on; axis([-82,-73,27,36]) [x,y]=ll2xy(la,lo,31.25,-80.75); bat(i).x=x; bat(i).y=y; bat(i).z=ZZ; bat(i).battype='gridded'; end