%
% Diagnostics on Two-Box Model
%
% Global average
  Ts = (Ts1+Ts2)/2;
  Ft = (Ft1+Ft2)/2;
  
% Final time solution
Ta1_end = Ta1(end),
Ta2_end = Ta2(end),
To1_end = To1(end),
To2_end = To2(end),
Ft1_end = Ft1(end),
Ft2_end = Ft2(end),
Ft_end = (Ft1_end+Ft2_end)/2,
Ha_end = Fa(end)*pi*(RADIUS^2)/PW;
Ho_end = Fo(end)*pi*(RADIUS^2)/PW;
Ts1_end = Ts1(end),
Ts2_end = Ts2(end),
Ts_end = (Ts1_end+Ts2_end)/2,
EPSA1_end = EPSA1,
EPSA2_end = EPSA2,
