Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4432)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4433)
@@ -5168,5 +5168,5 @@
 }
 /*}}}*/
-/*FUNCTION Tria::GetParameterValue( {{{1*/
+/*FUNCTION Tria::GetParameterValue(double* pvalue,Node* node1,Node* node2,double gauss_seg,int enumtype) {{{1*/
 void Tria::GetParameterValue(double* pvalue,Node* node1,Node* node2,double gauss_seg,int enumtype){
 
@@ -5182,4 +5182,5 @@
 
 	/*go through 3 nodes (all nodes for tria) and identify 1st and 2nd nodes: */
+	ISSMASSERT(nodes);
 	for(i=0;i<numnodes;i++){
 		if (node1==nodes[i]) grid1=i;
Index: /issm/trunk/src/c/objects/Loads/Numericalflux.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Numericalflux.cpp	(revision 4432)
+++ /issm/trunk/src/c/objects/Loads/Numericalflux.cpp	(revision 4433)
@@ -23,5 +23,5 @@
 	this->inputs=NULL;
 	this->parameters=NULL;
-	this->helements=NULL;
+	this->helement=NULL;
 	this->hnodes=NULL;
 }
@@ -111,6 +111,6 @@
 
 	/*Hooks: */
-	this->hnodes=new Hook(numericalflux_node_ids,num_nodes);
-	this->helements=new Hook(numericalflux_elem_ids,num_elems);
+	this->hnodes  =new Hook(numericalflux_node_ids,num_nodes);
+	this->helement=new Hook(numericalflux_elem_ids,1); // take only the first element for now
 
 	//intialize  and add as many inputs per element as requested: 
@@ -126,5 +126,5 @@
 	delete inputs;
 	this->parameters=NULL;
-	delete helements;
+	delete helement;
 	delete hnodes;
 }
@@ -144,5 +144,5 @@
 	printf("   analysis_type: %s\n",EnumAsString(analysis_type));
 	hnodes->DeepEcho();
-	helements->DeepEcho();
+	helement->DeepEcho();
 	printf("   parameters\n");
 	parameters->DeepEcho();
@@ -186,5 +186,5 @@
 	/*Marshall hooks: */
 	hnodes->Marshall(&marshalled_dataset);
-	helements->Marshall(&marshalled_dataset);
+	helement->Marshall(&marshalled_dataset);
 
 	/*Marshall inputs: */
@@ -208,5 +208,5 @@
 		+sizeof(analysis_type)
 		+hnodes->MarshallSize()
-		+helements->MarshallSize()
+		+helement->MarshallSize()
 		+inputs->MarshallSize()
 		+sizeof(int); //sizeof(int) for enum type
@@ -229,5 +229,5 @@
 	/*demarshall hooks: */
 	hnodes=new Hook(); hnodes->Demarshall(&marshalled_dataset);
-	helements=new Hook(); helements->Demarshall(&marshalled_dataset);
+	helement=new Hook(); helement->Demarshall(&marshalled_dataset);
 	
 	/*demarshall inputs: */
@@ -270,5 +270,5 @@
 	/*now deal with hooks and objects: */
 	numericalflux->hnodes=(Hook*)this->hnodes->copy();
-	numericalflux->helements=(Hook*)this->helements->copy();
+	numericalflux->helement=(Hook*)this->helement->copy();
 
 	return numericalflux;
@@ -284,5 +284,5 @@
 	 * datasets, using internal ids and offsets hidden in hooks: */
 	hnodes->configure(nodesin);
-	helements->configure(elementsin);
+	helement->configure(elementsin);
 
 	/*point parameters to real dataset: */
@@ -845,5 +845,5 @@
 
 	/*recover objects from hooks: */
-	tria=(Tria*)helements[0].deliverp();
+	tria=(Tria*)helement->delivers();
 	nodes=(Node**)hnodes->deliverp();
 
Index: /issm/trunk/src/c/objects/Loads/Numericalflux.h
===================================================================
--- /issm/trunk/src/c/objects/Loads/Numericalflux.h	(revision 4432)
+++ /issm/trunk/src/c/objects/Loads/Numericalflux.h	(revision 4433)
@@ -22,5 +22,5 @@
 		int analysis_type;
 
-		Hook*       helements;
+		Hook*       helement;
 		Hook*       hnodes;
 		Parameters *parameters;