function [] = plot_polarity_latlon(llatlon, siteData, quietwin, IARwin, day_or_frames, R, C, s)
% Function that plots the lat and lon position of a lightning strike
% matched to a B peak, differentiating between B peaks that have
% positive and negative polarity

sf3 = figure('Position', [0 0 s(3) s(4)]);%('Position', [0 0 s(3) s(4)]); % New figure

% Station location vecs
STN_Lat=siteData.geoLat; STN_Long=siteData.geoLon-360;

% Declare distance intervals in km
dl = 500:50:1000;
dl = dl';
% dl1 = 920;
% dl2 = 930;
% dl3 = 940;
% dl4 = 950;
% dl5 = 960;

% Declare empty variables for each bar set quadrant
Q1p = zeros(size(dl,1),1); Q1n = zeros(size(dl,1),1);
Q2p = zeros(size(dl,1),1); Q2n = zeros(size(dl,1),1);
Q3p = zeros(size(dl,1),1); Q3n = zeros(size(dl,1),1);
Q4p = zeros(size(dl,1),1); Q4n = zeros(size(dl,1),1);
% Q1ap = 0; Q1an = 0; Q1bp = 0; Q1bn = 0; Q1cp = 0; Q1cn = 0; Q1dp = 0; Q1dn = 0;
% Q2ap = 0; Q2an = 0; Q2bp = 0; Q2bn = 0; Q2cp = 0; Q2cn = 0; Q2dp = 0; Q2dn = 0;
% Q3ap = 0; Q3an = 0; Q3bp = 0; Q3bn = 0; Q3cp = 0; Q3cn = 0; Q3dp = 0; Q3dn = 0;
% Q4ap = 0; Q4an = 0; Q4bp = 0; Q4bn = 0; Q4cp = 0; Q4cn = 0; Q4dp = 0; Q4dn = 0;
ax1 = 1;

%c = categorical({'Q1-NE','Q2-NW','Q3-SW','Q4-SE','Total'});
c = categorical({'Q1-NE-Neg','Q1-NE-Pos','Q2-NW-Neg','Q2-NW-Pos','Q3-SW-Neg','Q3-SW-Pos','Q4-SE-Neg','Q4-SE-Pos'});


% Calculate i value for quietwin and Iarwin
IAR_per_int = 2*IARwin(1) : 2*IARwin(2);

if strcmp(day_or_frames,'day')
    per_int = 1:length(llatlon);
elseif strcmp(day_or_frames,'frames')
    per_int = [round(2*mean(quietwin/(100*3600))) IAR_per_int(2:end)];
end

