k=24; %!!!!!!!

[d, t, x, y, z] = textread([ 'VLI_' num2str(k) '1011.dat'], '%s %s %f %f %f'); %!!!!!!!

n=length(x);

dt=strcat(d,t);

dtc=datenum(dt,'yyyy-mmm-ddHH:MM:SS.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]), ylim([26800 27050]), grid on
ylabel('H_x (nT)')
%title(['KLOKOTOS, from ' datestr(dtc(1),'dd-mmm-yyyy HH:MM:SS') ' to ' datestr(dtc(n))]) %!!!!!!!
title(['VELIES, from ' datestr(dtc(1),'dd-mmm-yyyy HH:MM:SS') ' to ' datestr(dtc(n))]) %!!!!!!!
set(gca,'xtick',0:2:24)
hold off;

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

subplot(3,1,3)
plot(time,z,'.'), xlim([0 24]), ylim([36120 36220]), 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)