% show LISST_ST

function show_settling(lisst, metadata, txt)

ii = 1:lisst.nsettlingsamples:length(lisst.tm);
jj = ii+lisst.nsettlingsamples-1;
for iExpt = 1:lisst.nsettlingexperiments,
    figure('position',[117         332        1110         590])
    for iSize = 1:length(lisst.sizes),
        subplot(2,4,iSize)
        %semilogx(lisst.sttimes,lisst.conc(ii(iExpt):jj(iExpt),iSize))
        loglog(lisst.sttimes,lisst.conc(ii(iExpt):jj(iExpt),iSize))
        set(gca,'ylim',[1 10e5])
        set(gca,'xlim',[lisst.sttimes(1) lisst.sttimes(end)])
        xlabel('log(t), ms')
        ylabel('ul/l')
        text(0.1,0.9,sprintf('SV=%5.4f cm/s',lisst.svData(iExpt,iSize)), ...
            'units','normalized')
        title(sprintf('%5.2f um',lisst.sizes(iSize)))
        if iSize == 1,
            text(0.5,1.15, ...
                sprintf('Volume concentrations for 8 size classes from file %s, %s, experiment %d', ...
                lisst.procFile, txt, iExpt), 'units','normalized')
        end
    end
end