clear all
close all
format long g

data = csvread('sounding_processed.csv')

data_len = length(data)/2
 
%start_date = datenum('30-Jan-2014 13:05:00')

%for i = 1:data_len
%    date_time(i) = start_date;
%    start_date = addtodate(start_date, 5, 'minute');
%end




figure(1)
plot(data(1:data_len,7),data(1:data_len,6))
set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
set(gca,'XMinorTick','on')
datetick('x','dd/mmm HH:MM','keepticks')
axis([date_time(1) date_time(data_len) min(data(1:data_len,6))-0.5 max(data(1:data_len,6))+0.5])
title('Temperature Profile')
xlabel('Temperature')
ylabel('Height')

% figure(2)
% plot(date_time,data(1:data_len,7))
% set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
% set(gca,'XMinorTick','on')
% datetick('x','dd/mmm HH:MM','keepticks')
% axis([date_time(1) date_time(data_len) min(data(1:data_len,7))-1 max(data(1:data_len,7))+1])
% title('Temperature')
% xlabel('Date and Time')
% ylabel('Temperature (deg C)')
% 
% figure(3)
% plot(date_time,data(1:data_len,8))
% set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
% set(gca,'XMinorTick','on')
% datetick('x','dd/mmm HH:MM','keepticks')
% axis([date_time(1) date_time(data_len) 0 100])
% title('Humidity')
% xlabel('Date and Time')
% ylabel('Humidity (%)')
% 
% figure(4)
% plot(date_time,data(1:data_len,9))
% hold
% plot(date_time,data(1:data_len,10),'r--')
% set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
% set(gca,'XMinorTick','on')
% datetick('x','dd/mmm HH:MM','keepticks')
% legend('Wind Speed','Peak Wind Speed')
% title('Wind Speed')
% xlabel('Date and Time')
% ylabel('Wind Speed (m/s)')
% axis([date_time(1) date_time(data_len) min(data(1:data_len,10))-1 max(data(1:data_len,10))+1])
% 
% figure(5)
% plot(date_time,data(1:data_len,11),'*')
% set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
% set(gca,'XMinorTick','on')
% datetick('x','dd/mmm HH:MM','keepticks')
% axis([date_time(1) date_time(data_len) 0 360])
% title('Wind Direction')
% xlabel('Date and Time')
% ylabel('Wind Direction (deg)')
% 
% figure(6)
% plot(date_time,data(1:data_len,12)+800)
% set(gca,'Xtick',date_time(1):0.125000000931323*2:date_time(data_len))
% set(gca,'XMinorTick','on')
% datetick('x','dd/mmm HH:MM','keepticks')
% axis([date_time(1) date_time(data_len) min(data(1:data_len,12))+799 max(data(1:data_len,12))+801])
% title('Pressure')
% xlabel('Date and Time')
% ylabel('Pressure (hPa)')
