%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Project Name: RICO 2005 
%Purpose: Reads the aerosol_nohds.txt files and plots data
%
%Input Data: aerosol text files
%
%Last Modified: Kristen Rasmussen 6/24/05
%
%%%%%%%%%%READS THE TXT FILES %%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all;
fclose all;

addpath 'C:\rico'
addpath 'C:\rico\flux'

pathdata = 'C:\rico\aerosol';
outpath  = 'C:\rico\aerosol';
eval(['cd ', pathdata]);


x=load('C:\rico\aerosol\data\aerosol_nohds.txt');%read file with 1-min average data; 

jdy=x(:,2); %julian day at beginning of 5-min average
time=x(:,3); %time of sample (hhmmss)
int=x(:,4); %sample interval (minutes)
vol=x(:,5); %sample volume 
totcts=x(:,6); %total counts
sz1=x(:,7); %size 1 (microns)
cts1=x(:,8); %counts for size 1
sz2=x(:,9); %size 2 (microns)
cts2=x(:,10); %counts for size 2
sz3=x(:,11); %size 3 (microns)
cts3=x(:,12); %counts for size 3
sz4=x(:,13); %size 4 (microns)
cts4=x(:,14); %counts for size 4
sz5=x(:,15); %size 5 (microns)
cts5=x(:,16); %counts for size 5
sz6=x(:,17); %size 6 (microns)
cts6=x(:,18); %counts for size 6
convr=0.028*10.^5
%%%%%%%%%%%%%%Graphs%%%%%%%%%%%%%%%%%%%%%%%%

figure(1)
plot(totcts,'.');
set(gcf,'DefaultAxesFonts',10)
xlabel('Julian Day','FontName','times','FontSize',14);
ylabel('Aerosol Concentration per cc','FontName','times','FontSize',14);
title('RICO 05 R/V Seward Johnson Aerosol Concentration','FontName','times','FontSize',16);
% 
% figure(2)
% plot(jdy,org1);
% axis([19.54 19.7 0 50])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Rainrate from org1 (mm/hr)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Rainrate','FontName','times','FontSize',16);
% 
% figure(3)
% plot(jdy,ta);
% axis([9 24 20 30])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Air Temperature (C)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Air Temperature','FontName','times','FontSize',16);
% 
% figure(4)
% plot(jdy,ta);
% axis([19.54 19.7 22 26])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Air Temperature (C)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Air Temperature','FontName','times','FontSize',16);
% 
% figure(5)
% plot(jdy,urel);
% axis([9 24 0 20])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Relative Wind Speed (m/s)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Relative Wind Speed','FontName','times','FontSize',16);
% 
% figure(6)
% plot(jdy,urel);
% axis([19.54 19.7 0 12])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Relative Wind Speed (m/s)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Relative Wind Speed','FontName','times','FontSize',16);
% 
% figure(7)
% plot(jdy,reldir);
% axis([9 24 0 360])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Relative Wind Direction (degrees)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Relative Wind Direction','FontName','times','FontSize',16);
% 
% figure(8)
% plot(jdy,reldir);
% axis([19.54 19.7 0 360])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Relative Wind Direction (degrees)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Relative Wind Direction','FontName','times','FontSize',16);
% 
% figure(9)
% plot(jdy,head);
% axis([9 24 0 360])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Ship Heading (degrees)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Ship Heading','FontName','times','FontSize',16);
% 
% figure(10)
% plot(jdy,head);
% axis([19.54 19.7 0 360])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Ship Heading (degrees)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Ship Heading','FontName','times','FontSize',16);
% 
% figure(11)
% plot(jdy,qa);
% axis([9 24 12 18])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Air Specific Humidity (g/kg)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Air Specific Humidity','FontName','times','FontSize',16);
% 
% figure(12)
% plot(jdy,qa);
% axis([19.54 19.7 13 17])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Air Specific Humidity (g/kg)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Air Specific Humidity','FontName','times','FontSize',16);
% 
% figure(13)
% plot(jdy,rs);
% axis([9 24 0 1000])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Downward Solar Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Downward Solar Flux','FontName','times','FontSize',16);
% 
% figure(14)
% plot(jdy,rs);
% axis([19.54 19.7 0 1000])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Downward Solar Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Downward Solar Flux','FontName','times','FontSize',16);
% 
% figure(15)
% plot(jdy,rl);
% axis([9 24 300 500])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Downward IR Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Downward IR Flux','FontName','times','FontSize',16);
% 
% figure(16)
% plot(jdy,rl);
% axis([19.54 19.7 350 450])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Downward IR Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Downward IR Flux','FontName','times','FontSize',16);
% 
% figure(17)
% plot(jdy,Hsb);
% axis([9 24 0 150])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Bulk Sensible Heat Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Bulk Sensible Heat Flux','FontName','times','FontSize',16);
% 
% figure(18)
% plot(jdy,Hsb);
% axis([9 24 0 150])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Bulk Sensible Heat Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Bulk Sensible Heat Flux','FontName','times','FontSize',16);
% 
% figure(19)
% plot(jdy,Hlb);
% axis([9 24 0 700])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Bulk Latent Heat Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Bulk Latent Heat Flux','FontName','times','FontSize',16);
% 
% figure(20)
% plot(jdy,Hlb);
% axis([19.54 19.7 50 200])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Bulk Latent Heat Flux (W/m^2)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Bulk Latent Heat Flux','FontName','times','FontSize',16);
% 
% figure(21)
% plot(jdy,tsnk);
% axis([9 24 20 40])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Sea Snake Temperature (C)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Sea Snake Temperature','FontName','times','FontSize',16);
% 
% figure(22)
% plot(jdy,tsnk);
% axis([19.54 19.7 25 27])
% set(gcf,'DefaultAxesFonts',10)
% xlabel('Julian Day','FontName','times','FontSize',14);
% ylabel('Sea Snake Temperature (C)','FontName','times','FontSize',14);
% title('RICO 05 R/V Seward Johnson Sea Snake Temperature','FontName','times','FontSize',16);
% 





