Index: /issm/trunk/src/c/objects/Materials/Matice.cpp
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matice.cpp	(revision 8946)
+++ /issm/trunk/src/c/objects/Materials/Matice.cpp	(revision 8947)
@@ -574,4 +574,6 @@
 	/*Intermediaries*/
 	Element *element      = NULL;
+	Parameters* parameters= NULL;
+	int         dim;
 
 	/*Recover element*/
@@ -591,4 +593,20 @@
 					for (int i=0;i<3;i++) values[i]=vector[((Tria*)element)->nodes[i]->GetSidList()]; //use sid list, to index into serial oriented vector 
 					this->inputs->AddInput(new TriaVertexInput(name,values));
+					/*Special case for rheology B in 2D: Pourave land for this solution{{{1*/
+					if(name==RheologyBEnum){
+						/*Are we in 2D?:*/
+						if(element->Enum()==TriaEnum){
+							parameters=((Tria*)(element))->parameters;
+						}
+						else{
+							parameters=((Penta*)(element))->parameters;
+						}
+						parameters->FindParam(&dim,DimEnum);
+						if(dim==2){
+							/*Dupliacte rheology input: */
+							this->inputs->AddInput(new TriaVertexInput(RheologyBbarEnum,values));
+						}
+					}
+					/*}}}*/
 					return;
 
@@ -597,4 +615,7 @@
 		default: _error_("type %i (%s) not implemented yet",type,EnumToStringx(type));
 	}
+
+
+
 }
 /*}}}*/