Difference between revisions of "Team:Aix-Marseille/Collaborations"

(Progamming code)
(Progamming code)
Line 24: Line 24:
 
{{hidden
 
{{hidden
 
|Our computational model
 
|Our computational model
/*  
+
|/*  
* Program to run a model of bacteria in a fermenter  
+
| * Program to run a model of bacteria in a fermenter  
* with a 2 plasmid contention system.
+
| * with a 2 plasmid contention system.
*/
+
| */
 
+
|
#include <stdio.h>
+
|#include <stdio.h>
#include <stdlib.h>
+
|#include <stdlib.h>
#include <string.h>
+
|#include <string.h>
#include <stdint.h>
+
|#include <stdint.h>
#include <sys/time.h>
+
|#include <sys/time.h>
#include <gsl/gsl_rng.h>
+
|#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
+
|#include <gsl/gsl_randist.h>
#include <gsl/gsl_sf.h>
+
|#include <gsl/gsl_sf.h>
#include <math.h>
+
|#include <math.h>
 
+
|
#define MAX_POPULATION 2E6
+
|#define MAX_POPULATION 2E6
 
+
|
// Model parameters
+
| // Model parameters
//
+
| //
// Fermentation variables
+
| // Fermentation variables
double Sf    = 10.; // Substrate concentration (g/l) in feed solution
+
|double Sf    = 10.; // Substrate concentration (g/l) in feed solution
double S      =  0.19; // Current substrate concentration in fementer (g/l)
+
|double S      =  0.19; // Current substrate concentration in fementer (g/l)
double V      =  0.01;  // Simulation volume with about 10^6 bacteria in ml
+
|double V      =  0.01;  // Simulation volume with about 10^6 bacteria in ml
double Vtot  = 100.; // Volume of fermenter
+
|double Vtot  = 100.; // Volume of fermenter
 
double D      = 100.; // Dilution rate ml/hr
 
double D      = 100.; // Dilution rate ml/hr
 
double alpha  = 3.4e-11;// growth yield g of substrate needed for 10^6 cells
 
double alpha  = 3.4e-11;// growth yield g of substrate needed for 10^6 cells

Revision as of 09:15, 17 October 2016