fid=fopen('dataset.txt','r');
k=0;
tline = fgetl(fid);
tline = fgetl(fid);
tline = fgetl(fid);
while 1
  tline = fgetl(fid);
  if ~ischar(tline), break, end
    k=k+1;
    staId{k}=tline(1:14);
    tmp    = sscanf(tline(14:end),'%f');
    lon(k) = tmp(4);
    lat(k) = tmp(2);
    obs(k) = tmp(7);
    elev(k)= tmp(6);
    itc(k) = tmp(8);
    otc(k) = tmp(9);
    freeAnomaly(k)=tmp(10);
    bouguerAnomaly(k)=tmp(11);
end
fclose(fid);