function spe = outputTMMfS(Qext, Qsca, Qabs, WavlenRange, varMode, ...
    outputFile, outputData, outputPlot)
% outputTMMfS(Qext, Qsca, Qabs, WavlenRange, varMode, ...
%             outputFile, outputData, outputPlot)

%% Output
spe = WavlenRange';

outputN = height_y(outputData);
for j = 1:outputN
    Q = outputData(j,:);
    Qvalue = eval(Q);
    spe = [spe Qvalue.'];
end
save(outputFile, 'spe', '-ascii')

if isequal(outputPlot ,'yes')
    figure('Name','Q');
    for j = 1:outputN
        Q = outputData(j,:);
        Qvalue = eval(Q);
    
        subplot(outputN,1,j)
        if varMode=='w'
            plot(WavlenRange, Qvalue)
            xlabel('\lambda')
            ylabel(Q)
        else
            plot(fliplr(1./WavlenRange), fliplr(Qvalue))
            xlabel('\lambda^{-1}')
            ylabel(Q)
        end
    end
end