Index: /issm/trunk/src/m/solutions/ResetBoundaryConditions.m
===================================================================
--- /issm/trunk/src/m/solutions/ResetBoundaryConditions.m	(revision 4442)
+++ /issm/trunk/src/m/solutions/ResetBoundaryConditions.m	(revision 4442)
@@ -0,0 +1,32 @@
+function femmodel=ResetBoundaryConditions(femmodel,analyses_type,sub_analyses_type);
+%RESETBOUNDARYCONDITIONS - change boundary conditions of a model
+%
+%   Change boundary condition of a model using a vector solution from another model
+%
+%   Usage:
+%      femmodel=ResetBoundaryConditions(femmodel,analyses_type,sub_analyses_type)
+
+
+	%recover parameters common to all solutions
+	verbose=femmodel.parameters.Verbose;
+
+	displaystring(verbose,'\n%s',['updating boundary condition ...']);
+
+	ug=GetSolutionFromInputs(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
+	%set current analysis: 
+	femmodel=SetCurrentConfiguration(femmodel,analysis_type);
+
+	%For this analysis_type, free existing boundary condition vectors:
+	analysis_counter=femmodel.AnalysisCounter;
+
+	%global dof set
+	femmodel.m_yg{analysis_counter}=[];
+	%in the s-set
+	femmodel.m_ys{analysis_counter}=[];
+
+	%Now, duplicate ug (the solution vector) into the boundary conditions vector on the g-set
+	femmodel.m_yg{analysis_counter}=ug;
+
+	%Reduce from g to s set
+	femmodel.m_ys{analysis_counter}=Reducevectorgtos(femmodel.m_yg{analysis_counter},femmodel.m_nodesets{analysis_counter});