Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 26323)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 26324)
@@ -553,5 +553,5 @@
 
 					if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0.){
-						vec_constraint_nodes->SetValue(node->Sid(),1.0,INS_VAL);
+						vec_constraint_nodes->SetValue(node->Pid(),1.0,INS_VAL);
 					}
 				}
@@ -582,5 +582,5 @@
 				for(int in=0;in<numnodes;in++){
 					Node* node=element->GetNode(in);
-					if(constraint_nodes[node->Sid()]==1.){
+					if(constraint_nodes[node->Pid()]==1.){
 						isconnected = true;
 						break;
@@ -600,7 +600,7 @@
 						surface_input->GetInputValue(&surface,gauss);
 
-						if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0. && levelset>-300. && levelset<0. && constraint_nodes[node->Sid()]==0.){
+						if((surface_crevasse-surface>0. || crevassedepth-thickness>0.) && bed<0. && levelset>-300. && levelset<0. && constraint_nodes[node->Pid()]==0.){
 							local_nflipped++;
-							vec_constraint_nodes->SetValue(node->Sid(),1.0,INS_VAL);
+							vec_constraint_nodes->SetValue(node->Pid(),1.0,INS_VAL);
 						}
 					}
@@ -631,5 +631,5 @@
 				if(!node->IsActive()) continue;
 
-				if(constraint_nodes[node->Sid()]>0.){
+				if(constraint_nodes[node->Pid()]>0.){
 					node->ApplyConstraint(0,+1.);
 				}