Index: /issm/trunk/src/c/objects/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Penta.cpp	(revision 3400)
+++ /issm/trunk/src/c/objects/Penta.cpp	(revision 3401)
@@ -250,10 +250,11 @@
 
 	/*Update internal data if inputs holds new values: */
-	inputs->Recover("thickness",&this->properties.h[0],1,dofs,6,(void**)nodes);
-	inputs->Recover("surface",&this->properties.s[0],1,dofs,6,(void**)nodes);
-	inputs->Recover("bed",&this->properties.b[0],1,dofs,6,(void**)nodes);
-	inputs->Recover("drag",&this->properties.k[0],1,dofs,6,(void**)nodes);
-	inputs->Recover("melting",&this->properties.melting[0],1,dofs,6,(void**)nodes);
-	inputs->Recover("accumulation_param",&this->properties.accumulation[0],1,dofs,6,(void**)nodes);
+	if (this->properties.h)              inputs->Recover("thickness",&this->properties.h[0],1,dofs,6,(void**)nodes);
+	if (this->properties.s)              inputs->Recover("surface",&this->properties.s[0],1,dofs,6,(void**)nodes);
+	if (this->properties.b)              inputs->Recover("bed",&this->properties.b[0],1,dofs,6,(void**)nodes);
+	if (this->properties.k)              inputs->Recover("drag",&this->properties.k[0],1,dofs,6,(void**)nodes);
+	if (this->properties.melting)        inputs->Recover("melting",&this->properties.melting[0],1,dofs,6,(void**)nodes);
+	if (this->properties.accumulation)   inputs->Recover("accumulation",&this->properties.accumulation[0],1,dofs,6,(void**)nodes);
+	if (this->properties.geothermalflux) inputs->Recover("geothermalflux",&this->properties.geothermalflux[0],1,dofs,6,(void**)nodes);
 
 	//Update material if necessary