////////////////////////////////////////////////////////////////////////////
// 
//   gro 
// 
//   Copyright (c) 2011-2012 Eric Klavins, University of Washington
//   For more information, email klavins@uw.edu
// 
//   This program is free software; you can redistribute it and/or
//   modify it under the terms of the GNU General Public License
//   as published by the Free Software Foundation; either version 2
//   of the License, or (at your option) any later version.
// 
//   This program is distributed in the hope that it will be useful,
//   but WITHOUT ANY WARRANTY; without even the implied warranty of
//   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//   GNU General Public License for more details.
// 
//   You should have received a copy of the GNU General Public License
//   along with this program; if not, write to the Free Software
//   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
// 
// 

// Coded by Shelly Jang.

include gro

// Initialize the simulation parameters
k1 := 5;      // signal diffusion rate
k2 := 0.05;   // signal degradation rate
k3 := 0.1;    // signal emission rate
k4 := 0.02;   // cell death rate; Figure 6 shows the parameter variation simulations of k4. 

// Define a signal molecule
ahl := signal( k1, k2 ); 

program rk0() := {

  true  : {
    emit_signal( ahl, k3 );	// signal is emitted at every dt interval
  };

  rate( k4 * get_signal( ahl ) ) : { 
    die();		      	// cell death
  };

};

program main() := {

  ecoli ( [], program rk0() );	

};