input_dir = '/home/eleanor/data/etopo2/';
filename = 'faroeice_4727.xyz';
%filename='wholeoce_2944.xyz';

lat_lim = [56 67];
lon_lim = [-20 0];
%lat_lim = [-90 90];
%lon_lim = [-180 180];
if abs(lon_lim(1)-lon_lim(2))>40
x=input('Are you sure?');
end
[lat,lon,depth] = read_xyz(input_dir,filename,lat_lim,lon_lim);

if(0)
     [LAT LON]=meshgrid(lat,lon);
m_proj('UTM','long',[-18 -6],'lat',[60 66]);

m_elev('contourf',[500:500:6000]);
m_grid('box','fancy','tickdir','in');
m_pcolor(LON,LAT,depth);
shading flat

h=m_line([-10.5 -10.5],[62.3 62.7]);
set(h,'linewidth',2)
set(h,'color','k')
h=m_line([-10.5 -9.5],[62.7 62.7]);
set(h,'linewidth',2)
set(h,'color','k')
h=m_line([-10.5 -9.5],[62.3 62.3]);
set(h,'linewidth',2)
set(h,'color','k')
h=m_line([-9.5 -9.5],[62.3 62.7]);
set(h,'linewidth',2)
set(h,'color','k')
end