%  function [Fa] = calc_Fa(Ta1,Ta2,To1,To2)
%
    function [Fa] = calc_Fa(Ta1,Ta2,Ts1,Ts2)
 
 %Read parameters
 TwoBox_Param
 
% Atm. circulation intensity
 PSIa = calc_Psia(Ts1,Ts2); %in kg s-1
  
% Heat flux
 Fa = PSIa * CAheat * (Ta1-Ta2) / (pi*RADIUS^2); %in Wm-2
 
