Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5908)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 5909)
@@ -3979,94 +3979,4 @@
 }
 /*}}}*/
-/*FUNCTION Penta::GetLocalDofList {{{1*/
-int*  Penta::GetLocalDofList(int approximation_enum,int setenum){
-
-	int  i,j,count,numdof,numgdof;
-	int  ndof_list[NUMVERTICES];
-	int  ngdof_list_cumulative[NUMVERTICES];
-	int *doflist = NULL;
-
-	/*Get number of dofs per node, and total for this given set*/
-	numdof=0;
-	numgdof=0;
-	for(i=0;i<NUMVERTICES;i++){
-
-		/*Cumulative list= number of dofs before node i*/
-		ngdof_list_cumulative[i]=numgdof;
-
-		/*Number of dofs for node i for given set and for the g set*/
-		ndof_list[i]=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		numgdof    +=nodes[i]->GetNumberOfDofs(approximation_enum,GsetEnum);
-		numdof     +=ndof_list[i];
-	}
-
-	if(numdof){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numdof*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			nodes[i]->GetLocalDofList(&doflist[count],approximation_enum,setenum);
-			count+=ndof_list[i];
-		}
-
-		/*We now have something like: [0 1 0 2 1 2]. Offset by gsize, to get something like: [0 1 2 4 6 7]:*/
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			for(j=0;j<ndof_list[i];j++){
-				doflist[count+j]+=ngdof_list_cumulative[i];
-			}
-			count+=ndof_list[i];
-		}
-	}
-	else doflist=NULL;
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Penta::GetGlobalDofList {{{1*/
-int* Penta::GetGlobalDofList(int approximation_enum,int setenum){
-
-	int  i,numdof,count;
-	int  ndof_list[NUMVERTICES];
-	int *doflist = NULL;
-
-	/*First, figure out size of doflist: */
-	numdof=0;
-	for(i=0;i<NUMVERTICES;i++){
-		ndof_list[i]=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		numdof+=ndof_list[i];
-	}
-
-	if(numdof){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numdof*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			nodes[i]->GetDofList(&doflist[count],approximation_enum,setenum);
-			count+=ndof_list[i];
-		}
-	}
-	else doflist=NULL;
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Penta::GetNumberOfDofs {{{1*/
-int Penta::GetNumberOfDofs(int approximation_enum,int setenum){
-
-	/*output: */
-	int numberofdofs=0;
-
-	for(int i=0;i<NUMVERTICES;i++){
-		numberofdofs+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-	}
-
-	return numberofdofs;
-}
-/*}}}*/
 /*FUNCTION Penta::GetParameterListOnVertices(double* pvalue,int enumtype) {{{1*/
 void Penta::GetParameterListOnVertices(double* pvalue,int enumtype){
Index: /issm/trunk/src/c/objects/Elements/Penta.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.h	(revision 5908)
+++ /issm/trunk/src/c/objects/Elements/Penta.h	(revision 5909)
@@ -167,7 +167,4 @@
 		void	  GetDofList(int** pdoflist,int approximation_enum,int setenum);
 		void	  GetDofList1(int* doflist);
-		int*    GetLocalDofList(int approximation_enum,int setenum);
-		int*    GetGlobalDofList(int approximation_enum,int setenum);
-		int GetNumberOfDofs(int approximation_enum,int setenum);
 		int     GetElementType(void);
 		void    GetParameterListOnVertices(double* pvalue,int enumtype);
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5908)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 5909)
@@ -4793,94 +4793,4 @@
 }
 /*}}}*/
-/*FUNCTION Tria::GetLocalDofList {{{1*/
-int*  Tria::GetLocalDofList(int approximation_enum,int setenum){
-
-	int  i,j,count,numdof,numgdof;
-	int  ndof_list[NUMVERTICES];
-	int  ngdof_list_cumulative[NUMVERTICES];
-	int *doflist = NULL;
-
-	/*Get number of dofs per node, and total for this given set*/
-	numdof=0;
-	numgdof=0;
-	for(i=0;i<NUMVERTICES;i++){
-
-		/*Cumulative list= number of dofs before node i*/
-		ngdof_list_cumulative[i]=numgdof;
-
-		/*Number of dofs for node i for given set and for the g set*/
-		ndof_list[i]=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		numgdof    +=nodes[i]->GetNumberOfDofs(approximation_enum,GsetEnum);
-		numdof     +=ndof_list[i];
-	}
-
-	if(numdof){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numdof*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			nodes[i]->GetLocalDofList(&doflist[count],approximation_enum,setenum);
-			count+=ndof_list[i];
-		}
-		
-		/*We now have something like: [0 1 0 2 1 2]. Offset by gsize, to get something like: [0 1 2 4 6 7]:*/
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			for(j=0;j<ndof_list[i];j++){
-				doflist[count+j]+=ngdof_list_cumulative[i];
-			}
-			count+=ndof_list[i];
-		}
-	}
-	else doflist=NULL;
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Tria::GetGlobalDofList {{{1*/
-int* Tria::GetGlobalDofList(int approximation_enum,int setenum){
-
-	int  i,numdof,count;
-	int  ndof_list[NUMVERTICES];
-	int *doflist = NULL;
-
-	/*First, figure out size of doflist: */
-	numdof=0;
-	for(i=0;i<NUMVERTICES;i++){
-		ndof_list[i]=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		numdof+=ndof_list[i];
-	}
-
-	if(numdof){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numdof*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<NUMVERTICES;i++){
-			nodes[i]->GetDofList(&doflist[count],approximation_enum,setenum);
-			count+=ndof_list[i];
-		}
-	}
-	else doflist=NULL;
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Tria::GetNumberOfDofs {{{1*/
-int Tria::GetNumberOfDofs(int approximation_enum,int setenum){
-	
-	/*output: */
-	int numberofdofs=0;
-
-	for(int i=0;i<NUMVERTICES;i++){
-		numberofdofs+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-	}
-
-	return numberofdofs;
-}
-/*}}}*/
 /*FUNCTION Tria::GetDofList1 {{{1*/
 void  Tria::GetDofList1(int* doflist){
Index: /issm/trunk/src/c/objects/Elements/Tria.h
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.h	(revision 5908)
+++ /issm/trunk/src/c/objects/Elements/Tria.h	(revision 5909)
@@ -156,7 +156,4 @@
 		void	  GetDofList(int** pdoflist,int approximation_enum,int setenum);
 		void	  GetDofList1(int* doflist);
-		int*    GetLocalDofList(int approximation_enum,int setenum);
-		int*    GetGlobalDofList(int approximation_enum,int setenum);
-		int     GetNumberOfDofs(int approximation_enum,int setenum);
 		void    GetParameterListOnVertices(double* pvalue,int enumtype);
 		void    GetParameterListOnVertices(double* pvalue,int enumtype,double defaultvalue);
Index: /issm/trunk/src/c/objects/Loads/Icefront.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5908)
+++ /issm/trunk/src/c/objects/Loads/Icefront.cpp	(revision 5909)
@@ -742,134 +742,4 @@
 }
 /*}}}*/
-/*FUNCTION Icefront::GetInternalDofList {{{1*/
-int* Icefront::GetInternalDofList(int approximation_enum,int setenum){
-
-	int i,j;
-	int numberofdofs=0;
-	int count=0;
-	int count2=0;
-	int type;
-	int numberofnodes=2;
-	int ndof;
-
-	/*output: */
-	int* doflist=NULL;
-
-	/*recover type: */
-	inputs->GetParameterValue(&type,TypeEnum);
-
-	/*Some checks for debugging*/
-	ISSMASSERT(nodes);
-		
-	/*How many nodes? :*/
-	if(type==MacAyeal2dIceFrontEnum || type==MacAyeal3dIceFrontEnum)
-	 numberofnodes=2;
-	else 
-	 numberofnodes=4;
-	
-	/*Figure out size of doflist: */
-	for(i=0;i<numberofnodes;i++){
-		numberofdofs+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-	}
-
-	if(numberofdofs){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numberofdofs*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<numberofnodes;i++){
-			nodes[i]->GetDofList(doflist+count,approximation_enum,setenum);
-			count+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		}
-		
-		/*We now have something like: [0 1 0 1]. Offset by gsize, to get something like:
-		 * [0 1 2 3]: */
-		count=0;
-		count2=0;
-		for(i=0;i<numberofnodes;i++){
-			ndof=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-			for(j=count;j<count+ndof;j++)doflist[j]+=count2;
-			count+=ndof;
-			count2+=nodes[i]->GetNumberOfDofs(approximation_enum,GsetEnum);
-		}
-	}
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Icefront::GetExternalDofList{{{1*/
-int*  Icefront::GetExternalDofList(int approximation_enum,int setenum){
-
-	int i,j;
-	int numberofdofs=0;
-	int count=0;
-	int type;
-	int numberofnodes=2;
-
-	/*output: */
-	int* doflist=NULL;
-
-	/*recover type: */
-	inputs->GetParameterValue(&type,TypeEnum);
-
-	/*Some checks for debugging*/
-	ISSMASSERT(nodes);
-		
-	/*How many nodes? :*/
-	if(type==MacAyeal2dIceFrontEnum || type==MacAyeal3dIceFrontEnum)
-	 numberofnodes=2;
-	else 
-	 numberofnodes=4;
-	
-	/*Figure out size of doflist: */
-	for(i=0;i<numberofnodes;i++){
-		numberofdofs+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-	}
-
-	if(numberofdofs){
-		/*Allocate: */
-		doflist=(int*)xmalloc(numberofdofs*sizeof(int));
-
-		/*Populate: */
-		count=0;
-		for(i=0;i<numberofnodes;i++){
-			nodes[i]->GetDofList(doflist+count,approximation_enum,setenum);
-			count+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-		}
-	}
-
-	return doflist;
-}
-/*}}}*/
-/*FUNCTION Icefront::GetNumberOfDofs {{{1*/
-int Icefront::GetNumberOfDofs(int approximation_enum,int setenum){
-	
-	
-	int i;
-	int numberofdofs=0;
-	int type;
-	int numberofnodes=2;
-
-	/*recover type: */
-	inputs->GetParameterValue(&type,TypeEnum);
-
-	/*Some checks for debugging*/
-	ISSMASSERT(nodes);
-		
-	/*How many nodes? :*/
-	if(type==MacAyeal2dIceFrontEnum || type==MacAyeal3dIceFrontEnum)
-	 numberofnodes=2;
-	else 
-	 numberofnodes=4;
-	
-	/*Figure out size of doflist: */
-	for(i=0;i<numberofnodes;i++){
-		numberofdofs+=nodes[i]->GetNumberOfDofs(approximation_enum,setenum);
-	}
-
-	return numberofdofs;
-}
-/*}}}*/
 /*FUNCTION Icefront::GetSegmentNormal {{{1*/
 void Icefront:: GetSegmentNormal(double* normal,double xyz_list[4][3]){
Index: /issm/trunk/src/c/objects/Loads/Icefront.h
===================================================================
--- /issm/trunk/src/c/objects/Loads/Icefront.h	(revision 5908)
+++ /issm/trunk/src/c/objects/Loads/Icefront.h	(revision 5909)
@@ -82,7 +82,4 @@
 		void  CreatePVectorDiagnosticStokes(Vec pg);
 		void  GetDofList(int** pdoflist,int approximation_enum,int setenum);
-		int*  GetInternalDofList(int approximation_enum,int setenum);
-		int*  GetExternalDofList(int approximation_enum,int setenum);
-		int   GetNumberOfDofs(int approximation_enum,int setenum);
 		void GetSegmentNormal(double* normal,double xyz_list[2][3]);
 		void GetQuadNormal(double* normal,double xyz_list[4][3]);