%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Project Name: RICO 2005 
%Purpose: Reads the flux_30hf_rico_05a_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);%air specific humidity, g/kg (z=13 m)
cts2=x(:,10);%downward solar flux, W/m^2 (ETL units)
sz3=x(:,11);%downward IR flux, W/m^2 (ETL units)
cts3=x(:,12);%rainrate, mm/hr (ETL STI optical rain gauge, uncorrected)
sz4=x(:,13);%ship speed, m/s (ETL gps)
cts4=x(:,14);%ship heading, deg clockwise rel north (ETL gyrocompass)
sz5=x(:,15);%relative wind speed, m/s (ETL sonic)
cts5=x(:,16);%relative wind direction (from), deg clockwise rel ship's bow(ETL sonic)
sz6=x(:,17);%latitude, DDMM.MMMM (ETL GPS)
cts6=x(:,18);%longitude, DDDMM.MMMM (ETL GPS)
conv=0.028*10.^6
%%%%%%%%%%%%%%Graphs%%%%%%%%%%%%%%%%%%%%%%%%

figure(1)
plot(,jdy,totcts/conv,'.');
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(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);
%