Index: /issm/trunk/src/c/CostFunctionx/CostFunctionx.cpp
===================================================================
--- /issm/trunk/src/c/CostFunctionx/CostFunctionx.cpp	(revision 3716)
+++ /issm/trunk/src/c/CostFunctionx/CostFunctionx.cpp	(revision 3717)
@@ -31,5 +31,5 @@
 
 	/*add surface area to elements :*/
-	elements->AddInput(S,SurfaceAreaEnum);
+	elements->UpdateInputsFromVector(&S,SurfaceAreaEnum,ConstantEnum);
 	
 	/*Compute gradients: */
Index: /issm/trunk/src/c/DataSet/DataSet.cpp
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3716)
+++ /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3717)
@@ -1413,21 +1413,4 @@
 }
 /*}}}*/
-/*FUNCTION DataSet::AddInput{{{1*/
-void  DataSet::AddInput(double value, int enum_type){
-
-	vector<Object*>::iterator object;
-	Element* element=NULL;
-
-	for ( object=objects.begin() ; object < objects.end(); object++ ){
-
-		if(EnumIsElement((*object)->Enum())){
-
-			element=(Element*)(*object);
-			element->AddInput(value,enum_type);
-		}
-		else ISSMERROR("%s%i%s"," object with id: ",(*object)->Id()," is not an element, in a function that deals only with elements!");
-	}
-}
-/*}}}*/
 /*FUNCTION DataSet::UpdateVertexPositions{{{1*/
 void  DataSet::UpdateVertexPositions(double* thickness,double* bed){
Index: /issm/trunk/src/c/DataSet/DataSet.h
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.h	(revision 3716)
+++ /issm/trunk/src/c/DataSet/DataSet.h	(revision 3717)
@@ -73,5 +73,4 @@
 		void  UpdateInputsFromVector(double* vector,int name, int type);
 		void  UpdateInputsFromSolution(double* solution,int analysis_type,int sub_analysis_type);
-		void  AddInput(double value, int enum_type);
 		void  PenaltyCreateKMatrix(Mat Kgg,double kmax,int analysis_type,int sub_analysis_type);
 		void  PenaltyCreatePVector(Vec pg,double kmax,int analysis_type,int sub_analysis_type);
Index: /issm/trunk/src/c/Dux/Dux.cpp
===================================================================
--- /issm/trunk/src/c/Dux/Dux.cpp	(revision 3716)
+++ /issm/trunk/src/c/Dux/Dux.cpp	(revision 3717)
@@ -33,5 +33,5 @@
 
 	/*add surface area to elements :*/
-	elements->AddInput(S,SurfaceAreaEnum);
+	elements->UpdateInputsFromVector(&S,SurfaceAreaEnum,ConstantEnum);
 	
 	/*Get size of matrix: */
Index: /issm/trunk/src/c/Misfitx/Misfitx.cpp
===================================================================
--- /issm/trunk/src/c/Misfitx/Misfitx.cpp	(revision 3716)
+++ /issm/trunk/src/c/Misfitx/Misfitx.cpp	(revision 3717)
@@ -31,5 +31,5 @@
 
 	/*add surface area to elements :*/
-	elements->AddInput(S,SurfaceAreaEnum);
+	elements->UpdateInputsFromVector(&S,SurfaceAreaEnum,ConstantEnum);
 	
 
Index: /issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp
===================================================================
--- /issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp	(revision 3716)
+++ /issm/trunk/src/c/PenaltyConstraintsx/RiftConstraints.cpp	(revision 3717)
@@ -308,5 +308,5 @@
 
 	/*feed max_penetration to inputs: */
-	loads->AddInput(max_penetration,MaxPenetrationEnum);
+	loads->UpdateInputsFromVector(&max_penetration,MaxPenetrationEnum,ConstantEnum);
 }
 
Index: /issm/trunk/src/c/Qmux/DakotaResponses.cpp
===================================================================
--- /issm/trunk/src/c/Qmux/DakotaResponses.cpp	(revision 3716)
+++ /issm/trunk/src/c/Qmux/DakotaResponses.cpp	(revision 3717)
@@ -230,9 +230,12 @@
 			Vec u_g=NULL;
 			double* u_g_double=NULL;
+			double* vx=NULL;
+			double* vy=NULL;
+			double* vz=NULL;
 			double* fit=NULL;
 
 			/*retrieve active fem model: */
 			model->FindParam(&isstokes,IsStokesEnum);
