%
% create_model_simple
%

% Variables
global reaction_reactants;
global reaction_products;
global reaction_rates;
global num_reactions;
global max_reactions;
global species_names;
global num_species_to_add_after_cell_division;

num_species_to_add_after_cell_division = 0;

max_reactions = 100;
num_reactions = 0;

reaction_reactants = {''};
reaction_products  = {''};
reaction_rates     = {''};
species_names      = {};

%reactions      = zeros(max_reactions,5,5);  % for now, max reactants = 5 and max products = 5
%reaction_rates = zeros(max_reactions,1);

% create the reactions for each cell
add_reaction({'P_lux'},        {'P_lux' 'LuxI'}, 0.3);
add_reaction({'LuxI'},         {}, 1.0);
add_reaction({'Alive'},        {'Alive' 'Volume'}, 10.0);
add_reaction({'create_P_lux'}, {'P_lux'}, 1.0);

% special instructions for dealing with DNA replication after cell division
add_after_cell_division('create_P_lux', 10);