% read and plot 'ONeill_soil_temperatures.dat'

clear all

[header, data] = mhdrload('ONeill_soil_temperatures.dat')

hour = data(:,1)
min = data(:,2)

% plot time series

time = hour + min/60    % decimal hours

temp_time_ser_025 = data(:,3)   % extract data for one depth

figure(1)
clf

plot(time, temp_time_ser_025)

% plot profile

depth = [ 2.5 5 10 20 40 ]  % cm

temp_profile_0435 = data(1, 3:7) % extract data for one time

figure(2)
clf

plot(temp_profile_0435, depth)
set(gca,'YDir','reverse')
