% routine to extract data from lat/lon plot with ginput
function [id]=getidll(lo,la);
[lom,lam]=ginput(2);
id=find(lo>lom(1)&lo<lom(2));
plot(lo(id),la(id))
pause(2);
id2=find(la(id)>lam(1)&la(id)<lam(2));
plot(lo(id(id2)),la(id(id2)))
pause(2);
id=id(id2);