Changeset 23084

Timestamp:

08/10/18 16:16:47 (7 years ago)

Author:

Mathieu Morlighem

Message:

BUG: gauss point has to be at the base to pick up drag since drag is 0 on all layers except base

File:

• issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp (modified) (1 diff)

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

@@ -475 +475 @@
                 /*Deal with basal velocities*/
                 if(element->IsOnBase()){
+                        /*Make sure to now push the gauss point to the base*/
+                        delete gauss; gauss=element->NewGauss();
+                        gauss->GaussVertex(nodedown);
                         switch(frictionlaw){
                                 case 1: