Index: /issm/trunk/src/c/objects/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Penta.cpp	(revision 877)
+++ /issm/trunk/src/c/objects/Penta.cpp	(revision 878)
@@ -1099,15 +1099,19 @@
 	inputs->Recover("melting",&melting[0],1,dofs,6,(void**)nodes);
 	inputs->Recover("accumulation_param",&accumulation[0],1,dofs,6,(void**)nodes);
-	
+
 	//Update material if necessary
 	if(inputs->Recover("temperature_average",&temperature_list[0],1,dofs,6,(void**)nodes)){
-		temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
-		B_average=Paterson(temperature_average);
-		matice->SetB(B_average);
+		if(matice){
+			temperature_average=(temperature_list[0]+temperature_list[1]+temperature_list[2]+temperature_list[3]+temperature_list[4]+temperature_list[5])/6.0;
+			B_average=Paterson(temperature_average);
+			matice->SetB(B_average);
+		}
 	}
 	
 	if(inputs->Recover("B",&B_list[0],1,dofs,6,(void**)nodes)){
-		B_average=(B_list[0]+B_list[1]+B_list[2]+B_list[3]+B_list[4]+B_list[5])/6.0;
-		matice->SetB(B_average);
+		if(matice){
+			B_average=(B_list[0]+B_list[1]+B_list[2]+B_list[3]+B_list[4]+B_list[5])/6.0;
+			matice->SetB(B_average);
+		}
 	}