function [lat,lon,depth] = read_xyz(input_dir,filename,lat_lim,lon_lim);

%input_dir = '/home/eleanor/data/etopo2/';
%filename = 'faroeice_4727.xyz';
mydata = csvread([input_dir,filename],0,0);

[N c]=size(mydata);

lon = unique(mydata(:,1));
lat = unique(mydata(:,2));

for n=1:N

   ilon = find(lon==mydata(n,1));
   ilat = find(lat==mydata(n,2));
   depth(ilon,ilat) = mydata(n,3);
end

if 0
lat_range = find(lat>=min(lat_lim)&lat<=max(lat_lim));
lon_range = find(lon>=min(lon_lim)&lon<=max(lon_lim));
surf(lon(lon_range),lat(lat_range),depth(lon_range,lat_range)');
shading flat
end