%plotMeteorStatsNY.m
%
%Program to create a plot comparing the meteorological statistics from the
%buoys near the NY Bight study area
%
%Soupy Alexander, 2/7/2

meteorStatsNY

figure(1)
subplot(2,3,1)
hold on
for index = 1:5;
    plot(index,mean_windu(1,index),'.')
end
    plot(6,mean_windu(1,6),'*')


subplot(2,3,4)
hold on
for index = 1:5;
    plot(index,std_windu(1,index),'.')
end
    plot(6,std_windu(1,6),'*')
axis([0 7 0 10])


subplot(2,3,2)
hold on
for index = 1:5;
    plot(index,mean_windv(1,index),'.')
end
    plot(6,mean_windv(1,6),'*')


subplot(2,3,5)
hold on
for index = 1:5;
    plot(index,std_windv(1,index),'.')
end
    plot(6,std_windv(1,6),'*')


subplot(2,3,3)
hold on
for index = 1:5;
    plot(index,mean_windsp(1,index),'.')
end
    plot(6,mean_windsp(1,6),'*')

subplot(2,3,6)
hold on
for index = 1:5;
    plot(index,std_windsp(1,index),'.')
end
    plot(6,std_windsp(1,6),'*')


figure(2)
subplot(2,3,1)
hold on
for index = 1:5;
    plot(index,mean_taux(1,index),'.')
end
    plot(6,mean_taux(1,6),'*')


subplot(2,3,4)
hold on
for index = 1:5;
    plot(index,std_taux(1,index),'.')
end
    plot(6,std_taux(1,6),'*')
axis([0 7 0 10])


subplot(2,3,2)
hold on
for index = 1:5;
    plot(index,mean_tauy(1,index),'.')
end
    plot(6,mean_tauy(1,6),'*')


subplot(2,3,5)
hold on
for index = 1:5;
    plot(index,std_tauy(1,index),'.')
end
    plot(6,std_tauy(1,6),'*')


subplot(2,3,3)
hold on
for index = 1:5;
    plot(index,mean_tauxy(1,index),'.')
end
    plot(6,mean_tauxy(1,6),'*')

subplot(2,3,6)
hold on
for index = 1:5;
    plot(index,std_tauxy(1,index),'.')
end
    plot(6,std_tauxy(1,6),'*')


figure(3)
subplot(2,2,1)
hold on
for index = 1:5;
    plot(index,mean_airtemp(1,index),'.')
end
plot(6,mean_airtemp(1,6),'*')
    
subplot(2,2,2)
hold on
for index = 1:5;
    plot(index,mean_sstemp(1,index),'.')
end
plot(6,mean_sstemp(1,6),'*')

subplot(2,2,3)
hold on
for index = 1:5;
    plot(index,std_airtemp(1,index),'.')
end
plot(6,std_airtemp(1,6),'*')

subplot(2,2,4)
hold on
for index = 1:5;
    plot(index,std_sstemp(1,index),'.')
end
plot(6,std_sstemp(1,6),'*')

meteorStatsNY

figure(1)
subplot(2,3,1)
hold on
for index = 1:5;
    plot(index,mean_windu(1,index),'r.')
end
    plot(6,mean_windu(1,6),'r*')
axis([0 7 -6 6])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[-6:2:6])
ylabel('WU422 (cm/s)')
title('Wind U')

subplot(2,3,4)
hold on
for index = 1:5;
    plot(index,std_windu(1,index),'r.')
end
    plot(6,std_windu(1,6),'r*')
axis([0 7 0 10])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:10])
ylabel('SD WU422 (cm/s)')
title('SD Wind U')

subplot(2,3,2)
hold on
for index = 1:5;
    plot(index,mean_windv(1,index),'r.')
end
    plot(6,mean_windv(1,6),'r*')
axis([0 7 -6 6])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[-6:2:6])
ylabel('WV422 (cm/s)')
title('Wind V')

subplot(2,3,5)
hold on
for index = 1:5;
    plot(index,std_windv(1,index),'r.')
end
    plot(6,std_windv(1,6),'r*')
