Changeset 27579

Timestamp:

02/15/23 13:29:08 (2 years ago)

Author:

Mathieu Morlighem

Message:

BUG: melt perturbation not implemented yet for depth-dependent melt

File:

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

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

@@ -654 +654 @@
         Input* gmb_input        = element->GetInput(BasalforcingsGroundediceMeltingRateEnum);  _assert_(gmb_input);
         Input* fmb_input        = element->GetInput(BasalforcingsFloatingiceMeltingRateEnum);  _assert_(fmb_input);
-        Input* fmb_pert_input   = element->GetInput(BasalforcingsPerturbationMeltingRateEnum); _assert_(fmb_pert_input);
+        //Input* fmb_pert_input   = element->GetInput(BasalforcingsPerturbationMeltingRateEnum); _assert_(fmb_pert_input);
         Input* gllevelset_input = element->GetInput(MaskOceanLevelsetEnum);              _assert_(gllevelset_input);
         Input* ms_input         = element->GetInput(SmbMassBalanceEnum);                       _assert_(ms_input);
@@ -682 +682 @@
                 gmb_input->GetInputValue(&gmb,gauss);
                 fmb_input->GetInputValue(&fmb,gauss);
-                fmb_pert_input->GetInputValue(&fmb_pert,gauss);
+                //fmb_pert_input->GetInputValue(&fmb_pert,gauss);
                 gllevelset_input->GetInputValue(&gllevelset,gauss);
                 thickness_input->GetInputValue(&thickness,gauss);
@@ -696 +696 @@
                 else if(melt_style==NoMeltOnPartiallyFloatingEnum){
                         if (phi<0.00000001){
-                                mb=fmb+fmb_pert;
+                                mb=fmb;//+fmb_pert;
                         }
                         else mb=gmb;