Index: /issm/trunk/src/c/objects/IoModel.cpp
===================================================================
--- /issm/trunk/src/c/objects/IoModel.cpp	(revision 3812)
+++ /issm/trunk/src/c/objects/IoModel.cpp	(revision 3813)
@@ -112,4 +112,5 @@
 	xfree((void**)&this->my_bordervertices);
 	xfree((void**)&this->penaltypartitioning);
+	xfree((void**)&this->singlenodetoelementconnectivity);
 	#endif
 }
@@ -363,5 +364,4 @@
 	this->numpenalties=0;
 	this->penalties=NULL;
-	this->penaltypartitioning=NULL;
 
 	/*!basal: */
@@ -383,4 +383,5 @@
 	this->my_bordervertices=NULL;
 	this->penaltypartitioning=NULL;
+	this->singlenodetoelementconnectivity=NULL;
 }
 /*}}}*/
Index: /issm/trunk/src/c/objects/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 3812)
+++ /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 3813)
@@ -85,4 +85,6 @@
 	/*hooks: */
 	pengrid_node_id=index+1;
+	ISSMASSERT(iomodel->singlenodetoelementconnectivity);
+	ISSMASSERT(index>=0 && index<iomodel->numberofvertices);
 	pengrid_element_id=iomodel->singlenodetoelementconnectivity[index];
 	pengrid_matpar_id=iomodel->numberofelements+1; //refers to the constant material parameters object