source: issm/oecreview/Archive/26740-27031/ISSM-27014-27015.diff@ 27032

../trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp

@@ -681 +681 @@
                                 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 +719 @@
                                 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 +760 @@
                                 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 +862 @@
                                                 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);