Index: /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 20152)
+++ /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 20153)
@@ -179,5 +179,34 @@
 }/*}}}*/
 void           SealevelriseAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/
-	_error_("not implemented yet");
+	
+	IssmDouble *deltaS  = NULL;
+	IssmDouble *S  = NULL;
+	int*        sidlist = NULL;
+	int         numvertices;
+	
+	numvertices= element->GetNumberOfVertices();
+	sidlist=xNew<int>(numvertices);
+	
+	element->GetVerticesSidList(sidlist);
+
+	deltaS = xNew<IssmDouble>(numvertices);
+	for(int i=0;i<numvertices;i++){
+		deltaS[i]=solution[sidlist[i]];
+	}
+
+	S = xNew<IssmDouble>(numvertices);
+	element->GetInputListOnVertices(S,SealevelriseSEnum,0);
+
+	/*Add deltaS to S:*/
+	for (int i=0;i<numvertices;i++)S[i]+=deltaS[i];
+
+	/*Add S back into inputs: */
+	element->AddInput(SealevelriseSEnum,S,P1Enum);
+
+	/*Free ressources:*/
+	xDelete<int>(sidlist);
+	xDelete<IssmDouble>(deltaS);
+	xDelete<IssmDouble>(S);
+
 }/*}}}*/
 void           SealevelriseAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/
Index: /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp	(revision 20152)
+++ /issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp	(revision 20153)
@@ -59,5 +59,5 @@
 
 		/*get results into elements:*/
-		InputUpdateFromVectorx(femmodel,Sg,SealevelriseSEnum,VertexSIdEnum);
+		InputUpdateFromSolutionx(femmodel,Sg);
 
 		if(save_results){