annualplot.m

data=load('IPYnoonMH.dat');
day = data(:,1);
dayno = data(:,2);
month = data(:,3);
year = data(:,4);
slt = data(:,5);
ut = data(:,6);
gdalt = data(:,7);
f107=data(:,8)./1e-22;
fbar = data(:,9)./1e-22;
nel = data(:,10);
nelm = data(:,11);
% ...  te, ti, ,,,

ii=find(nel>9.5 & nel<13);

plot (dayno(ii),nel(ii),'ro')

xlabel('day number of the year')
ylabel('electron density')