?clear?all;
ncload（‘etopo2.nc‘）;
sum=0;
for?i=8377:9421????
a=lon（i1）;
for?j=2728:3931
b=lat（j1）;
c=topo（ji）;
d=[a?b?c];
aaa=size（d）;
pp（sum+1:sum+aaa:）=d;
sum=sum+aaa;
end?
end
save?etopo.dat?-ascii?pp
b=find（pp（:1）>111.48?&?pp（:1）<112.26355?&pp（:2）>21.085217?&pp（:2）<24.14）;
pp（b:）=[];
c=find（pp（:1）>112.26355?&?pp（:1）<113.0711608?&pp（:2）>21.496067?&pp（:2）<24.14）;
pp（c:）=[];
d=find（pp（:1）>113.0711608?&?pp（:1）<114.4086373&pp（:2）>21.76037?&pp（:2）<24.14）;