%matlab time series of flight through model
     clf
  %   load 1CH2neutfoxA999 -ascii;
 %tim=1CH2neutfoxA999 (:,1);
 %den=1CH2neutfoxA999 (:,2);
     load C2H2neutfoxA999 -ascii;
 tim=C2H2neutfoxA999  (:,1);
 den1=C2H2neutfoxA999  (:,2);
     load C2H3CNneutfoxA999 -ascii;
 tim=C2H3CNneutfoxA999   (:,1);
 den2=C2H3CNneutfoxA999   (:,2);
     load C2H3neutfoxA999 -ascii;
 tim=C2H3neutfoxA999 (:,1);
 den3=C2H3neutfoxA999  (:,2);
     load C2H4neutfoxA999 -ascii;
 tim=C2H4neutfoxA999  (:,1);
 den4=C2H4neutfoxA999  (:,2);
semilogy(tim,den,'b-x')
%     plot(tim,den,'b-','linewidth',2)
     axis([-180,180,-50,20])
     xlabel ('Time in Days of Year 2004','fontsize',16)
     ylabel ('Density cm^-^3','fontsize',16)
     title('TA Neutral Number Densities','fontsize',16)
hold on
plot(tim,den1,'m-','linewidth',2)
hold on
plot(tim,den2,'r-','linewidth',2)
hold on
plot(tim,den3,'g-','linewidth',2)
hold on
plot(tim,den4,'k-','linewidth',1.5)

     legend('1CH2','C2H2','C2H3CN','C2H3','C2H4')