-			model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum);
+			model->FindParam(&ismacayealpattyn,IsMacAyealPattynEnum);
 			model->FindParam(&ishutter,IsHutterEnum);
 
@@ -258,11 +261,11 @@
 			if(!found)ISSMERROR(" could not find velocity to compute misfit");
 
-			SplitSolutionVector(&vx,&vy,&vz,u_g,numberofnodes*numberofdofspernode);
+			SplitSolutionVector(&vx,&vy,&vz,u_g,numberofnodes,numberofdofspernode);
 
 			/*Add to inputs: */
-			femmodel->elements->AddInput(vx,VxEnum);
-			femmodel->elements->AddInput(vy,VyEnum);
-			femmodel->elements->AddInput(vz,VzEnum);
-			femmodel->elements->AddInput(fit[0],FitEnum);
+			femmodel->elements->UpdateInputsFromVector(vx,VxEnum,VertexEnum);
+			femmodel->elements->UpdateInputsFromVector(vy,VyEnum,VertexEnum);
+			femmodel->elements->UpdateInputsFromVector(vz,VzEnum,VertexEnum);
+			femmodel->elements->UpdateInputsFromVector(fit[0],FitEnum,ConstantEnum);
 
 			/*Compute misfit: */
@@ -300,5 +303,5 @@
 			/*retrieve active fem model: */
 			model->FindParam(&isstokes,IsStokesEnum);
-			model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum);
+			model->FindParam(&ismacayealpattyn,IsMacAyealPattynEnum);
 			model->FindParam(&ishutter,IsHutterEnum);
 
Index: /issm/trunk/src/c/objects/Elements/Beam.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Beam.cpp	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Beam.cpp	(revision 3717)
@@ -749,3 +749,8 @@
 }
 /*}}}1*/
-
+/*FUNCTION Beam::UpdateInputsFromVector(double* vector, int name, int type);{{{1*/
+void  Beam::UpdateInputsFromVector(double* vector, int name, int type){
+	ISSMERROR(" not supported yet!");
+}
+/*}}}*/
+
Index: /issm/trunk/src/c/objects/Elements/Beam.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Beam.h	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Beam.h	(revision 3717)
@@ -53,5 +53,5 @@
 		Object* copy();
 		void  SetClone(int* minranks);
-		void    UpdateInputsFromVector(double* vector, int name, int type){ISSMERROR("Not implemented yet!");}
+		void    UpdateInputsFromVector(double* vector, int name, int type);
 		void    UpdateInputsFromSolution(double* solution, int analysis_type, int sub_analysis_type);
 
@@ -89,5 +89,4 @@
 		void  GetJacobianDeterminant(double* pJdet,double* z_list, double gauss_coord);
 		double MassFlux(double* segment,double* ug);
-		void AddInput(double value, int enum_type){ISSMERROR("not supporte yet!");}
 		/*}}}*/
 
Index: /issm/trunk/src/c/objects/Elements/Element.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Element.h	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Element.h	(revision 3717)
@@ -45,5 +45,4 @@
 		virtual void   ComputeStrainRate(Vec eps,     int analysis_type,int sub_analysis_type)=0;
 		virtual double MassFlux(double* segment,double* ug)=0;
-		virtual void   AddInput(double value, int enum_type)=0;
 
 		/*Implementation: */
Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 3717)
@@ -4458,2 +4458,7 @@
 }
 /*}}}*/
+/*FUNCTION Penta::UpdateInputsFromVector(double* vector, int name, int type);{{{1*/
+void  Penta::UpdateInputsFromVector(double* vector, int name, int type){
+	ISSMERROR(" not supported yet!");
+}
+/*}}}*/
Index: /issm/trunk/src/c/objects/Elements/Penta.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.h	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Penta.h	(revision 3717)
@@ -141,5 +141,4 @@
 		void  UpdateInputsFromSolutionBalancedvelocities( double* solution,int analysis_type,int sub_analysis_type);
 		void  UpdateInputsFromVector(double* vector, int name, int type);
-		void  AddInput(double value, int enum_type){ISSMERROR("not supported yet!");}
 
 		/*}}}*/
Index: /issm/trunk/src/c/objects/Elements/Sing.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Sing.cpp	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Sing.cpp	(revision 3717)
@@ -544,2 +544,7 @@
 }
 /*}}}1*/
