%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%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(:,); %time of sample (hhmmss)
dir=x(:,3);%true wind direction (from), deg (clockwise rel north)
tsnk=x(:,4);%sea snake temperature, C (0.05 m depth)
tsg=x(:,5);%tsg, tsnk by default
sal=x(:,6);%NA
ta=x(:,7);%air temperature, C (z=15.5 m)
qse=x(:,8);%sea surface specific humidity, g/kg
qa=x(:,9);%air specific humidity, g/kg (z=13 m)
rs=x(:,10);%downward solar flux, W/m^2 (ETL units)
rl=x(:,11);%downward IR flux, W/m^2 (ETL units)
org1=x(:,12);%rainrate, mm/hr (ETL STI optical rain gauge, uncorrected)
ushp=x(:,13);%ship speed, m/s (ETL gps)
head=x(:,14);%ship heading, deg clockwise rel north (ETL gyrocompass)
urel=x(:,15);%relative wind speed, m/s (ETL sonic)
reldir=x(:,16);%relative wind direction (from), deg clockwise rel ship's bow(ETL sonic)
Lat=x(:,17);%latitude, DDMM.MMMM (ETL GPS)
Lon=x(:,18);%longitude, DDDMM.MMMM (ETL GPS)
zts=x(:,19);%depth for bulk flux Ts reference, =0.05 when snake is used
sig_u=x(:,20);%std dev of ship speed, m/s (>.2 indicates maneuver)
Taub=x(:,21);%bulk wind stress along mean wind, N/m^2
Hsb=x(:,22);%bulk sensible heat flux, W/m^2
Hlb=x(:,23);%bulk latent heat flux, W/m^2 (includes Webb et al. correction)
Hrain=x(:,24);%rain heat flux, W/m^2 as per Gosnell et al
Ta_im=x(:,25);% IMET air temp, C (15 m)
Qa_im=x(:,26);% IMET air specific humidity, g/kg (15 m)
U_im=x(:,27);% IMET true wind speed, m/s (15 m)
Dir_im=x(:,28);% IMET true wind direction, deg (15 m)
Psp_im=x(:,29);% IMET solar flux, w/m^2

%%%%%%%%%%%%%%Graphs%%%%%%%%%%%%%%%%%%%%%%%%

figure(1)
plot(jdy,org1);
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);