
%filename = idm_path;
%fname = idm_file;
fid = fopen('/hao/aim2/dmsp/dmsp_idm_4s/9801/l1398015a.txt','r');
instr_date = fgets(fid);
A = textscan(fid,'%f %f %f %f %f %f %f %f %f %f %f %f %f','HeaderLines',0,'ReturnOnError',1);
A = cell2mat(A);
[m,n]=size(A);
while m < 20275
    B = textscan(fid,'%f %f %f %f %f %f %f %f %f %f %f %f %f',...
        'HeaderLines',m+1,'ReturnOnError',1);
    B = cell2mat(B);
    [i,j]=size(B);
    if i<2
        i=i+1;
    end
    newArray = zeros((m+i),13);
    newArray(1:m,:) = A;
    newArray(m+1:end,:) = B;
    clear A
    clear B
    A = newArray;
    clear newArray
    [m,n]=size(A);
    disp(m)
end
data = A;
