# HG changeset patch
# User Michael Guravage <michael.guravage@cwi.nl>
# Date 2011-08-26 17:04:13
# Node ID e730ce1f2d31ab502530c4e41706559b31493d5a
# Parent  9a4aa6470987b8147bd875bd9779eea3f6072727

Added second reaction-diffusion protocol: Fick's laws of diffusion.

 --
 user: Michael Guravage <michael.guravage@cwi.nl>
 branch 'default'
 changed src/protocols/MyDiffusionModel/mydiffusionmodel.cpp

diff --git a/src/protocols/MyDiffusionModel/mydiffusionmodel.cpp b/src/protocols/MyDiffusionModel/mydiffusionmodel.cpp
--- a/src/protocols/MyDiffusionModel/mydiffusionmodel.cpp
+++ b/src/protocols/MyDiffusionModel/mydiffusionmodel.cpp
@@ -31,7 +31,7 @@
 
 QString Mydiffusionmodel::ModelID(void) {
   // specify the name of your model here
-  return QString( "Diffusion dependent coloriing" );
+  return QString( "Fick's Laws of diffusion" );
 }
 
 // return the number of chemicals your model uses
@@ -65,6 +65,14 @@ void Mydiffusionmodel::CellHouseKeeping(
 
 void Mydiffusionmodel::CelltoCellTransport(Wall *w, double *dchem_c1, double *dchem_c2) {
   // add biochemical transport rules here
+
+  // Passive fluxes (Fick's law)
+  for (int c=0;c<NChem();c++) {
+    if (w->C1()->BoundaryPolP() || w->C2()->BoundaryPolP()) return;
+    double phi = w->Length() * ( par->D[c] ) * ( w->C2()->Chemical(c) - w->C1()->Chemical(c) );
+    dchem_c1[c] += phi;
+    dchem_c2[c] -= phi;
+  }
 }
 
 void Mydiffusionmodel::WallDynamics(Wall *w, double *dw1, double *dw2) {