%PlotStatsNY
%
%Program to plot the statistics calculated for the NY Bight in Stats_Curr_Temp.
%Modifed from Ben Guitterez's "plotstats.m".  This program treats the 5 months
%of the NY Bight data as months 1-5, then relabels them to Dec.-Apr.
%
%Soupy Alexander, 12/19/2001

if hit_depth_type == 's';
    ymin = -150;
    ymax = 150;
else
    ymin = -100;
    ymax = 100;
end

ymin_temp = 0;
ymax_temp = 20;

monthvector = [12 1 2 3 4];

subplot(2,2,1)
plot(value2Index(monthvector,ut),umean,'b',value2Index(monthvector,ut),umean,'bo');
hold on
plot(value2Index(monthvector,ut),ustd,'g',value2Index(monthvector,ut),ustd,'go');
plot(value2Index(monthvector,ut),umax,'r',value2Index(monthvector,ut),umax,'ro');
plot(value2Index(monthvector,ut),umin,'m',value2Index(monthvector,ut),umin,'mo');
plot(value2Index(monthvector,yt),uymean,'b.','MarkerSize',20)
plot(value2Index(monthvector,yt),uystd,'g.','MarkerSize',20)
axis([1 5 ymin ymax])
set(gca,'xtick',[1:5])
set(gca,'xticklabel',monthvector);
title('Eastward Velocity Statistics')
ylabel('cm/s')



subplot(2,2,2)
plot(value2Index(monthvector,vt),vmean,'b',value2Index(monthvector,vt),vmean,'bo');
hold on
plot(value2Index(monthvector,vt),vstd,'g',value2Index(monthvector,vt),vstd,'go');
plot(value2Index(monthvector,vt),vmax,'r',value2Index(monthvector,vt),vmax,'ro');
plot(value2Index(monthvector,vt),vmin,'m',value2Index(monthvector,vt),vmin,'mo');
plot(value2Index(monthvector,yt),vymean,'b.','MarkerSize',20)
plot(value2Index(monthvector,yt),vystd,'g.','MarkerSize',20)
axis([1 5 ymin ymax])
set(gca,'xtick',[1:5])
set(gca,'xticklabel',monthvector);
title('Northward Velocity Statistics')
ylabel('cm/s')



subplot(2,2,3)
plot(value2Index(monthvector,cst),csmean,'b',value2Index(monthvector,cst),csmean,'bo');
hold on
plot(value2Index(monthvector,cst),csstd,'g',value2Index(monthvector,cst),csstd,'go')
plot(value2Index(monthvector,cst),csmax,'r',value2Index(monthvector,cst),csmax,'ro')
plot(value2Index(monthvector,cst),csmin,'m',value2Index(monthvector,cst),csmin,'mo')
plot(value2Index(monthvector,yt),csymean,'b.','MarkerSize',20)
plot(value2Index(monthvector,yt),csystd,'g.','MarkerSize',20)
axis([1 5 0 ymax])
set(gca,'xtick',[1:5])
set(gca,'xticklabel',monthvector);
title('Current Speed Statistics')
ylabel('cm/s')




subplot(2,2,4)
plot(value2Index(monthvector,tt),tmean,'b',value2Index(monthvector,tt),tmean,'bo');
hold on
plot(value2Index(monthvector,tt),tstd,'g',value2Index(monthvector,tt),tstd,'go')
plot(value2Index(monthvector,tt),tmax,'r',value2Index(monthvector,tt),tmax,'ro')
plot(value2Index(monthvector,tt),tmin,'m',value2Index(monthvector,tt),tmin,'mo')
plot(value2Index(monthvector,tyt),tymean,'b.','MarkerSize',20)
plot(value2Index(monthvector,tyt),tystd,'g.','MarkerSize',20)
axis([1 5 ymin_temp ymax_temp])
set(gca,'xtick',[1:5])
set(gca,'xticklabel',monthvector);
title('Temperature Statistics')
ylabel('deg. C')