Index: /issm/trunk-jpl/src/c/classes/Contours.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Contours.h	(revision 24107)
+++ /issm/trunk-jpl/src/c/classes/Contours.h	(revision 24108)
@@ -22,11 +22,8 @@
 
 	/*intermediary: */
-	int                  nprof;
-	int                 *profnvertices = NULL;
-	doubletype         **pprofx        = NULL;
-	doubletype         **pprofy        = NULL;
-
-	/*output: */
-	Contours *domain = NULL;
+	int          nprof;
+	int         *profnvertices = NULL;
+	doubletype **pprofx        = NULL;
+	doubletype **pprofy        = NULL;
 
 	/*If domainname is an empty string, return empty dataset*/
@@ -35,9 +32,9 @@
 	}
 	else{
-		ExpRead(&nprof,&profnvertices,&pprofx, &pprofy, NULL,domainname);
+		ExpRead<doubletype>(&nprof,&profnvertices,&pprofx, &pprofy, NULL,domainname);
 	}
 
 	/*now create dataset of contours: */
-	domain=new Contours();
+	Contours *domain=new Contours();
 
 	for(int i=0;i<nprof;i++){
Index: /issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp	(revision 24107)
+++ /issm/trunk-jpl/src/c/modules/SystemMatricesx/SystemMatricesx.cpp	(revision 24108)
@@ -120,4 +120,25 @@
 	//Kff->AllocationInfo();
 	//Kfs->AllocationInfo();
+	if(analysisenum==GLheightadvectionAnalysisEnum){
+		PetscViewer viewer;      PetscViewerASCIIOpen(PETSC_COMM_WORLD,"KFF.m",&viewer);      PetscViewerPushFormat(viewer,PETSC_VIEWER_ASCII_MATLAB);
+		MatView(Kff->pmatrix->matrix,viewer);
+		PetscViewerPopFormat(viewer);
+		PetscViewerDestroy(&viewer);
+	for(int e=0;e<femmodel->elements->Size();e++){
+		Element* element=xDynamicCast<Element*>(femmodel->elements->GetObjectByOffset(e));
+		int numdofs;
+		for(int i=0;i<element->GetNumberOfNodes();i++){
+			numdofs=element->nodes[i]->GetNumberOfDofs(NoneApproximationEnum,FsetEnum);
+			int  *doflist = xNew<int>(numdofs);
+			element->nodes[i]->GetDofList(doflist,NoneApproximationEnum,FsetEnum);
+			for(int j=0;j<numdofs;j++){
+				if(doflist[j]==11796-1){
+					printf("node: %i, element: %i, element->nodes[i]->Sid: %i \n", i, element->id, element->nodes[i]->Sid());
+				}
+			}
+			xDelete<int>(doflist);
+		}
+	}
+	}
 
 	/*cleanu up and assign output pointers: */
Index: /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_glads_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_glads_nonlinear.cpp	(revision 24107)
+++ /issm/trunk-jpl/src/c/solutionsequences/solutionsequence_glads_nonlinear.cpp	(revision 24108)
@@ -77,5 +77,5 @@
 
 		/*Converged if inner loop converged in one solution*/
-		//if(count_in==1) converged_out = true;
+		if(count_in==1) converged_out = true;
 
 		/*Increase count: */