clear, %close all
tic
clc
%for k = 1:24
antenna_choices = {'11','21','22'};
channel_choices = {'AWGN','AWGN+Baseline','AWGN+Rayleigh','AWGN+RayleighBlock'};
channel_choices = {'AWGN+Rayleigh'};
global seed
set_no = 0;
for antenna_choice = 1:length(antenna_choices)
    for channel_choice = 1:length(channel_choices)
        clearvars -except channel_choice channel_choices antenna_choice antenna_choices seed set_no
        set_no = set_no + 1;
        disp(['Starting Set No: ' , num2str(set_no)])

        seed = 1;
        mastermfile = mfilename;
        extra_comment = [mfilename ' extra comment here'];
        % this variable is for logging an extra comment when the results are
        % written on disk: Write the aim of simulation in that variable

        %% 1.1 Common simulation parameters for each burst
        handles_common.fft = '1024'; % Choices: {2048, 1024}
        handles_common.EsNo = [0:1:10]; % EbNo vector
        %handles_common.dectype = 'LLR';
        handles_common.trial_no = [20]; % 5 5 10 10 10 20 20 20 50 50 50 100]; %[ones(1,10)*50 ones(1,10)*300]; [50 50 50 100 100 500]; %[1     1     4     7    12    20]; %randint(size(handles_common.EsNo,1), size(handles_common.EsNo,2), [2 6]); % no. of trials at each EbNo
        %handles_common.trial_no = 1; %randint(size(handles_common.EsNo,1), size(handles_common.EsNo,2), [2 6]); % no. of trials at each EbNo
        %handles_common.trial_no = ones(size(handles_common.EsNo)) *

        %% 1.2 Common channel and antenna configurations
        handles_common.antenna = antenna_choices{antenna_choice}; % 'SISO', 'MIMO - 2x1', 'MIMO- 2x2','11','21','22' are all valid
        demapchoice = 1;
        handles_common.channel = channel_choices{channel_choice}; %'AWGN+Rayleigh'; % 'AWGN + Baseline', 'AWGN + No Baseline',  'Baseline + No AWGN', 'AWGN + Rayleigh','No Channel'
        %handles_common.channel = 'nochannel'; % 'AWGN + Baseline', 'AWGN + No Baseline',  'Baseline + No AWGN', 'AWGN + Rayleigh','No Channel'
        handles_common.v = 3;
        handles_common.scenario = 3; % 1,2,3,4
        handles_common.timedomain = 0;
        handles_common.scatterplot = 0;
        handles_common.sameconf = 0;

        %% 1.3 Plot parameters
        plot_param.log = 1; % 0 or 1
        plot_param.hold = 1; % 0 or 1

        %% 2.1 Process variables for each burst

        burst_codes_dl = {'C17','D16'}; %,'D16'}; %,num2str(jjj)],['D',num2str(jjj)], ['E',num2str(jjj)]};
        burst_codes_ul = {'C1'}; % {'A1'};

        %% Burst frequency and OFDM symbol allocations

        DL_burst(1,:) = [0, 22, 0, 10];
        DL_burst(2,:) = [0, 22, 10, 10];
        DL_burst(3,:) = [0, 22, 20, 10];
        DL_burst(4,:) = [0, 22, 15, 0];
        DL_burst(5,:) = [0, 22, 20, 0];
        DL_burst(6,:) = [0, 22, 25, 0];


        % Original burst configuration -- copy and paste from block if needed
        %     DL_burst(1,:) = [0, 22, 0, 5];
        %     DL_burst(2,:) = [0, 22, 5, 5];
        %     DL_burst(3,:) = [0, 22, 10, 5];
        %     DL_burst(4,:) = [0, 22, 15, 5];
        %     DL_burst(5,:) = [0, 22, 20, 5];
        %     DL_burst(6,:) = [0, 22, 25, 5];

        UL_burst(1,:) = [0, 9, 0, 5];
        UL_burst(2,:) = [0, 9, 5, 5];
        UL_burst(3,:) = [0, 9, 10, 5];
        UL_burst(4,:) = [0, 9, 15, 5];
        UL_burst(5,:) = [0, 9, 20, 5];
        UL_burst(6,:) = [0, 9, 25, 5];
        scenarios_common
        % after the variables are defined, perform the common tasks to complete
        % the simulation of a scenario
    end
end

toc