%  function [Fo] = calc_Fo(PSIo,To1,To2)
%
    function [Fo] = calc_Fo(PSIo,To1,To2)
 
 %Read parameters
 TwoBox_Param
 
% Oceanic circulation intensity
 [PSIa = calc_Psia(Ts1,Ts2);
 PSIo = PSIfrac * PSIa; 
 
% Heat flux
 Fo = PSIo * CPO * (To1-To2) / (pi*RADIUS^2); %in Wm-2