Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 4450)
+++ /issm/trunk/src/c/Makefile.am	(revision 4451)
@@ -1077,5 +1077,4 @@
 					./modules/Qmux/Qmux.h\
 					./modules/Qmux/Qmux.cpp\
-					./modules/Qmux/DakotaResponses.cpp\
 					./modules/Qmux/DakotaMPI_Bcast.cpp\
 					./modules/Qmux/DakotaFree.cpp\
Index: /issm/trunk/src/c/modules/DakotaResponsesx/DakotaResponsesx.cpp
===================================================================
--- /issm/trunk/src/c/modules/DakotaResponsesx/DakotaResponsesx.cpp	(revision 4450)
+++ /issm/trunk/src/c/modules/DakotaResponsesx/DakotaResponsesx.cpp	(revision 4451)
@@ -35,38 +35,38 @@
 
 		if(strcmp(response_descriptor,"min_vel")==0){
-			MinVelx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MinVelx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_vel")==0){
-			MaxVelx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxVelx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"min_vx")==0){
-			MinVxx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MinVxx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_vx")==0){
-			MaxVxx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxVxx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_abs_vx")==0){
-			MaxAbsVxx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxAbsVxx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"min_vy")==0){
-			MinVyx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MinVyx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_vy")==0){
-			MaxVyx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxVyx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_abs_vy")==0){
-			MaxAbsVyx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxAbsVyx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"min_vz")==0){
-			MinVzx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MinVzx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_vz")==0){
-			MaxVzx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxVzx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"max_abs_vz")==0){
-			MaxAbsVzx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			MaxAbsVzx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"misfit")==0){
-			Misfitx( &femmodel_response, femmodel->elements,femmodel->nodes, femmodel->vertices, femmodel->loads, femmodel->materials, femmodel->parameters);
+			Misfitx( &femmodel_response, elements,nodes, vertices, loads, materials, parameters);
 		}
 		else if(strcmp(response_descriptor,"mass_flux")==0){
@@ -76,8 +76,8 @@
 			
 			/*retrieve qmu_mass_flux_segments: */
-			femmodel->parameters->FindParam(&segments,&num_segments,&dummy,QmuMassFluxSegmentsEnum);
+			parameters->FindParam(&segments,&num_segments,&dummy,QmuMassFluxSegmentsEnum);
 
 			/*call mass flux module: */
-			MassFluxx(&femmodel_response,femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,segments,num_segments);
+			MassFluxx(&femmodel_response,elements,nodes,vertices,loads,materials,parameters,segments,num_segments);
 
 			/*Free ressources:*/
