function torque_tc = interp_engine_w_torque_converter(omega_tc)

% [TORQUE_TC] = ENGINE_W_TORQUE_CONVERTER (OMEGA_TC)
%   computes the torque converter torque, TORQUE_TC [Nm], as a function of
%   the  rotational velocity of the torque converter shaft, OMEGA_TC [rad/s].
%   This model is a simplified version of ENGINE_W_TORQUE_CONVERTER;  it
%   avoids the root-finding problem (which is very computationally
%   expensive) through interpolation.  The interpolation table was
%   determined by solving ENGINE_W_TORQUE_CONVERTER for 200 equally spaced
%   values for OMEGA_TC between 0 and 700.
%
%   INPUT:
%          omega_tc: speed at the output of the torque converter in rad/s
%   OUTPUT:
%          torque_tc: 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
%%

% hardcode the interpolation data here.
omega_tc_data = linspace(0,700,200);
torque_tc_data =  [ 552.7146
  546.4182
  540.2245
  534.1356
  528.1529
  522.2776
  516.5105
  510.8828
  505.5013
  500.2961
  495.1778
  490.0581
  484.8500
  479.4678
  473.9188
  468.2835
  462.5974
  456.8925
  451.1979
  445.5400
  439.9230
  434.2889
  428.6575
  423.0582
  417.5203
  412.0725
  406.7424
  401.5384
  396.4499
  391.4686
  386.5858
  381.7924
  377.0790
  372.4406
  367.9176
  363.5083
  359.1963
  354.9721
  350.8386
  346.8006
  342.8451
  338.9511
  335.1351
  331.4125
  327.7971
  324.3010
  320.9363
  317.7142
  314.6146
  311.6148
  308.6945
  305.8354
  303.0212
  300.2422
  297.5137
  294.8461
  292.2474
  289.7241
  287.2817
  284.9245
  282.6825
  280.5556
  278.5070
  276.5047
  274.5219
  272.5354
  270.5260
  268.5026
  266.4863
  264.4837
  262.4990
  260.5348
  258.5929
  256.6631
  254.5136
  252.4207
  250.8425
  249.7904
  248.9338
  248.5016
  248.5631
  248.6832
  248.7950
  248.8975
  248.9900
  249.0720
  249.1431
  249.2030
  249.2513
  249.2879
  249.3126
  249.3255
  249.3264
  249.3153
  249.2922
  249.2571
  249.2099
  249.1504
  249.0786
  248.9942
  248.8970
  248.7869
  248.6636
  248.5268
  248.3764
  248.2121
  248.0337
  247.8409
  247.6337
  247.4116
  247.1746
  246.9224
  246.6547
  246.3712
  246.0718
  245.7560
  245.4236
  245.0741
  244.7073
  244.3234
  243.9231
  243.5050
  243.0677
  242.6095
  242.1292
  241.6254
  241.0970
  240.5429
  239.9620
  239.3531
  238.7150
  238.0447
  237.3411
  236.6037
  235.8314
  235.0229
  234.1767
  233.2913
  232.3647
  231.3951
  230.3804
  229.3182
  228.2062
  227.0419
  225.8227
  224.5457
  223.2079
  221.8064
  220.3378
  218.7988
  217.1859
  215.4954
  213.7171
  211.8444
  209.8756
  207.8079
  205.6384
  203.3645
  200.9849
  198.5008
  195.9212
  193.2666
  190.5139
  187.6600
  184.7069
  181.6568
  178.5099
  175.2620
  171.9009
  168.4061
  164.7514
  160.9330
  156.9703
  145.8915
  114.0246
   88.1667
   65.4435
   48.7037
    0.5147
   -0.3144
   -8.1803
  -25.3398
  -49.6650
  -79.0279
 -111.3005
 -144.3551
 -176.0635
 -204.2979
 -226.9304
 -242.7132
 -255.2494
 -266.4491
 -277.2727
 -287.5812
 -297.4536
 -306.9754
 -316.2320
 -325.3089
 -334.2916
 -343.2565];

% use the interpolation table to determine the torque 
% spline interpolation is used to obtain a smooth response 
torque_tc = interp1(omega_tc_data, torque_tc_data, omega_tc, 'spline');