Index: ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27014)
+++ ../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 27015)
@@ -681,7 +681,7 @@
 				b_input->GetInputValue(&bed,gauss);
 				sl_input->GetInputValue(&sealevel,gauss);
 				dis_input->GetInputValue(&distance,gauss);
-				if(thickness<min_thickness && bed<sealevel && distance<mig_max*dt){
+				if(thickness<min_thickness && bed<sealevel && fabs(distance)<mig_max*dt){
 					node->ApplyConstraint(0,+1.);
 				}
 				else {
@@ -719,7 +719,7 @@
 				b_input->GetInputValue(&bed,gauss);
 				sl_input->GetInputValue(&sealevel,gauss);
 				dis_input->GetInputValue(&distance,gauss);
-				if(thickness<min_thickness && bed<sealevel && distance<mig_max*dt){
+				if(thickness<min_thickness && bed<sealevel && fabs(distance)<mig_max*dt){
 					node->ApplyConstraint(0,+1.);
 				}
 				else {
@@ -760,7 +760,7 @@
 				dis_input->GetInputValue(&distance,gauss);
 				hab_fraction_input->GetInputValue(&hab_fraction,gauss);
 
-				if(thickness<((rho_water/rho_ice)*(1+hab_fraction)*-water_depth) && distance<mig_max*dt){
+				if(thickness<((rho_water/rho_ice)*(1+hab_fraction)*-water_depth) && fabs(distance)<mig_max*dt){
 					node->ApplyConstraint(0,+1.);
 				}
 				else {
@@ -862,7 +862,7 @@
 						thickness_input->GetInputValue(&thickness,gauss);
 						surface_input->GetInputValue(&surface,gauss);
 
-                  /*FIXME: not sure about levelset<0. && levelset>-mig_max*dt! SHould maybe be distance<mig_max*dt*/
+                  /*FIXME: not sure about levelset<0. && fabs(levelset)>-mig_max*dt! SHould maybe be distance<mig_max*dt*/
                   if((surface_crevasse>surface || crevassedepth>crevasse_threshold*thickness) && bed<0. && levelset<0. && levelset>-mig_max*dt && constraint_nodes[node->Lid()]==0.){
 							local_nflipped++;
 							vec_constraint_nodes->SetValue(node->Pid(),1.0,INS_VAL);