+/*FUNCTION Sing::UpdateInputsFromVector(double* vector, int name, int type);{{{1*/
+void  Sing::UpdateInputsFromVector(double* vector, int name, int type){
+	ISSMERROR(" not supported yet!");
+}
+/*}}}*/
Index: /issm/trunk/src/c/objects/Elements/Sing.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Sing.h	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Sing.h	(revision 3717)
@@ -53,6 +53,6 @@
 		void  SetClone(int* minranks);
 		int   MyRank();
-		void    UpdateInputsFromVector(double* vector, int name, int type){ISSMERROR("Not implemented yet!");}
-		void    UpdateInputsFromSolution(double* solution, int analysis_type, int sub_analysis_type);
+		void  UpdateInputsFromVector(double* vector, int name, int type);
+		void  UpdateInputsFromSolution(double* solution, int analysis_type, int sub_analysis_type);
 
 		/*}}}*/
@@ -84,5 +84,4 @@
 		double CostFunction(int,int);
 		double MassFlux(double* segment,double* ug);
-		void AddInput(double value, int enum_type){ISSMERROR("not supporte yet!");}
 		/*}}}*/
 
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 3717)
@@ -5149,2 +5149,7 @@
 }
 /*}}}*/
+/*FUNCTION Tria::UpdateInputsFromVector(double* vector, int name, int type);{{{1*/
+void  Tria::UpdateInputsFromVector(double* vector, int name, int type){
+	ISSMERROR(" not supported yet!");
+}
+/*}}}*/
Index: /issm/trunk/src/c/objects/Elements/Tria.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.h	(revision 3716)
+++ /issm/trunk/src/c/objects/Elements/Tria.h	(revision 3717)
@@ -50,5 +50,5 @@
 		int   MyRank();
 		void  SetClone(int* minranks);
-		void    UpdateInputsFromVector(double* vector, int name, int type){ISSMERROR("Not implemented yet!");}
+		void    UpdateInputsFromVector(double* vector, int name, int type);
 		void    UpdateInputsFromSolution(double* solution, int analysis_type, int sub_analysis_type);
 
@@ -127,5 +127,4 @@
 		void  UpdateInputsFromSolutionBalancedthickness2( double* solution,int analysis_type,int sub_analysis_type);
 		void  UpdateInputsFromSolutionBalancedvelocities( double* solution,int analysis_type,int sub_analysis_type);
-		void  AddInput(double value, int enum_type){ISSMERROR("not supporte yet!");}
 
 		/*}}}*/
Index: /issm/trunk/src/c/parallel/ControlInitialization.cpp
===================================================================
--- /issm/trunk/src/c/parallel/ControlInitialization.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/ControlInitialization.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void ControlInitialization(Model* model, ParameterInputs* inputs){
@@ -54,5 +54,5 @@
 	model->FindParam(&dim,DimEnum);
 	model->FindParam(&ishutter,IsHutterEnum);
-	model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum);
+	model->FindParam(&ismacayealpattyn,IsMacAyealPattynEnum);
 	model->FindParam(&isstokes,IsStokesEnum);
 	model->FindParam(&numberofnodes,NumberOfNodesEnum);
Index: /issm/trunk/src/c/parallel/ControlRestart.cpp
===================================================================
--- /issm/trunk/src/c/parallel/ControlRestart.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/ControlRestart.cpp	(revision 3717)
@@ -5,5 +5,5 @@
 
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void ControlRestart(Model* model,double* param_g){
Index: /issm/trunk/src/c/parallel/ProcessResults.cpp
===================================================================
--- /issm/trunk/src/c/parallel/ProcessResults.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/ProcessResults.cpp	(revision 3717)
@@ -115,5 +115,5 @@
 	model->FindParam(&ishutter,IsHutterEnum);
 	model->FindParam(&isstokes,IsStokesEnum);
-	model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum);
+	model->FindParam(&ismacayealpattyn,IsMacAyealPattynEnum);
 
 	/*Recover femmodels first: */
