%
% run_sim
%

% Global variables
global num_reactions;
global num_species;
global num_external_species;
global species_names;
global external_species_names;
global max_cells;
global max_alive_cell_index;
global add_species_after_cell_division_index;
global add_species_after_cell_division_amount;
global num_species_to_add_after_cell_division;
global reaction_reactants;
global reaction_reactants_index;
global reaction_reactants_num;
global reaction_products;
global reaction_rates;
global num_reactions;
global external_reaction_reactants;
global external_reaction_products;
global external_reaction_rates;
global num_external_reactions;
global KILLER_THRESHOLD;
global DIVIDE_THRESHOLD;
global num_external_input_reactions;
global external_input_reactions;
global reaction_has_external_inputs;
global curr_time;
global cell_population;
global num_alive_cells;

% Various simulation definitions
global REGRESSION;         % Monte Carlo Simulation Times
global END_TIME;
global RECORD_STEP;
global PLASMID_COPIES;
global CUR_SPECIES_VALUES;
global EXTERNAL_CUR_SPECIES_VALUES;
global REC_SPECIES;
global EXTERNAL_REC_SPECIES;
global AVALS;
global EXTERNAL_AVALS;
global CELL_INTERNAL_AVALS;
global CELL_EXTERNAL_INPUT_AVALS;
global REC_time;

stochastic_sim