 % Assume data loaded...
 
 what_to_match = 33;
 
 % Choose subset of data...
 
 % ok = find((measured(:, 33) > 0.001) & (measured(:, 33) < 0.002) & (measured(:, 34) > 0.002) & (measured(:, 34) < 0.005));
 ok = 1:length(measured(:, what_to_match));
 
 % Compute surrogate
 
 [surrogate, index, weight, residual] = surrogate_LISST(measured(ok, 1:32), measured(ok, what_to_match), 0);
 
 % Plot surrogate vs actual, and time series plots...
 
 figure(1)
 plot(measured(ok, what_to_match), surrogate, '.');
 figure(2)
 
 subplot(2, 1, 1);
 hold off
 plot(measured(ok, 46), measured(ok, what_to_match), '-'); 
 hold on;
 plot(measured(ok, 46), surrogate, 'r-'); 
 hold off
 
 subplot(2, 1, 2);
 hold off
 plot(measured(ok, 46), measured(ok, what_to_match), '-'); 
 hold on;
 plot(measured(ok, 46), surrogate, 'r-'); 
 hold off
 
% ------------------Index of variables...
%    '1'    'L 1.4um'    '2'    'L 1.7um'    '3'     'L 2.0um'    '4'    'L 2.3um'    '5'    'L 2.7um'    '6'    'L 3.2um'
%    '7'    'L 3.8um'    '8'    'L 4.4um'    '9'     'L 5.2um'    '10'    'L 6.1um'
%    '11'    'L 7.2um'   '12'   'L 8.4um'    '13'    'L 10um'     '14'    'L 12um'
%    '15'    'L 14um'    '16'   'L 17um'     '17'    'L 19um'     '18'    'L 23um'
%    '19'    'L 27um'    '20'   'L 32um'     '21'    'L 38um'     '22'    'L 44um'
%    '23'    'L 52um'    '24'   'L 62um'     '25'    'L 73um'     '26'    'L 86um'
%    '27'    'L 101um'   '28'   'L 119um'    '29'    'L 141um'    
%    '30'    'L 166um'   '31'   'L 196um'    '32'    'L 231um'

%    '33'    'Fluor 700nm'  '34'    'OBS 420nm'  '35'    'OBS 700nm' 
%    '36'    'Biolum'       '37'    'Color angle'
%    '38'    'T'            '39'    'S'          '40'    'Density'   
%    '41'    'O2'           '42'    'Nitrate'
%    '43'    'Transmission' '44'    'Pressure'   '45'    'Depth'