function [torque_engine, torque_load] = torque_converter(omega_engine, omega_load)

% [TORQUE_ENGINE, TORQUE_LOAD] = TORQUE_CONVERTER (OMEGA_ENGINE, OMEGA_LOAD)
%   computes the engine torque, TORQUE_ENGINE [Nm] and the load torque,
%   TORQUE_LOAD [Nm] as a function of engine and load speeds, OMEGA_ENGINE
%   and OMEGA_LOAD, both in rad/s. 
%
%   INPUT:
%          omega_engine: speed of the engine in rad/s
%          omega_load: speed at the output of the torque converter in rad/s
%
%   OUTPUT:
%          torque_engine: torque at produced by the engine in Nm
%          torque_load: torque at the output of the torque conververter in Nm

%
%   Assumptions:
%     - Although qualitatively correct, this function is based entirely
%       on fictitious data.

%%
%% ME 2016
%% Fall 2007
%%
%% AUTHOR:  Chris Paredis
%%

% define the vectors with (experimental) data
sr_axis=[0.000,0.098,0.196,0.289,0.379,0.465,0.500,0.543,0.612,0.676,...
        0.730,0.777,0.815,0.834,0.848,0.889,0.917,0.935,0.946,0.955,0.958,0.965,...
        0.974,0.987,0.995,0.996,0.997,1.001,1.055,1.070,1.111,1.176,1.294,1.544,...
        2.095,3.816]; 
tr_axis=[2.306,2.144,1.998,1.842,1.688,1.544,1.488,1.422,1.318,1.228,...
        1.152,1.088,1.028,1.000,0.990,0.990,0.990,0.990,0.990,0.990,0.990,0.990,...
        0.990,0.990,0.990,0.990,0.990,1.100,1.100,1.100,1.100,1.100,1.100,1.100,...
        1.100,1.100];
C_axis=[(-0.00451386747055463),(-0.00449754738459139),(-0.00443314419611169),...
        (-0.00431351074122866),(-0.00417914948909407),(-0.00397777464756128),(-0.00388903734640203),...
        (-0.0037780496964034),(-0.00356220648753515),(-0.00331230679772083),(-0.00306022025408675),...
        (-0.00279288677201309),(-0.00253368936382028),(-0.00242277084148496),(-0.00231096588569689),...
        (-0.00186751985415466),(-0.00152749662232769),(-0.00125077653551252),(-0.00103873941176525),...
        (-0.000704569889647544),(-0.000567878375277169),(-0.000432115260711337),...
        (-0.000293351138282709),(-0.000150655253559483),(-7.65218326576329E-005),...
        (-4.94543833941436E-006),(-1.00927313049273E-007),1.00927313049273E-007,...
        0.00054036792233659,0.000616754617750885,0.000779061152125976,...
        0.00100124713347082,0.00129019431543173,0.00161342904709868,...
        0.00197186033473214,0.00237397126625058];

% compute the speed ratio from which to compute the torque ratio and C
% value (some sort of impedance value)
speed_ratio = omega_load ./ max(omega_engine, eps);
torque_ratio = interp1(sr_axis, tr_axis, speed_ratio, 'pchip');
C = interp1(sr_axis, C_axis, speed_ratio, 'pchip');

% compute the actual torque values based on the turbine characteristic of
% the torque converter
torque_engine = -C*omega_engine^2;
torque_load = torque_engine .* torque_ratio;

return;