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

Timestamp:

08/06/15 13:31:32 (10 years ago)

Author:

seroussi

Message:

BUG: fixed mismip melting rate and added NR test

File:

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

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

@@ -1664 +1664 @@
         this->GetInputListOnVertices(bed,BedEnum);
         for(int i=0;i<numvertices;i++){
-                if(base[i]<upperdepthmelt)  values[i]=0;
-                else values[i]=meltratefactor*(tanh(base[i]-bed[i])/thresholdthickness)*(upperdepthmelt-bed[i]);
+                if(base[i]>upperdepthmelt){
+                        values[i]=0;
+                }
+                else{
+                        values[i]=meltratefactor*tanh((base[i]-bed[i])/thresholdthickness)*(upperdepthmelt-base[i]);
+                }
         }