Index: /issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.cpp	(revision 15861)
+++ /issm/trunk-jpl/src/c/classes/Constraints/SpcTransient.cpp	(revision 15862)
@@ -101,10 +101,10 @@
 void SpcTransient::ConstrainNode(Nodes* nodes,Parameters* parameters){
 
-	Node* node=NULL;
-	IssmDouble time=0;
-	int    i;
-	IssmDouble alpha=-1;
-	IssmDouble value;
-	bool   found=false;
+	Node       *node  = NULL;
+	IssmDouble  time  = 0.;
+	int         i;
+	IssmDouble  alpha = -1.;
+	IssmDouble  value;
+	bool        found = false;
 
 	/*Chase through nodes and find the node to which this SpcTransient applys: */
@@ -140,4 +140,5 @@
 		/*Apply or relax constraint: */
 		if(xIsNan<IssmDouble>(value)){
+			printf("-------------- file: SpcTransient.cpp line: %i\n",__LINE__); 
 			node->RelaxConstraint(dof);
 		}
Index: /issm/trunk-jpl/src/c/classes/Node.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15861)
+++ /issm/trunk-jpl/src/c/classes/Node.cpp	(revision 15862)
@@ -393,5 +393,5 @@
 	/*Dof should be added to the f-set, and taken out of the s-set:*/
 	DofInFSet(dof-1);
-	this->indexing.svalues[dof-1]=NAN;
+	this->indexing.svalues[dof-1]=0.;
 }
 /*}}}*/
@@ -454,6 +454,9 @@
 	_assert_(dof<this->indexing.gsize);
 
-	this->indexing.f_set[dof]=0; //n splits into f (for which we solve) and s (single point constraints)
-	this->indexing.s_set[dof]=1;
+	if(this->indexing.f_set[dof] == 1){
+		this->indexingupdate = true;
+		this->indexing.f_set[dof]=0; //n splits into f (for which we solve) and s (single point constraints)
+		this->indexing.s_set[dof]=1;
+	}
 }
 /*}}}*/
@@ -464,6 +467,9 @@
 	 * to a fixed value during computations. */
 
-	this->indexing.f_set[dof]=1; 
-	this->indexing.s_set[dof]=0;
+	if(this->indexing.f_set[dof] == 0){
+		this->indexingupdate = true;
+		this->indexing.f_set[dof]=1; 
+		this->indexing.s_set[dof]=0;
+	}
 }
 /*}}}*/
