%clearvars -except momold
clear data
clear mom
clear error

nu=1:2:15;
N=50*2.^[0:5];
M=1:5;
nv=length(vars);

delimiter = ' ';

formatSpec = '%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%f%[^\n\r]';
formatSpec2 = '%f%f%f%f%f%[^\n\r]';

itest=0;
for g=1:length(nu)
    for n=1:length(N)
        for m=1:length(M)
           
            itest=itest+1;            
            
            for i=1:nv
                if i~=nv
                    folder=[testset,'/nu',num2str(g),...
                        '_N',num2str(n),'_M',num2str(m),'/'];
                else
                    folder=[testsetb,'/nu',num2str(g),...
                        '_N',num2str(n),'_M',num2str(m),'/'];
                end
                filename = [folder,vars{i},'.txt'];
 
                fileID = fopen(filename,'r');

                dataArray = textscan(fileID, formatSpec, 'Delimiter', delimiter, ...
                    'MultipleDelimsAsOne', true,  'ReturnOnError', false);
                fclose(fileID);
                
                for j=1:10
                    mom(j).c(:,itest,i)=NaN;
                    mom(j).r(:,itest,i)=NaN;
                    mom(j).f(:,itest,i)=NaN;
            
                    mom(j).c(1:length(dataArray{j})-1,itest,i) = dataArray{j+3}(2:end);
                    mom(j).r(1:length(dataArray{j})-1,itest,i) = dataArray{j+3+16}(2:end);
                    
                    mom(j).c(mom(j).c<=0)=NaN;
                    mom(j).r(mom(j).r<=0)=NaN;
                    mom(j).f = mom(j).c+mom(j).r;
                end
           
                 if mom(7).r(1,itest,1)>.05e-13 %check for improper initial conditions
                    display('large error at beginning')
                 end

                info.nu(itest)=g;
                info.N(itest)=n;
                info.M(itest)=m;
            end
        end
    end
end
ntest=itest;