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