function [] = plotSWIFT(SWIFT);

% wind and wave plot
figure(1), clf, n = 4; 

ax(1) = subplot(n,1,1);
plot( [SWIFT.time],[SWIFT.windspd],'bx') 
datetick
ylabel([ 'windppd [m/s]'])
set(gca,'Ylim',[0 20])

ax(2) = subplot(n,1,2);
plot( [SWIFT.time],[SWIFT.sigwaveheight],'g+') 
datetick
ylabel([ 'waveheight [m]'])
set(gca,'Ylim',[0 7])


ax(3) = subplot(n,1,3);
plot( [SWIFT.time],[SWIFT.peakwaveperiod],'g+') 
datetick
ylabel([ 'waveperiod [s]'])
set(gca,'Ylim',[0 12])

ax(4) = subplot(n,1,4);
plot([SWIFT.time],[SWIFT.winddirT],'bx'), hold on, 
plot([SWIFT.time],[SWIFT.peakwavedirT],'g+'), hold on
datetick
ylabel('directions [^\circ T]')
set(gca,'Ylim',[0 360])
set(gca,'YTick',[0 180 360])


linkaxes(ax,'x')
set(gca,'XLim',[min([SWIFT.time]) max([SWIFT.time])] )
print -dpng windandwaves.png

% temperature plot
figure(2), clf, n = 3; 
 
tax(1) = subplot(n,1,1); 
plot( [SWIFT.time],[SWIFT.airtemp],'g+') 
datetick
ylabel(['airtemp [C]'])
set(gca,'Ylim',[-10 15])


tax(2) = subplot(n,1,2); 
plot( [SWIFT.time],[SWIFT.watertemp],'g+') 
datetick
ylabel(['watertemp [C]'])
set(gca,'Ylim',[-2 18])


tax(3) = subplot(n,1,3); 
plot( [SWIFT.time],[SWIFT.salinity],'g+') 
datetick
ylabel(['salinity [PSU]'])
set(gca,'Ylim',[20 35])


linkaxes(tax,'x')
set(gca,'XLim',[min([SWIFT.time]) max([SWIFT.time])])
print -dpng tempandsalinity.png
 
% turbulence plot
figure(3), clf
for ai = 1:length(SWIFT), 

        plot(SWIFT(ai).uplooking.tkedissipationrate,SWIFT(ai).uplooking.z), hold on
        set(gca,'YDir','reverse')
        ylabel('z [m]')
        xlabel('\epsilon [W/kg]')
        set(gca,'Ylim',[0 0.7])
        set(gca,'Xlim',[0 1e-2])
        
end
print -dpng turbulence.png
        
% wave spectra plot
figure(4), clf
for ai = 1:length(SWIFT), 

        loglog(SWIFT(ai).wavespectra.freq,SWIFT(ai).wavespectra.energy), hold on
        xlabel('freq [Hz]')
        ylabel('Energy [m^2/Hz')
end
print -dpng wavespectra.png
 
% downlooking velocity profile plot
figure(5), clf
for ai = 1:length(SWIFT), 
        
        plot(SWIFT(ai).downlooking.velocityprofile,SWIFT(ai).downlooking.z), hold on
        set(gca,'YDir','reverse')
        ylabel('z [m]')
        xlabel('velocity [m/s]')
end
print -dpng velocityprofiles.png

% drift plot
figure(6), %clf
%quiver(lon,lat,dlondt,dlatdt,1), hold on
quiver([SWIFT.lon],[SWIFT.lat],[SWIFT.driftspd].*sind([SWIFT.driftdirT]),[SWIFT.driftspd].*cosd([SWIFT.driftdirT]),1,'r'), hold on
xlabel('longitude'), ylabel('latitude')
axlims = axis;
%quiver(axlims(1) +(axlims(2)-axlims(1))./10, axlims(3)+(axlims(4)-axlims(3))./10, .1, 0,0 );
%text(axlims(1) +(axlims(2)-axlims(1))./9, axlims(3)+(axlims(4)-axlims(3))./8,'0.1 m/s')
plot([SWIFT.lon],[SWIFT.lat],'bo','markersize',2), hold on
%plot(lon(length(lon)),lat(length(lon)),'r.','markersize',20), hold on
print -dpng drift.png


