function unitCompare(culture)

times = {'_13DIV_19to20hr.mat', '_14DIV_36to37hr.mat', '_15DIV_64to65hr.mat', '_16DIV_108to109hr.mat'};

fileName = ['C:\Users\USER\Documents\mingfaifong\haruki17898\matfiles\' culture times{1}];
unitData = load(fileName, 'unit');
unitNum = max(unitData.unit);

fh = figure('units','normalized','outerposition',[0 0 .5 1]);

for unitID = 1:unitNum
    figure(fh)
    yLL = [];
    yUL = [];
    for i = 1:length(times)
        fileName = ['C:\Users\USER\Documents\mingfaifong\haruki17898\matfiles\' culture times{i}];
        load(fileName);
        subplot(length(times), 1, i);
        plot(waveform(:,unit==unitID));
        %ylims = ylim;
%         yLL = [yLL ylims(1)];
%         yUL = [yUL ylims(2)];
        yLL = [yLL min(min(waveform(:,unit==unitID)))];
        yUL = [yUL max(max(waveform(:,unit==unitID)))];
        title(['Unit Number ' num2str(unitID)]);
    end
    yMin = min(yLL);
    yMax = max(yUL);
    if size([yMin yMax])==[1 2]
        for i = 1:length(times)
            subplot(length(times), 1, i);
            ylim([yMin yMax]);
        end
    end
        print('-dpng',['C:\Users\USER\Documents\mingfaifong\haruki17898\unitcheck\' culture '_unit' num2str(unitID) '.png']);
end 

end