for k=04 %!!!!!!!
   
[ye, mo, da, ho, mi, se, x, y, z, d1, d2] = textread([ '2012110' num2str(k) '000000.pc'], '%s %s %s %s %s %s %f %f %f %f %f'); %!!!!!!!

end

n=length(x);

dt=strcat(ye,mo,da,ho,mi,se);

dtc=datenum(dt,'yyyymmddHHMMSS.FFF');

time=(dtc-dtc(1))*24;

figure('Position',[100,100,900,600],...
'PaperPosition',[1,1,10,6.5],...
'PaperOrientation','Landscape');

subplot(3,1,1)
plot(time,x,'.'), xlim([0 24]), grid on%ylim([-100 50]), grid on
ylabel('H_x (nT)')
title(['DIONYSOS, from ' datestr(dtc(1),'dd-mm-yyyy HH:MM:SS') ' to ' datestr(dtc(n),'dd-mm-yyyy HH:MM:SS')]) %!!!!!!!
set(gca,'xtick',0:2:24)
hold off;

subplot(3,1,2)
plot(time,y,'.'), xlim([0 24]), grid on
ylabel('H_y (nT)')
set(gca,'xtick',0:2:24)
hold off;

subplot(3,1,3)
plot(time,z,'.'), xlim([0 24]), grid on
xlabel('Time (hours)')
ylabel('H_z (nT)')
set(gca,'xtick',0:2:24)

set(gca,'LineWidth',1)
set(findobj('Type','line'),'Linewidth',2)