From 2007.igem.org

Cell auto

!include "gbviewOutput.mgs" ;;

N := 30;;

gbf grille = <nord, est ; 30 * nord, 30 * est> and record Bact = { DAP:float, DAPe:float } and record BactG = Bact + { bactG } and record BactS = Bact + { bactS} ;;

fun iota_v(v,n) = map((\x.v), iota(n,seq:())) ;; fun new_grid(v) = iota_v(iota_v(v,N), N) following |nord>, |est> ;;

fun majBactG[DAPexport = 0.2, DAPimport = 0.2, DAPconso= 1, SupDiff = 2.4, InfDiff = 2.02 ](x) =

	if ( ( x.DAP > InfDiff ) & ( x.DAP < SupDiff ) & (random(1) < 0.8) ) then

{DAP = x.DAP , DAPe = x.DAPe, bactS} else let dap = if x.DAP - DAPconso > 0 then x.DAP - DAPconso else 0 fi and dape= x.DAPe in if dap > dape then

                   x + {DAP = dap - DAPexport * dap, DAPe = dape + DAPexport * dap }

else x + {DAP = dap + DAPimport * dape, DAPe = dape - DAPimport * dape }

                 fi
       fi 

fun majBactS[DAPexport = 0.2, DAPimport = 0.2, DAPprod=10](x) = let dap = x.DAP + DAPprod and dape= x.DAPe in if dap > dape then

                   x + {DAP = dap - DAPexport * dap, DAPe = dape + DAPexport * dap }

else x + {DAP = dap + DAPimport * dape, DAPe = dape - DAPimport * dape }

                 fi

trans evol[Delta_t=0.1,DAPeDiff=1, DAPeDegrad = 0.2] = {

 (* production of DAP *)
 x:BactS => ( 

let d = neighborsfold(

     (\y.\acc.( DAPeDiff*Delta_t*(y.DAPe-x.DAPe) + acc)),
     x.DAPe,
     x)
   in majBactS( x + {DAPe = if ( d - DAPeDegrad * x.DAPe) > 0 then d - DAPeDegrad * x.DAPe else 0 fi } )
     );

 (* Diffirentiation and Consumption *)
x:BactG => ( 

let d = neighborsfold(

     (\y.\acc.( DAPeDiff*Delta_t*(y.DAPe-x.DAPe) + acc)),
     x.DAPe,
     x)
   in majBactG(x + {DAPe = if ( d - DAPeDegrad * x.DAPe) > 0 then  d - DAPeDegrad * x.DAPe else 0 fi } )
   

 )

} ;;

g := new_grid({DAP=6.0, DAPe=3.0, bactG}) ; 0 ;; g := set_gbfpos(g, (random(30)*|nord> + random(30)*|est>), {DAP = 6.0, DAPe=3.0, bactS}) ;; g := set_gbfpos(g, (random(30)*|nord> + random(30)*|est>), {DAP = 6.0, DAPe=3.0, bactS}) ;; g := set_gbfpos(g, (random(30)*|nord> + random(30)*|est>), {DAP = 6.0, DAPe=3.0, bactS}) ;; g := set_gbfpos(g, (random(30)*|nord> + random(30)*|est>), {DAP = 6.0, DAPe=3.0, bactS}) ;; evol[iter=1000,interlude=GBVexport((\v.("DAPe = " + (12 * v.DAPe ) + if BactS(v) then ", bactS= 255" else ", bactG = 255" fi )))](g) ;;