Changeset 23037

Timestamp:

08/01/18 11:13:54 (7 years ago)

Author:

tpelle

Message:

CHG: Modified GL height for basal melt plume parameterization

File:

• issm/trunk-jpl/src/c/analyses/GLheightadvectionAnalysis.cpp (modified) (4 diffs)

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

@@ -87 +87 @@
         //Input* vx_input=element->GetInput(BaseSlopeXEnum); _assert_(vx_input);
         //Input* vy_input=element->GetInput(BaseSlopeYEnum); _assert_(vy_input);
-        Input* vx_input=element->GetInput(VxEnum); _assert_(vx_input);
+   Input* vx_input=element->GetInput(VxEnum); _assert_(vx_input);
         Input* vy_input=element->GetInput(VyEnum); _assert_(vy_input);
 
@@ -105 +105 @@
 
                 /*Diffusion */
-                if(sqrt(vx*vx+vy*vy)<10./31536000.){
-                        IssmPDouble kappa = 1.e-4;
+                if(sqrt(vx*vx+vy*vy)<1000./31536000.){
+                        IssmPDouble kappa = -10;
                         for(int i=0;i<numnodes;i++){
                                 for(int j=0;j<numnodes;j++){
@@ -113 +113 @@
                         }
                 }
-
                 /*Advection: */
                 vx_input->GetInputValue(&vx,gauss);
@@ -119 +118 @@
                 for(int i=0;i<numnodes;i++){
                         for(int j=0;j<numnodes;j++){
-                                Ke->values[i*numnodes+j] += D_scalar*(vx*dbasis[0*numnodes+j]*basis[i] + vy*dbasis[1*numnodes+j]*basis[i]);
+                                Ke->values[i*numnodes+j] += (D_scalar*(vx*dbasis[0*numnodes+j]*basis[i] + vy*dbasis[1*numnodes+j]*basis[i]))*1e-2;
                         }
                 }