source: issm/oecreview/Archive/17984-18295/ISSM-18190-18191.diff

../trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp

@@ -2941 +2941 @@
 
         }
 
+        /*Transform Coordinate System*/
+        element->TransformStiffnessMatrixCoord(Ke,cs_list);
+
+        /*The pressure augmentation should not be transformed*/
         MatrixMultiply(BtBUzawa,pnumdof,numdof,1,
                                 BtBUzawa,pnumdof,numdof,0,
                                 &Ke->values[0],1);
 
-        /*Transform Coordinate System*/
-        element->TransformStiffnessMatrixCoord(Ke,cs_list);
-
         /*Clean up and return*/
         delete gauss;
         xDelete<IssmDouble>(xyz_list);
@@ -3722 +3723 @@
         }
 
         /*Transform coordinate system*/
-        element->TransformLoadVectorCoord(pe,cs_list);
+        //element->TransformLoadVectorCoord(pe,cs_list); Do not transform pressure augmentation
 
         /*Clean up and return*/
         delete gauss;