clear all; close all; clc

%% Guardamos todos los nombres de los txt de la carpeta "SOHO_2012" en una estructura 'Nombres'

path='C:\Users\win8.1\Documents\Fcien\Astronomia\Profe\CTE I 2020\2020\Practica 5\SOHO_2012';
type='*.txt';
TXT=dir(fullfile(path,type));
display('ARCHIVOS: ')
for j=1:length(TXT)
    Nombres{j}=[TXT(j).name];
    display([num2str(j) '. ' Nombres{j}])
end

%% Usamos la funcion 'leerSOHO.m' para extraer los datos de velocidad de cada uno de esos txt

% La funcion leerSOHO.m devuelve vectores v, a y m, solo necesitamos las velocidades:
Velocidad(1)=0;
for j=1:length(TXT)
    [v{j} a m]=leerSOHO([Nombres{j}]);
    k=1;
    for i=length(Velocidad):length(Velocidad)+length(v{j})-1
        Velocidad(i)=v{j}(k); % Vector que contendra las velocidades de todas las CME de 2012
        k=k+1;
    end
end

%% Contamos la cantidad de CME con velocidades en distintos rangos
paso=150;
RangoVel=[0:paso:max(Velocidad)+paso*2];
for i=1:length(RangoVel)-1
    n=0;
for j=1:length(Velocidad)
    elemento=Velocidad(j);
    if  ((elemento < RangoVel(i+1))&&(elemento > RangoVel(i)))
        N(i)=n+1;
        n=N(i);
    end
end
end

RangoVel(1)=[]; 
RangoVel(length(RangoVel))=[];

%% Obtenemos la curva

plot(RangoVel,N,'*r')
xlabel('Velocidad (km/s)')
ylabel('Numero de eventos')
title('Cantidad de CME durante el 2012 en funcion de su velocidad')

% Especificamos limites y escala en el eje Y para una mejor comparacion
ylim([0 1000])
set(gca, 'YScale', 'log')
hold on
x=[0:0.1:3000];
y=10*ones(1,length(x));
plot(x,y,'-k')
hold on
y=100*ones(1,length(x));
plot(x,y,'-k')

y=[0:1000];
x=500*ones(1,length(y));
plot(x,10.^y,'-k')
x=1000*ones(1,length(y));
plot(x,10.^y,'-k')
x=2000*ones(1,length(y));
plot(x,10.^y,'-k')