%  function [Fa] = calc_Fa(PSIa,Ta1,Ta2)
%
    function [Fa] = calc_Fa(PSIa,Ta1,Ta2)
 
 %Read parameters
 TwoBox_Param
 
% Heat flux
 Fa = PSIa * CAheat * (Ta1-Ta2) / (pi*RADIUS^2); %in Wm-2
 
