diff --git a/src/build_models/meinhardtplugin.cpp b/src/build_models/meinhardtplugin.cpp
--- a/src/build_models/meinhardtplugin.cpp
+++ b/src/build_models/meinhardtplugin.cpp
@@ -41,90 +41,86 @@ void MeinhardtPlugin::OnDivide(ParentInf
 
 void MeinhardtPlugin::SetCellColor(CellBase *c, QColor *color) { 
 
-	if (fpclassify(c->Chemical(0))==FP_NAN) {
-		// somehow the function isnan doesn't work properly on my system... SuSE Linux
-		// 10.0 64-bits (isnan seems not be implemented using fpclassify).
-		MyWarning::warning("Whoops! Numerical instability!!");
-		color->setNamedColor("red");
-	} else {
-		double range_min = 0.;//, range_max = 1.;
-		if (c->Chemical(0)<range_min) {
-			MyWarning::warning("Whoops! Numerical instability!!");
-			color->setNamedColor("blue");
-		} else {
-			color->setRgb(c->Chemical(1)/(1+c->Chemical(1)) * 255.,(c->Chemical(0)/(1+c->Chemical(0)) * 255.),(c->Chemical(3)/(1+c->Chemical(3)) *255.) );
-		}
-		
-	}
+  if (fpclassify(c->Chemical(0))==FP_NAN) {
+    // somehow the function isnan doesn't work properly on my system... SuSE Linux
+    // 10.0 64-bits (isnan seems not be implemented using fpclassify).
+    MyWarning::warning("Whoops! Numerical instability!!");
+    color->setNamedColor("red");
+  } else {
+    double range_min = 0.;//, range_max = 1.;
+    if (c->Chemical(0)<range_min) {
+      MyWarning::warning("Whoops! Numerical instability!!");
+      color->setNamedColor("blue");
+    } else {
+      color->setRgb(c->Chemical(1)/(1+c->Chemical(1)) * 255.,(c->Chemical(0)/(1+c->Chemical(0)) * 255.),(c->Chemical(3)/(1+c->Chemical(3)) *255.) );
+    }
+  }
 }
 
 
 
 void MeinhardtPlugin::CellHouseKeeping(CellBase *c) {
 	
-	if (c->Area() > par->rel_cell_div_threshold * c->BaseArea() ) {
-		c->Divide();
-	}
+  if (c->Area() > par->rel_cell_div_threshold * c->BaseArea() ) {
+    c->Divide();
+  }
 	
-    // cell expansion is inhibited by substrate (chem 3)
-	if (!par->constituous_expansion_limit || c->NCells()<par->constituous_expansion_limit) {
-		c->EnlargeTargetArea(par->cell_expansion_rate );
-	} else {
-		if (c->Chemical(0)<0.5) {
-			double tmp;
-			c->EnlargeTargetArea((tmp=(1.-par->vessel_inh_level*c->Chemical(3))*par->cell_expansion_rate /* + c->Chemical(4)*/)<0?0:tmp); 
-		} else {
-			c->EnlargeTargetArea(par->vessel_expansion_rate);
-		}
-	} 
-    
+  // cell expansion is inhibited by substrate (chem 3)
+  if (!par->constituous_expansion_limit || c->NCells()<par->constituous_expansion_limit) {
+    c->EnlargeTargetArea(par->cell_expansion_rate );
+  } else {
+    if (c->Chemical(0)<0.5) {
+      double tmp;
+      c->EnlargeTargetArea((tmp=(1.-par->vessel_inh_level*c->Chemical(3))*par->cell_expansion_rate /* + c->Chemical(4)*/)<0?0:tmp); 
+    } else {
+      c->EnlargeTargetArea(par->vessel_expansion_rate);
+    }
+  } 
 }
 
 void MeinhardtPlugin::CelltoCellTransport(Wall *w, double *dchem_c1, double *dchem_c2) {
 	
-	// No flux boundaries for all chemicals, except activator: boundary is sink
+  // No flux boundaries for all chemicals, except activator: boundary is sink
   if (w->C1()->BoundaryPolP() || w->C2()->BoundaryPolP()) {
 		
-		if (w->C1()->BoundaryPolP()) {
-			dchem_c2[1] -=  w->Length() * ( par->D[1] ) * ( w->C2()->Chemical(1) );
-			} else {
-			dchem_c1[1] -=  w->Length() * ( par->D[1] ) * ( w->C1()->Chemical(1) );
-			} 
-		return;
-		}
+    if (w->C1()->BoundaryPolP()) {
+      dchem_c2[1] -=  w->Length() * ( par->D[1] ) * ( w->C2()->Chemical(1) );
+    } else {
+      dchem_c1[1] -=  w->Length() * ( par->D[1] ) * ( w->C1()->Chemical(1) );
+    } 
+    return;
+  }
 	
 
-    // Passive fluxes (Fick's law)
-    for (int c=0;c<NChem();c++) {
-		double phi = w->Length() * ( par->D[c] ) * ( w->C2()->Chemical(c) - w->C1()->Chemical(c) );
-		dchem_c1[c] += phi; 
-		dchem_c2[c] -= phi;
-	}
-   
-
+  // Passive fluxes (Fick's law)
+  for (int c=0;c<NChem();c++) {
+    double phi = w->Length() * ( par->D[c] ) * ( w->C2()->Chemical(c) - w->C1()->Chemical(c) );
+    dchem_c1[c] += phi; 
+    dchem_c2[c] -= phi;
+  }
 }
 
 void MeinhardtPlugin::WallDynamics(Wall *w, double *dw1, double *dw2) {
   w = NULL;
-	for (int c = 0;c<NChem();c++) {
-		dw1[c] = 0.; dw2[c] = 0.;
-    }
+  for (int c = 0;c<NChem();c++) {
+    dw1[c] = 0.; dw2[c] = 0.;
+  }
 }
 
 void MeinhardtPlugin::CellDynamics(CellBase *c, double *dchem) {
 
-	double Y = c->Chemical(0);
-	double A = c->Chemical(1);
-	double H = c->Chemical(2);
-	double S = c->Chemical(3);
-
+  double Y = c->Chemical(0);
+  double A = c->Chemical(1);
+  double H = c->Chemical(2);
+  double S = c->Chemical(3);
 	
-	
-    dchem[0] = ( par->d * A - par->e * Y + Y*Y/(1 + par->f * Y*Y ) );
-	dchem[1] = ( par->c * A*A*S/H - par->mu * A + par->rho0*Y );
-	dchem[2] = ( par->c * A*A*S - par->nu*H + par->rho1*Y );
-	dchem[3] = ( par->c0 - par->gamma*S - par->eps * Y * S );
+  dchem[0] = ( par->d * A - par->e * Y + Y*Y/(1 + par->f * Y*Y ) );
+  dchem[1] = ( par->c * A*A*S/H - par->mu * A + par->rho0*Y );
+  dchem[2] = ( par->c * A*A*S - par->nu*H + par->rho1*Y );
+  dchem[3] = ( par->c0 - par->gamma*S - par->eps * Y * S );
 }
 
 
 Q_EXPORT_PLUGIN2(meinhardtplugin, MeinhardtPlugin)
+
+/* finis */