for k=24:25 %!!!!!!!
   
[yep{k}, mop{k}, dap{k}, hop{k}, mip{k}, sep{k}, xp{k}, yp{k}, zp{k} d1p{k} d2p{k}] = textread([ '201210' num2str(k) '000000.fl'], '%s %s %s %s %s %s %f %f %f %f %f'); %!!!!!!!

end

% d=[dp{1};dp{2};dp{3};dp{4};dp{5};dp{6};dp{7};dp{8};dp{9};dp{10};dp{11};...
% dp{12};dp{13};dp{14};dp{15};dp{16};dp{17};dp{18};dp{19};dp{20};dp{21};...
% dp{22};dp{23};dp{24};dp{25};dp{26};dp{27};dp{28};dp{29};dp{30};dp{31}]; %!!!!!!!
% 
% t=[tp{1};tp{2};tp{3};tp{4};tp{5};tp{6};tp{7};tp{8};tp{9};tp{10};tp{11};...
% tp{12};tp{13};tp{14};tp{15};tp{16};tp{17};tp{18};tp{19};tp{20};tp{21};...
% tp{22};tp{23};tp{24};tp{25};tp{26};tp{27};tp{28};tp{29};tp{30};tp{31}]; %!!!!!!!
% 
% x=[xp{1};xp{2};xp{3};xp{4};xp{5};xp{6};xp{7};xp{8};xp{9};xp{10};xp{11};...
% xp{12};xp{13};xp{14};xp{15};xp{16};xp{17};xp{18};xp{19};xp{20};xp{21};...
% xp{22};xp{23};xp{24};xp{25};xp{26};xp{27};xp{28};xp{29};xp{30};xp{31}]; %!!!!!!!
% 
% y=[yp{1};yp{2};yp{3};yp{4};yp{5};yp{6};yp{7};yp{8};yp{9};yp{10};yp{11};...
% yp{12};yp{13};yp{14};yp{15};yp{16};yp{17};yp{18};yp{19};yp{20};yp{21};...
% yp{22};yp{23};yp{24};yp{25};yp{26};yp{27};yp{28};yp{29};yp{30};yp{31}]; %!!!!!!!
% 
% z=[zp{1};zp{2};zp{3};zp{4};zp{5};zp{6};zp{7};zp{8};zp{9};zp{10};zp{11};...
% zp{12};zp{13};zp{14};zp{15};zp{16};zp{17};zp{18};zp{19};zp{20};zp{21};...
% zp{22};zp{23};zp{24};zp{25};zp{26};zp{27};zp{28};zp{29};zp{30};zp{31}]; %!!!!!!!

ye=[yep{24};yep{25}]; %!!!!!!!

mo=[mop{24};mop{25}]; %!!!!!!!

da=[dap{24};dap{25}]; %!!!!!!!

ho=[hop{24};hop{25}]; %!!!!!!!

mi=[mip{24};mip{25}]; %!!!!!!!

se=[sep{24};sep{25}]; %!!!!!!!

x=[xp{24};xp{25}]; %!!!!!!!

y=[yp{24};yp{25}]; %!!!!!!!

z=[zp{24};zp{25}]; %!!!!!!!

n=length(x);

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

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

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

dst=[9 6 3 1 2 6 2 2 -1 -1 6 13 18 16 12 9 8 8 26 27 22 31 3 -55 -103 ... 
    -134 -116 -119 -122 -123 -122 -112 -100 -84 -73 -58 -50 -48 -47 -47 ... 
    -49 -49 -54 -48 -44 -37 -31 -31];

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

subplot(3,1,1)
plot(1:48,dst,'o-'), xlim([0 time(n)]), grid on
ylabel('Dst index (nT)')
title(['DIONYSOS, from ' datestr(dtc(1),'dd-mm-yyyy HH:MM:SS') ' to ' datestr(dtc(n))]) %!!!!!!!
set(gca,'xtick',0:2:48)
hold off;

subplot(3,1,2)
plot(time,x,'.'), xlim([0 time(n)]), ylim([-100 100]), grid on
ylabel('H_x (nT)')
set(gca,'xtick',0:2:48)
hold off;

subplot(3,1,3)
plot(time,y,'.'), xlim([0 time(n)]), grid on
xlabel('Time (hours)')
ylabel('H_y (nT)')
set(gca,'xtick',0:2:48)
% hold off;
% 
% subplot(4,1,4)
% plot(time,z-mean(z),'.'), xlim([0 time(n)]), grid on
% xlabel('Time (hours)')
% ylabel('H_z (nT)')
% set(gca,'xtick',0:2:48)

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