error('this should not be tried! Takes a while');
sti='d:/Data/ibcao/';
file='IBCAO_beta_lonlat.asc';
cd(sti)
load(file)
lim=[-10 10 59 84];

ii=find(...
IBCAO_beta_lonlat(:,1)>=lim(1)&...
IBCAO_beta_lonlat(:,1)<=lim(2)&...
IBCAO_beta_lonlat(:,2)>=lim(3)&...
IBCAO_beta_lonlat(:,2)<=lim(4));

IBCAO_gin=IBCAO_beta_lonlat(ii,:);

clear IBCAO_beta_lonlat

save IBCAO_gin IBCAO_gin

lim=[-2 6 64 68];

ii=find(...
IBCAO_gin(:,1)>=lim(1)&...
IBCAO_gin(:,1)<=lim(2)&...
IBCAO_gin(:,2)>=lim(3)&...
IBCAO_gin(:,2)<=lim(4));

IBCAO_owsm=IBCAO_gin(ii,:);

clear IBCAO_gin

save IBCAO_owsm IBCAO_owsm


fig IBCAO-test; clf;
m_proj('utm','lon',lim(1:2),'lat',lim(3:4))
m_plot(IBCAO_owsm(:,1),IBCAO_owsm(:,2),'.','markersize',2);
m_grid;

lim=[-2 6 64.0 68.0];
fig ETOPO5-test;clf
m_proj('utm','lon',lim(1:2),'lat',lim(3:4));
[elev,long,lat]=m_tbase(lim);
m_plot(long,lat,'.r','markersize',2);
m_grid;