ssbin2mlab(filename, desired_variable)

function [data, plot_title] = ssbin2mlab(filename, desired_variable)

% [data, plot_title] = ssbin2mlab(filename, desired_variable) reads a 
% Binary Smartspice RAW file 'filename' and extract the variable(s) of
% interest from 'desired_variable' using the file's own variable index 
% format and place the data in column into the variable 'data' and also 
% return a title into 'plot_title'.
%
% Input argument 'desired_variable' can be omitted if all variables are
% desired. If 'desired_variable' is set to [] (the empty matrix), the
% function will only return the index of all the variables from that file
% and no data is read.
%
% Simon Damphousse, July 09.

if (nargin == 1);
    desired_variable = 1;
end

%open the file
fid = fopen(filename,'r');
if fid == -1
	error('file not opened');
end
% display the index
if isempty(desired_variable)
    go=0;
    impr=0;
    disp(' ');
    disp('The variables are:')
    disp(' ');
    while (go == 0)    
        line = fgetl(fid);
        w = strfind(line,'Variables:');
        if w == 1,
            impr=1;
            line = line(11:max(size(line)));
        end
        w = strfind(line,'Binary:');
        if w == 1,
            impr=0;
            go=1;
        end
        if impr == 1
            
            disp(line);
        end
    end
    data=0;
    plot_title = 'noplot';

% Loading data
else
    go=0;
    while (go == 0)
        line = fgetl(fid);
        %find the name of the circuit
        w = strfind(line,'Input deck file name:');
        if ~isempty(w)
            id1 = strfind(line,'\');
            id1e = id1(max(size(id1)));
            id2 = strfind(line,'.cir');
            circuit = line(id1e+1:id2-1);
        end
        %find the plot name.
        w = strfind(line,'Plotname:');
        if ~isempty(w)
            id1 = strfind(line,':');
            id2 = max(size(line));
            name = line(id1+1:id2);
        end
        %find the number of variables
        w = strfind(line,'No. Variables:');
        if ~isempty(w);
            id1 = strfind(line,':');
            id2 = max(size(line));
            nv = str2num(line(id1+1:id2));
        end
        %find the number of points
        w = strfind(line,'No. Points:');
        if ~isempty(w);
            id1 = strfind(line,':');
            id2 = max(size(line));
            np = str2num(line(id1+1:id2));
        end
        %Find if complex or real data
        w = strfind(line,'Flags:');
        if ~isempty(w);
            id1 = strfind(line,':');
            id2 = max(size(line));
            cplx = strfind(line(id1+1:id2),'complex');
            if ~isempty(cplx)
                cplx=2;
                crs='the data type is complex';
            else
                cplx=1;
                crs='the data type is real';
            end
        end
        
        go = strcmp(line,'Binary:');
    end

    if (nargin == 1);
        desired_variable = [0:nv-1];
    end

    disp(['number of variables selected/total: ',num2str(max(size(desired_variable))), '/', num2str(nv)]);
    disp(['number of points: ',num2str(np)]);
    disp(crs);

    plot_title = [circuit,': ', name]; 
    % reading all data on the file
    donnee = double(fread(fid, nv*np*cplx, 'double'));
    dim = max(size(donnee));
    data = [];
    %Sorting data and create matrix
    if cplx == 1,
        for i=(desired_variable)+1,
            data = [data , donnee(i:nv:dim)];
        end
    end
    if cplx == 2,
        for i=(desired_variable)+1,
            data = [data , donnee((2*i-1):nv*cplx:dim)+j*donnee((2*i):nv*cplx:dim)];
        end
    end
    
end
fclose(fid);