Changeset 23827

Timestamp:

04/07/19 09:12:23 (6 years ago)

Author:

Mathieu Morlighem

Message:

CHG: calve only if bed<sealevel

File:

• issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp (modified) (3 diffs)

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

@@ -699 +699 @@
         IssmDouble  rho_ice,rho_water;
         IssmDouble  bed,water_depth;
-        IssmDouble  levelset,grounded_ls;
+        IssmDouble  levelset,sealevel;
 
         femmodel->parameters->FindParam(&calvinglaw,CalvingLawEnum);
@@ -715 +715 @@
                         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 +725 @@
 
                                 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.);
                                 }