thermal = getresult('HEAT','thermal');
T = thermal.T;
t = thermal.t;
T_matrix = matrix(1,length(t));
num_grid = length(T)/length(t);

for (i = 1:length(t)){
T_matrix(1,i) = max(T(1:num_grid,i,1));
}

plot(t, T_matrix-T_matrix(1));