Index: /issm/trunk/src/c/DataSet/DataSet.cpp
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3965)
+++ /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3966)
@@ -549,55 +549,4 @@
 
 /*Objects methods*/
-/*FUNCTION DataSet::ComputeBasalStress{{{1*/
-void DataSet::ComputeBasalStress(Vec sigma_b,int analysis_type,int sub_analysis_type){
-
-	vector<Object*>::iterator object;
-	Element* element=NULL;
-
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		if(EnumIsElement((*object)->Enum())){
-
-			element=(Element*)(*object);
-			element->ComputeBasalStress(sigma_b,analysis_type,sub_analysis_type);
-		}
-	}
-
-}
-/*}}}*/
-/*FUNCTION DataSet::ComputePressure{{{1*/
-void DataSet::ComputePressure(Vec p_g,int analysis_type,int sub_analysis_type){
-
-	vector<Object*>::iterator object;
-	Element* element=NULL;
-
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		if(EnumIsElement((*object)->Enum())){
-
-			element=(Element*)(*object);
-			element->ComputePressure(p_g,analysis_type,sub_analysis_type);
-		}
-	}
-
-}
-/*}}}*/
-/*FUNCTION DataSet::ComputeStrainRate{{{1*/
-void DataSet::ComputeStrainRate(Vec eps,int analysis_type,int sub_analysis_type){
-
-	vector<Object*>::iterator object;
-	Element* element=NULL;
-
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		if(EnumIsElement((*object)->Enum())){
-
-			element=(Element*)(*object);
-			element->ComputeStrainRate(eps,analysis_type,sub_analysis_type);
-		}
-	}
-
-}
-/*}}}*/
 /*FUNCTION DataSet::Configure{{{1*/
 void DataSet::Configure(DataSet* elements,DataSet* loads, DataSet* nodes, DataSet* vertices, DataSet* materials,Parameters* parameters){
Index: /issm/trunk/src/c/DataSet/DataSet.h
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.h	(revision 3965)
+++ /issm/trunk/src/c/DataSet/DataSet.h	(revision 3966)
@@ -91,7 +91,4 @@
 		void  InputExtrude(int enum_type);
 		int   DeleteObject(Object* object);
-		void  ComputeBasalStress(Vec sigma_b,int analysis_type,int sub_analysis_type);
-		void  ComputePressure(Vec p_g,int analysis_type,int sub_analysis_type);
-		void  ComputeStrainRate(Vec eps,int analysis_type,int sub_analysis_type);
 		void  FieldExtrude(Vec field,double* field_serial,char* field_name, int collapse);
 		void  InputToResult(Mat* psolution,int enum_type);
Index: /issm/trunk/src/c/modules/ComputeBasalStressx/ComputeBasalStressx.cpp
===================================================================
--- /issm/trunk/src/c/modules/ComputeBasalStressx/ComputeBasalStressx.cpp	(revision 3965)
+++ /issm/trunk/src/c/modules/ComputeBasalStressx/ComputeBasalStressx.cpp	(revision 3966)
@@ -13,8 +13,9 @@
 			int analysis_type,int sub_analysis_type){
 
+	/*Intermediary*/
 	int i;
-
-	int  found=0;
-	double numberofelements;
+	int      found=0;
+	double   numberofelements;
+	Element* element=NULL;
 
 	/*output: */
@@ -29,10 +30,13 @@
 
 	/*Get elements configured: */
-	elements->Configure(elements,loads, nodes,vertices, materials,parameters);
-	nodes->Configure(elements,loads, nodes,vertices, materials,parameters);
+	elements->  Configure(elements,loads, nodes,vertices, materials,parameters);
+	nodes->     Configure(elements,loads, nodes,vertices, materials,parameters);
 	parameters->Configure(elements,loads, nodes,vertices, materials,parameters);
 
-	/*Call on dataset driver: */
-	elements->ComputeBasalStress(sigma,analysis_type,sub_analysis_type);
+	/*Compute basal stress for each element: */
+	for (i=0;i<elements->Size();i++){
+		element=(Element*)elements->GetObjectByOffset(i);
+		element->ComputeBasalStress(sigma,analysis_type,sub_analysis_type);
+	}
 
 	/*Assemble vector: */
Index: /issm/trunk/src/c/modules/ComputePressurex/ComputePressurex.cpp
===================================================================
--- /issm/trunk/src/c/modules/ComputePressurex/ComputePressurex.cpp	(revision 3965)
+++ /issm/trunk/src/c/modules/ComputePressurex/ComputePressurex.cpp	(revision 3966)
@@ -13,5 +13,8 @@
 			int analysis_type,int sub_analysis_type){
 
+	/*intermediary*/
 	int numberofnodes;
+	int i;
+	Element* element=NULL;
 
 	/*output: */
@@ -30,6 +33,9 @@
 	parameters->Configure(elements,loads,nodes,vertices,materials,parameters);
 
-	/*Call on dataset driver: */
-	elements->ComputePressure(p_g,analysis_type,sub_analysis_type);
+	/*Compute pressure for each element: */
+	for (i=0;i<elements->Size();i++){
+		element=(Element*)elements->GetObjectByOffset(i);
+		element->ComputePressure(p_g,analysis_type,sub_analysis_type);
+	}
 
 	/*Assemble vector: */
Index: /issm/trunk/src/c/modules/ComputeStrainRatex/ComputeStrainRatex.cpp
===================================================================
--- /issm/trunk/src/c/modules/ComputeStrainRatex/ComputeStrainRatex.cpp	(revision 3965)
+++ /issm/trunk/src/c/modules/ComputeStrainRatex/ComputeStrainRatex.cpp	(revision 3966)
@@ -13,8 +13,9 @@
 			int analysis_type,int sub_analysis_type){
 
+	/*Intermediary*/
 	int i;
-
-	int  found=0;
-	double numberofelements;
+	int      found=0;
+	double   numberofelements;
+	Element* element=NULL;
 
 	/*output: */
@@ -28,10 +29,13 @@
 
 	/*Get elements configured: */
-	elements->Configure(elements,loads, nodes,vertices, materials,parameters);
-	nodes->Configure(elements,loads, nodes,vertices, materials,parameters);
+	elements->  Configure(elements,loads, nodes,vertices, materials,parameters);
+	nodes->     Configure(elements,loads, nodes,vertices, materials,parameters);
 	parameters->Configure(elements,loads, nodes,vertices, materials,parameters);
 
-	/*Call on dataset driver: */
-	elements->ComputeStrainRate(eps,analysis_type,sub_analysis_type);
+	/*Compute basal stress for each element: */
+	for (i=0;i<elements->Size();i++){
+		element=(Element*)elements->GetObjectByOffset(i);
+		element->ComputeStrainRate(eps,analysis_type,sub_analysis_type);
+	}
 
 	/*Assemble vector: */