Changeset 27538

Timestamp:

01/24/23 08:28:00 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added perturbation to linear melt calculation

File:

• issm/trunk-jpl/src/c/classes/Elements/Element.cpp (modified) (3 diffs)

issm/trunk-jpl/src/c/classes/Elements/Element.cpp

@@ -2535 +2535 @@
         IssmDouble deepwaterel,upperwaterel,deepwatermelt,upperwatermelt;
         IssmDouble base[MAXVERTICES];
+        IssmDouble perturbation[MAXVERTICES];
         IssmDouble values[MAXVERTICES];
         IssmDouble time;
@@ -2546 +2547 @@
 
         this->GetInputListOnVertices(&base[0],BaseEnum);
+        this->GetInputListOnVertices(&perturbation[0],BasalforcingsPerturbationMeltingRateEnum);
         for(int i=0;i<NUM_VERTICES;i++){
                 if(base[i]>=upperwaterel){
@@ -2557 +2559 @@
                         values[i]=deepwatermelt*alpha+(1.-alpha)*upperwatermelt;
                 }
+
+                /*Add perturbation*/
+                values[i] += perturbation[i];
         }