close all


path='C:\Program Files\MATLAB71\work';
control_data=[];


for ddd=11;
    dd=num2str(ddd);
    if ddd<10;
       dd=['0' num2str(ddd)];	

    for hhh=0:23;
        hr=num2str(hhh)
        if hhh<10;
           hr=['0' num2str(hhh)];	

    for mmm=7;
        mm=num2str(mmm);
        if mmm<10;
        dd=['0' num2str(mmm)];	
        

    for yyy=7;
        yy=num2str(yyy);
        if yyy<10;
           yy=['0' num2str(yyy)];	

    for minu=0
        min=num2str(minu);
        if  minu<10;
            min=['0' num2str(minu)];	

control_imp=importda([path, '\O3_control\FastO3\',yy,'-',mm,'-',dd,'\control_' , yy, '-' , mm, '-',dd,'_',hr ,min, '.lvm']);
control_data=[control_data;control_imp]	;

                end;
            end;
        end;
    end;
end;


time=control_data(1,:);
sf=control_data(2,:);
NOf=control_data(3,:);
AF=control_data(4,:);
PMT=control_data(5,:);
RV=control_data(6,:);
caseT=control_data(7,:);
P=control_data(8,:);
ppb=control_data(9,:);
count=control_data(10,:);
PSL=control_data(11,:);
drier=control_data(12,:);
ML3=control_data(13,:);


figure; plot(time, PMT, '-'); grid; xlbel('time / s'); ylabel('T /°C'); title('PMT');
figure; plot(time, RV, '-'); grid; xlbel('time / s'); ylabel('T / °C'); title('reaction vessel');
figure; plot(time, caseT, '-'); grid; xlbel('time / s'); ylabel('T / °C'); title('case temperature');
figure; plot(time, ppb, '-'); grid; xlbel('time / s'); ylabel('ppb'); title('O3 Monitor, ppb');
figure; plot(time, count, '-'); grid; xlbel('time / s'); ylabel('counts / S'); title('fast O3, counts');
figure; plot(time, ML3, '-'); grid; xlbel('time / s'); ylabel('ppb'); title('ML#3');
figure; plot(time, P, '-'); grid; xlbel('time / s'); ylabel('pressure / torr'); title('Pressure');