%function [  ] = plotCPFdynamics( data )


% --------------------------------
% CPF Dynamics output as follows:
% --------------------------------
% Row 1 --- Time (sec)
% Row 2 --- Current Velocity (m/s)
% Row 3 --- CPF Position (m)
% Row 4 --- CPF Velocity (m/s)
% Row 5 --- CPF Acceleration (m/s/s)
% Row 6 --- CPF Acceleration (mg)
% Row 7 --- Reynolds Number 



figure(1), clf

        
 % Column 1
        subplot(3,2,1)
        plot(data(:,1), data(:,3));
        xlabel('Time [sec]');
        ylabel('CPF Position [m]');
        title('CPF Position');
        
        subplot(3,2,3)
        plot(data(:,1), data(:,4));
        xlabel('Time [sec]');
        ylabel('CPF Velocity [m/s]');
        title('CPF Velocity');
        
        subplot(3,2,5)
        plot(data(:,1), data(:,5));
        xlabel('Time [sec]');
        ylabel('CPF Acceleration [m/s/s]');
        title('CPF Acceleration');
        
        
% Column 2     
        subplot(3,2,2)
        plot(data(:,1), data(:,7));
        xlabel('Time [sec]');
        ylabel('Reynolds Number');
        title('Reynolds Number (need Re > 10^4)');
        
        
        subplot(3,2,4)
        plot(data(:,1), data(:,2));
        xlabel('Time [sec]');
        ylabel('Current Velocity [m/s]');
        title('Current Velocity Profile');
        
        subplot(3,2,6)
        plot(data(:,1), data(:,6));
        xlabel('Time [sec]');
        ylabel('CPF Acceleration [mg]');
        title('CPF Acceleration');
        
        
        
        
        
%end

