load /data5/jmanning/bathy/gom15
id=find(lon>-69&lon<-67);
id2=find(lat>40.5&lat<41.5);
x=lon(id);y=lat(id2);z=z(id2,id);
xg=ones(239,1)*x;
y=y';
yg=y*ones(1,479);
xgom=reshape(xg,239*479,1);
ygom=reshape(yg,239*479,1);
dgom=reshape(z,239*479,1);
xgom=xgom(1:4:length(xgom));
ygom=ygom(1:4:length(ygom));
dgom=dgom(1:4:length(dgom));
save /home3/ocn/jmanning/atsee/w_depthsf xgom ygom dgom
