Index: /issm/trunk/src/c/Container/Inputs.cpp
===================================================================
--- /issm/trunk/src/c/Container/Inputs.cpp	(revision 4236)
+++ /issm/trunk/src/c/Container/Inputs.cpp	(revision 4237)
@@ -401,5 +401,5 @@
 	/*Compute the 3d Strain Rate (6 components):
 	 *
-	 * epsilon=[exx eyy ezz exy exz euz]
+	 * epsilon=[exx eyy ezz exy exz eyz]
 	 */
 
@@ -460,5 +460,5 @@
 	/*Compute the 3d Blatter/PattynStrain Rate (5 components):
 	 *
-	 * epsilon=[exx eyy exy exz euz]
+	 * epsilon=[exx eyy exy exz eyz]
 	 *
 	 * with exz=1/2 du/dz
Index: /issm/trunk/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp
===================================================================
--- /issm/trunk/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp	(revision 4236)
+++ /issm/trunk/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp	(revision 4237)
@@ -38,4 +38,6 @@
 	IoModelFetchData(&iomodel->upperelements,NULL,NULL,iomodel_handle,"upperelements");
 	IoModelFetchData(&iomodel->lowerelements,NULL,NULL,iomodel_handle,"lowerelements");
+	IoModelFetchData(&iomodel->rheology_B,NULL,NULL,iomodel_handle,"rheology_B");
+	IoModelFetchData(&iomodel->rheology_n,NULL,NULL,iomodel_handle,"rheology_n");
 	
 	/*Create elements and materials: */
@@ -57,4 +59,6 @@
 	xfree((void**)&iomodel->upperelements);
 	xfree((void**)&iomodel->lowerelements);
+	xfree((void**)&iomodel->rheology_B);
+	xfree((void**)&iomodel->rheology_n);
 
 	/*Add new constrant material property tgo materials, at the end: */
Index: /issm/trunk/src/c/objects/Materials/Matice.cpp
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matice.cpp	(revision 4236)
+++ /issm/trunk/src/c/objects/Materials/Matice.cpp	(revision 4237)
@@ -231,4 +231,9 @@
 	}
 
+	/*Checks in debugging mode*/
+	ISSMASSERT(B>0);
+	ISSMASSERT(n>0);
+	ISSMASSERT(viscosity>0);
+
 	/*Return: */
 	*pviscosity=viscosity;
@@ -291,4 +296,9 @@
 		}
 	}
+
+	/*Checks in debugging mode*/
+	ISSMASSERT(B>0);
+	ISSMASSERT(n>0);
+	ISSMASSERT(viscosity3d>0);
 
 	/*Assign output pointers:*/
@@ -356,4 +366,9 @@
 	}
 
+	/*Checks in debugging mode*/
+	ISSMASSERT(B>0);
+	ISSMASSERT(n>0);
+	ISSMASSERT(viscosity3d>0);
+
 	/*Assign output pointers:*/
 	*pviscosity3d=viscosity3d;
@@ -405,4 +420,9 @@
 		viscosity_complement=4.5*pow((double)10,(double)17);
 	}
+
+	/*Checks in debugging mode*/
+	ISSMASSERT(B>0);
+	ISSMASSERT(n>0);
+	ISSMASSERT(viscosity_complement>0);
 		
 	/*Return: */