Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 23826)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 23827)
@@ -699,5 +699,5 @@
 	IssmDouble  rho_ice,rho_water;
 	IssmDouble  bed,water_depth;
-	IssmDouble  levelset,grounded_ls;
+	IssmDouble  levelset,sealevel;
 
 	femmodel->parameters->FindParam(&calvinglaw,CalvingLawEnum);
@@ -715,5 +715,6 @@
 			Gauss*   gauss    = element->NewGauss();
 			Input*   H_input  = element->GetInput(ThicknessEnum); _assert_(H_input);
-			Input*   gr_input = element->GetInput(MaskGroundediceLevelsetEnum); _assert_(gr_input);
+			Input*   b_input = element->GetInput(BedEnum); _assert_(b_input);
+			Input*   sl_input = element->GetInput(SealevelEnum); _assert_(sl_input);
 
 			/*Potentially constrain nodes of this element*/
@@ -724,6 +725,7 @@
 
 				H_input->GetInputValue(&thickness,gauss);
-				gr_input->GetInputValue(&grounded_ls,gauss);
-				if(thickness<min_thickness && grounded_ls<0.){
+				b_input->GetInputValue(&bed,gauss);
+				sl_input->GetInputValue(&sealevel,gauss);
+				if(thickness<min_thickness && bed<sealevel){
 					node->ApplyConstraint(0,+1.);
 				}