for i = per_int
    for j = 1:size(llatlon{i},1)
        for k = 1:size(dl,1)-1
            if llatlon{i}(j,1) > STN_Lat && llatlon{i}(j,2) > STN_Long
                if llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) > 0
                    Q1p(k,1) = Q1p(k,1) + 1;
                elseif llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) < 0
                    Q1n(k,1) = Q1n(k,1) + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) > 0
                    %                 Q1bp = Q1bp + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) < 0
                    %                 Q1bn = Q1bn + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) > 0
                    %                 Q1cp = Q1cp + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) < 0
                    %                 Q1cn = Q1cn + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) > 0
                    %                 Q1dp = Q1dp + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) < 0
                    %                 Q1dn = Q1dn + 1;
                end
                
            elseif llatlon{i}(j,1) > STN_Lat && llatlon{i}(j,2) < STN_Long
                if llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) > 0
                    Q2p(k,1) = Q2p(k,1) + 1;
                elseif llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) < 0
                    Q2n(k,1) = Q2n(k,1) + 1;
                    %             if llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) > 0
                    %                 Q2ap = Q2ap + 1;
                    %             elseif llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) < 0
                    %                 Q2an = Q2an + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) > 0
                    %                 Q2bp = Q2bp + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) < 0
                    %                 Q2bn = Q2bn + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) > 0
                    %                 Q2cp = Q2cp + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) < 0
                    %                 Q2cn = Q2cn + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) > 0
                    %                 Q2dp = Q2dp + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) < 0
                    %                 Q2dn = Q2dn + 1;
                end
                
            elseif llatlon{i}(j,1) < STN_Lat && llatlon{i}(j,2) < STN_Long
                if llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) > 0
                    Q3p(k,1) = Q3p(k,1) + 1;
                elseif llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) < 0
                    Q3n(k,1) = Q3n(k,1) + 1;
                    %             if llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) > 0
                    %                 Q3ap = Q3ap + 1;
                    %             elseif llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) < 0
                    %                 Q3an = Q3an + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) > 0
                    %                 Q3bp = Q3bp + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) < 0
                    %                 Q3bn = Q3bn + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) > 0
                    %                 Q3cp = Q3cp + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) < 0
                    %                 Q3cn = Q3cn + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) > 0
                    %                 Q3dp = Q3dp + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) < 0
                    %                 Q3dn = Q3dn + 1;
                end
                
            elseif llatlon{i}(j,1) < STN_Lat && llatlon{i}(j,2) > STN_Long
                if llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) > 0
                    Q4p(k,1) = Q4p(k,1) + 1;
                elseif llatlon{i}(j,3) >= dl(k) && llatlon{i}(j,3) < dl(k+1) && llatlon{i}(j,4) < 0
                    Q4n(k,1) = Q4n(k,1) + 1;
                    %             if llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) > 0
                    %                 Q4ap = Q4ap + 1;
                    %             elseif llatlon{i}(j,3) >= dl1 && llatlon{i}(j,3) < dl2 && llatlon{i}(j,4) < 0
                    %                 Q4an = Q4an + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) > 0
                    %                 Q4bp = Q4bp + 1;
                    %             elseif llatlon{i}(j,3) >= dl2 && llatlon{i}(j,3) < dl3 && llatlon{i}(j,4) < 0
                    %                 Q4bn = Q4bn + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) > 0
                    %                 Q4cp = Q4cp + 1;
                    %             elseif llatlon{i}(j,3) >= dl3 && llatlon{i}(j,3) < dl4 && llatlon{i}(j,4) < 0
                    %                 Q4cn = Q4cn + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) > 0
                    %                 Q4dp = Q4dp + 1;
                    %             elseif llatlon{i}(j,3) >= dl4 && llatlon{i}(j,3) < dl5 && llatlon{i}(j,4) < 0
                    %                 Q4dn = Q4dn + 1;
                end
            end
        end
    end
    
    % Bar graph for same quiet and IAR windows
    if strcmp(day_or_frames,'frames')
        %Qpnbins = [Q1n Q1p; Q2n Q2p; Q3n Q3p; Q4n Q4p; Q1n+Q2n+Q3n+Q4n Q1p+Q2p+Q3p+Q4p];
        %Qpnbins = [Q1an Q1bn Q1cn Q1dn; Q1ap Q1bp Q1cp Q1dp; Q2an Q2bn Q2cn Q2dn; Q2ap Q2bp Q2cp Q2dp; Q3an Q3bn Q3cn Q3dn; Q3ap Q3bp Q3cp Q3dp; Q4an Q4bn Q4cn Q4dn; Q4ap Q4bp Q4cp Q4dp];
        Qpnbins = [Q1n'; Q1p'; Q2n'; Q2p'; Q3n'; Q3p'; Q4n'; Q4p'];
        axs1(i)=subplot(R,C,ax1);
%         plot(Qpnbins(7,:))
%         hold on
%         plot(Qpnbins(8,:))
        bar(c,Qpnbins,'stacked');
        ax1 = ax1 + 1;
    end
    
end
% Bar graph for all day
if strcmp(day_or_frames,'day')
    %Qpnbins = [Q1n Q1p; Q2n Q2p; Q3n Q3p; Q4n Q4p; Q1n+Q2n+Q3n+Q4n Q1p+Q2p+Q3p+Q4p];
    %Qpnbins = [Q1an Q1bn Q1cn Q1dn; Q1ap Q1bp Q1cp Q1dp; Q2an Q2bn Q2cn Q2dn; Q2ap Q2bp Q2cp Q2dp; Q3an Q3bn Q3cn Q3dn; Q3ap Q3bp Q3cp Q3dp; Q4an Q4bn Q4cn Q4dn; Q4ap Q4bp Q4cp Q4dp];
    Qpnbins = [Q1n'; Q1p'; Q2n'; Q2p'; Q3n'; Q3p'; Q4n'; Q4p'];
    figure
%     plot(Qpnbins(7,:))
%     hold on
%     plot(Qpnbins(8,:))
    bar(c,Qpnbins,'stacked');
end
end