Index: /issm/trunk/src/c/objects/Loads/Riftfront.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Riftfront.cpp	(revision 4695)
+++ /issm/trunk/src/c/objects/Loads/Riftfront.cpp	(revision 4696)
@@ -35,4 +35,5 @@
 	/*data: */
 	int    riftfront_node_ids[2];
+	int    riftfront_elem_ids[2];
 	int    riftfront_matpar_id;
 	int    riftfront_type;
@@ -60,7 +61,10 @@
 	riftfront_node_ids[0]=iomodel->nodecounter+grid1;
 	riftfront_node_ids[1]=iomodel->nodecounter+grid2;
+	riftfront_elem_ids[0]=el1;
+	riftfront_elem_ids[1]=el2;
 	riftfront_matpar_id=iomodel->numberofelements+1; //matlab indexing
 
 	this->hnodes=new Hook(riftfront_node_ids,2);
+	this->helements=new Hook(riftfront_elem_ids,2);
 	this->hmatpar=new Hook(&riftfront_matpar_id,1);
 
@@ -102,4 +106,5 @@
 	this->parameters=NULL;
 	delete hnodes;
+	delete helements;
 	delete hmatpar;
 	delete helements;
@@ -120,4 +125,5 @@
 	printf("   analysis_type: %s\n",EnumAsString(analysis_type));
 	hnodes->DeepEcho();
+	helements->DeepEcho();
 	hmatpar->DeepEcho();
 	printf("   parameters\n");
@@ -168,4 +174,5 @@
 	/*Marshall hooks: */
 	hnodes->Marshall(&marshalled_dataset);
+	helements->Marshall(&marshalled_dataset);
 	hmatpar->Marshall(&marshalled_dataset);
 
@@ -199,4 +206,5 @@
 		+sizeof(penalty_lock)
 		+hnodes->MarshallSize()
+		+helements->MarshallSize()
 		+hmatpar->MarshallSize()
 		+inputs->MarshallSize()
@@ -231,4 +239,5 @@
 	/*demarshall hooks: */
 	hnodes=new Hook(); hnodes->Demarshall(&marshalled_dataset);
+	helements=new Hook(); helements->Demarshall(&marshalled_dataset);
 	hmatpar=new Hook(); hmatpar->Demarshall(&marshalled_dataset);
 	
@@ -271,4 +280,5 @@
 	/*now deal with hooks and objects: */
 	riftfront->hnodes=(Hook*)this->hnodes->copy();
+	riftfront->helements=(Hook*)this->helements->copy();
 	riftfront->hmatpar=(Hook*)this->hmatpar->copy();
 
@@ -285,4 +295,5 @@
 	 * datasets, using internal ids and offsets hidden in hooks: */
 	hnodes->configure(nodesin);
+	helements->configure(elementsin);
 	hmatpar->configure(materialsin);
 
@@ -298,4 +309,5 @@
 	 * datasets, using internal ids and offsets hidden in hooks: */
 	hnodes->configure(nodesin);
+	helements->configure(elementsin);
 	hmatpar->configure(materialsin);