Index: /issm/trunk-jpl/src/c/Container/Observations.cpp
===================================================================
--- /issm/trunk-jpl/src/c/Container/Observations.cpp	(revision 14910)
+++ /issm/trunk-jpl/src/c/Container/Observations.cpp	(revision 14911)
@@ -403,7 +403,7 @@
 	/*Solve the three linear systems*/
 #if _HAVE_GSL_
-	SolverxSeq(&GinvG0,Gamma,gamma0,n_obs); // Gamma^-1 gamma0
-	SolverxSeq(&Ginv1, Gamma,ones,n_obs);   // Gamma^-1 ones
-	SolverxSeq(&GinvZ, Gamma,obs,n_obs);    // Gamma^-1 Z
+	DenseGslSolve(&GinvG0,Gamma,gamma0,n_obs); // Gamma^-1 gamma0
+	DenseGslSolve(&Ginv1, Gamma,ones,n_obs);   // Gamma^-1 ones
+	DenseGslSolve(&GinvZ, Gamma,obs,n_obs);    // Gamma^-1 Z
 #else
 	_error_("GSL is required");
@@ -511,5 +511,5 @@
 		/*Compute weights*/
 #if _HAVE_GSL_
-		SolverxSeq(&weights,Green,obs,n_obs); // Green^-1 obs
+		DenseGslSolve(&weights,Green,obs,n_obs); // Green^-1 obs
 #else
 		_error_("GSL is required");
Index: /issm/trunk-jpl/src/c/modules/Shp2Expx/Shp2Expx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/Shp2Expx/Shp2Expx.cpp	(revision 14910)
+++ /issm/trunk-jpl/src/c/modules/Shp2Expx/Shp2Expx.cpp	(revision 14911)
@@ -4,4 +4,6 @@
 
 #include "./Shp2Expx.h"
+#include "../Ll2xyx/Ll2xyx.h"
+#include "../Xy2llx/Xy2llx.h"
 #include "../../shared/shared.h"
 #include "../../include/include.h"
Index: /issm/trunk-jpl/src/c/modules/Shp2Kmlx/Shp2Kmlx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/Shp2Kmlx/Shp2Kmlx.cpp	(revision 14910)
+++ /issm/trunk-jpl/src/c/modules/Shp2Kmlx/Shp2Kmlx.cpp	(revision 14911)
@@ -4,4 +4,5 @@
 
 #include "./Shp2Kmlx.h"
+#include "../Xy2llx/Xy2llx.h"
 #include "../../shared/shared.h"
 #include "../../include/include.h"
Index: /issm/trunk-jpl/src/c/toolkits/gsl/DenseGslSolve.cpp
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/gsl/DenseGslSolve.cpp	(revision 14910)
+++ /issm/trunk-jpl/src/c/toolkits/gsl/DenseGslSolve.cpp	(revision 14911)
@@ -23,4 +23,14 @@
 class Parameters;
 
+void DenseGslSolve(IssmPDouble** pX,IssmPDouble* A,IssmPDouble* B, int n){ /*{{{*/
+
+	/*Intermediary: */
+	IssmPDouble *X  = xNew<IssmPDouble>(n);
+	SolverxSeq(X,A,B,n);
+
+	/*allocate output pointers: */
+	*pX=X;
+}
+/*}}}*/
 void DenseGslSolve(/*output*/ IssmPDouble** px,/*stiffness matrix:*/ IssmPDouble* Kff, int Kff_M, int Kff_N, /*right hand side load vector: */ IssmPDouble* pf, int pf_M, Parameters* parameters){ /*{{{*/
 
Index: /issm/trunk-jpl/src/c/toolkits/gsl/gslincludes.h
===================================================================
--- /issm/trunk-jpl/src/c/toolkits/gsl/gslincludes.h	(revision 14910)
+++ /issm/trunk-jpl/src/c/toolkits/gsl/gslincludes.h	(revision 14911)
@@ -20,4 +20,5 @@
 class Parameters;
 
+void DenseGslSolve(IssmPDouble** pX,IssmPDouble* A,IssmPDouble* B, int n);
 void DenseGslSolve(IssmDouble** px, IssmDouble* Kff,int Kff_M, int Kff_N, IssmDouble* pf, int pf_M, Parameters* parameters);