%  function [Fo] = calc_Fo(PSIo,To1,To2)
%
    function [Fo] = calc_Fo(PSIo,To1,To2)
 
 %Read parameters
 TwoBox_Param
 
% Heat flux
 Fo = PSIo * CPO * (To1-To2) / (pi*RADIUS^2); %in Wm-2
 
