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