axis([0 7 0 10])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:10])
ylabel('SD WV422 (cm/s)')
title('SD Wind V')

subplot(2,3,3)
hold on
for index = 1:5;
    plot(index,mean_windsp(1,index),'r.')
end
    plot(6,mean_windsp(1,6),'r*')
axis([0 7 0 12])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:12])
ylabel('WS400 (cm/s)')
title('Wind Speed')

subplot(2,3,6)
hold on
for index = 1:5;
    plot(index,std_windsp(1,index),'r.')
end
    plot(6,std_windsp(1,6),'r*')
axis([0 7 0 10])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:10])
ylabel('SD WS400 (cm/s)')
title('SD Wind Speed')

figure(2)
subplot(2,3,1)
hold on
for index = 1:5;
    plot(index,mean_taux(1,index),'r.')
end
    plot(6,mean_taux(1,6),'r*')
axis([0 7 -1.4 1.4])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[-1.4:.2:1.4])
ylabel('TX440 (dynes/cm^2)')
title('Zonal Wind Stress')

subplot(2,3,4)
hold on
for index = 1:5;
    plot(index,std_taux(1,index),'r.')
end
    plot(6,std_taux(1,6),'r*')
axis([0 7 0 2.4])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:.2:2.4])
ylabel('SD TX440 (dynes/cm^2)')
title('SD Zonal Wind Stress')

subplot(2,3,2)
hold on
for index = 1:5;
    plot(index,mean_tauy(1,index),'r.')
end
    plot(6,mean_tauy(1,6),'r*')
axis([0 7 -1.4 1.4])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[-1.4:.2:1.4])
ylabel('TY441 (dynes/cm^2)')
title('Meridional Wind Stress')

subplot(2,3,5)
hold on
for index = 1:5;
    plot(index,std_tauy(1,index),'r.')
end
    plot(6,std_tauy(1,6),'r*')
axis([0 7 0 2.4])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:.2:2.4])
ylabel('SD TY441 (dynes/cm^2)')
title('SD Meridional Wind Stress')

subplot(2,3,3)
hold on
for index = 1:5;
    plot(index,mean_tauxy(1,index),'r.')
end
    plot(6,mean_tauxy(1,6),'r*')
axis([0 7 0 2.42])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:.2:2.4])
ylabel('Txy448 (dynes/cm^2)')
title('Total Wind Stress')

subplot(2,3,6)
hold on
for index = 1:5;
    plot(index,std_tauxy(1,index),'r.')
end
    plot(6,std_tauxy(1,6),'r*')
axis([0 7 0 2.42])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:.2:2.4])
ylabel('SD Txy448 (dynes/cm^2)')
title('SD Total Wind Stress')

figure(3)
subplot(2,2,1)
hold on
for index = 1:5;
    plot(index,mean_airtemp(1,index),'r.')
end
plot(6,mean_airtemp(1,6),'r*')
axis([0 7 0 12])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:12])
ylabel('AT_21 (Deg. C)')
title('Mean Air Temperature')
    
subplot(2,2,2)
hold on
for index = 1:5;
    plot(index,mean_sstemp(1,index),'r.')
end
plot(6,mean_sstemp(1,6),'r*')
axis([0 7 0 12])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:12])
ylabel('T_25 (Deg. C)')
title('Mean Sea Surface Temperature')

subplot(2,2,3)
hold on
for index = 1:5;
    plot(index,std_airtemp(1,index),'r.')
end
plot(6,std_airtemp(1,6),'r*')
axis([0 7 0 6])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:6])
ylabel('SD AT_21 (Deg. C)')
title('SD Air Temperature')

subplot(2,2,4)
hold on
for index = 1:5;
    plot(index,std_sstemp(1,index),'r.')
end
plot(6,std_sstemp(1,6),'r*')
axis([0 7 0 6])
set(gca,'xtick',[1:6]);
set(gca,'xticklabel',['Dec'; 'Jan'; 'Feb'; 'Mar'; 'Apr'; 'All';]);
set(gca,'ytick',[0:2:6])
ylabel('SD T_25 (Deg. C)')
title('SD Sea Surface Temperature')