Changeset 18784

Timestamp:

11/14/14 15:06:29 (10 years ago)

Author:

srebuffi

Message:

CHG : correct stabilization for calving

File:

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

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

@@ -272 +272 @@
                 /* Stabilization */
                 vel=0.;
-                for(i=0;i<dim;i++) vel+=v[i]*v[i];
+                for(i=0;i<dim;i++) vel+=w[i]*w[i];
                 vel=sqrt(vel)+1.e-14;
                 switch(stabilization){
@@ -298 +298 @@
                                 /* Streamline Upwinding */
                                 basalelement->ElementSizes(&hx,&hy,&hz);
-                                h=sqrt( pow(hx*v[0]/vel,2) + pow(hy*v[1]/vel,2) );
+                                h=sqrt( pow(hx*w[0]/vel,2) + pow(hy*w[1]/vel,2) );
                                 for(row=0;row<dim;row++) 
                                         for(col=0;col<dim;col++) 
-                                                D[row*dim+col] = D_scalar*h/(2.*vel)*v[row]*v[col];
+                                                D[row*dim+col] = D_scalar*h/(2.*vel)*w[row]*w[col];
 
                                 TripleMultiply(Bprime,dim,numnodes,1,