Changeset 25778

Timestamp:

11/23/20 15:54:26 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: limit channel cross section

File:

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

issm/trunk-jpl/src/c/classes/Loads/Channel.cpp

@@ -702 +702 @@
                 /*Solve ODE*/
                 this->S = ODE1(alpha,beta,this->Sold,dt,2);
-                if(this->S<0.) this->S = 0.;
                 _assert_(!xIsNan<IssmDouble>(this->S)); 
+
+                /*Constrain the cross section to be between 0 and 500 m^2*/
+                if(this->S<0.)   this->S = 0.;
+                if(this->S>500.) this->S = 500.;
 
                 count++;