Index: /issm/trunk/src/c/parallel/control_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/control_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/control_core.cpp	(revision 3717)
@@ -4,5 +4,5 @@
 
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void control_core(DataSet* results,Model* model, ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/convergence.cpp
===================================================================
--- /issm/trunk/src/c/parallel/convergence.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/convergence.cpp	(revision 3717)
@@ -4,5 +4,5 @@
 
 #include "../objects/objects.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void convergence(int* pconverged, Mat Kff,Vec pf,Vec uf,Vec old_uf,Parameters* parameters){
Index: /issm/trunk/src/c/parallel/diagnostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void diagnostic_core(DataSet* results,Model* model){
@@ -62,5 +62,5 @@
 	model->FindParam(&dim,DimEnum);
 	model->FindParam(&ishutter,IsHutterEnum);
-	model->FindParam(&ismacayealpattyn,IsMacayealPattynEnum);
+	model->FindParam(&ismacayealpattyn,IsMacAyealPattynEnum);
 	model->FindParam(&numberofnodes,NumberOfNodesEnum);
 	model->FindParam(&isstokes,IsStokesEnum);
Index: /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 3717)
@@ -6,5 +6,5 @@
 #include "../objects/objects.h"
 #include "../EnumDefinitions/EnumDefinitions.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void diagnostic_core_linear(Vec* pug,FemModel* fem,ParameterInputs* inputs,int analysis_type,int sub_analysis_type){
Index: /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/diagnostic_core_nonlinear.cpp	(revision 3717)
@@ -6,5 +6,5 @@
 #include "../objects/objects.h"
 #include "../EnumDefinitions/EnumDefinitions.h"
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/gradjcompute_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/gradjcompute_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/gradjcompute_core.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/objectivefunctionC.cpp
===================================================================
--- /issm/trunk/src/c/parallel/objectivefunctionC.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/objectivefunctionC.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/prognostic.cpp
===================================================================
--- /issm/trunk/src/c/parallel/prognostic.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/prognostic.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/prognostic2.cpp
===================================================================
--- /issm/trunk/src/c/parallel/prognostic2.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/prognostic2.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/prognostic2_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/prognostic2_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/prognostic2_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void prognostic2_core(DataSet* results,Model* model,ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/prognostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/prognostic_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/prognostic_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void prognostic_core(DataSet* results,Model* model,ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/slopecompute.cpp
===================================================================
--- /issm/trunk/src/c/parallel/slopecompute.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/slopecompute.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/slopecompute_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/slopecompute_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/slopecompute_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void slopecompute_core(DataSet* results,Model* model,ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/steadystate.cpp
===================================================================
--- /issm/trunk/src/c/parallel/steadystate.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/steadystate.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/steadystate_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/steadystate_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/steadystate_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void steadystate_core(DataSet* results,Model* model, ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/thermal.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/thermal.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/thermal_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/thermal_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void thermal_core(DataSet* results,Model* model, ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/thermal_core_nonlinear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal_core_nonlinear.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/thermal_core_nonlinear.cpp	(revision 3717)
@@ -6,5 +6,5 @@
 #include "../objects/objects.h"
 #include "../EnumDefinitions/EnumDefinitions.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void thermal_core_nonlinear(Vec* ptg,double* pmelting_offset,FemModel* fem,ParameterInputs* inputs,int analysis_type,int sub_analysis_type){
Index: /issm/trunk/src/c/parallel/transient.cpp
===================================================================
--- /issm/trunk/src/c/parallel/transient.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/transient.cpp	(revision 3717)
@@ -3,5 +3,5 @@
  */ 
 
-#include "../issm.h"
+#include "../modules.h"
 #include "./parallel.h"
 
Index: /issm/trunk/src/c/parallel/transient_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/transient_core.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/transient_core.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void transient_core(DataSet* results,Model* model, ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/transient_core_2d.cpp
===================================================================
--- /issm/trunk/src/c/parallel/transient_core_2d.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/transient_core_2d.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void transient_core_2d(DataSet* results,Model* model, ParameterInputs* inputs){
Index: /issm/trunk/src/c/parallel/transient_core_3d.cpp
===================================================================
--- /issm/trunk/src/c/parallel/transient_core_3d.cpp	(revision 3716)
+++ /issm/trunk/src/c/parallel/transient_core_3d.cpp	(revision 3717)
@@ -8,5 +8,5 @@
 #include "../EnumDefinitions/EnumDefinitions.h"
 #include "./parallel.h"
-#include "../issm.h"
+#include "../modules.h"
 
 void transient_core_3d(DataSet* results,Model* model, ParameterInputs* inputs){