%  function [PSIa,PSIo] = calc_Psi(Ts1,Ts2)
%
    function [Psia,Psio] = calc_Psi(Ts1,Ts2)
 
%Read parameters
  TwoBox_Param
   
%Dynamical parameterization
  Psia = KEFF * (Ts1-Ts2); %in kg/s
  Psio = Psia * PSIfrac; %in kg/s
