Changeset 22526

Timestamp:

03/12/18 08:30:55 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: better handling of anisotropic diffusion

File:

• issm/trunk-jpl/src/c/analyses/Balancethickness2Analysis.cpp (modified) (2 diffs)

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

@@ -105 +105 @@
                         vel = sqrt(vx*vx+vy*vy);
                 }
-                else if(vel<0.1/yts){
-                        vx = vx/vel*0.1;
-                        vy = vy/vel*0.1;
+                else if(vel<30./yts){
+                        vx = 0.;//vx/vel*0.01;
+                        vy = 0.;//vy/vel*0.01;
                         vel = sqrt(vx*vx+vy*vy);
+                        vel = 30./yts*500000.;
                 }
 
@@ -128 +129 @@
 ElementVector* Balancethickness2Analysis::CreatePVector(Element* element){/*{{{*/
 
+        return NULL;
         /*Intermediaries */
         IssmDouble  dhdt[2],mb[2],ms[2],Jdet;