Index: /issm/trunk/src/c/DataSet/DataSet.cpp
===================================================================
--- /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3495)
+++ /issm/trunk/src/c/DataSet/DataSet.cpp	(revision 3496)
@@ -902,10 +902,6 @@
 	for ( object=objects.begin() ; object < objects.end(); object++ ){
 
-		/*Check this is a object: */
-		if((*object)->Enum()==NodeEnum()){
-
-			dofobject=dynamic_cast<DofObject*>(*object);
-			dofobject->CreatePartition(partition);
-		}
+		dofobject=dynamic_cast<DofObject*>(*object);
+		dofobject->CreatePartition(partition);
 	}
 
Index: /issm/trunk/src/c/objects/FemModel.cpp
===================================================================
--- /issm/trunk/src/c/objects/FemModel.cpp	(revision 3495)
+++ /issm/trunk/src/c/objects/FemModel.cpp	(revision 3496)
@@ -14,5 +14,6 @@
 #include "./FemModel.h"
 
-/*Object constructors, destructors: {{{1:*/
+/*constructors/destructors*/
+/*FUNCTION FemModel::FemModel {{{1*/
 FemModel::FemModel(){
 
@@ -35,5 +36,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FemModel {{{1*/
 FemModel::FemModel(DataSet* femmodel_elements,DataSet* femmodel_nodes,DataSet* femmodel_vertices, DataSet* femmodel_constraints,DataSet* femmodel_loads,
 		DataSet* femmodel_materials,DataSet* femmodel_parameters, DofVec* femmodel_partition,DofVec* femmodel_tpartition,DofVec* femmodel_yg,
@@ -59,5 +61,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::~FemModel {{{1*/
 FemModel::~FemModel(){
 
@@ -81,5 +84,7 @@
 }
 /*}}}*/
-/* Object management: {{{1*/
+
+/*object functions*/
+/*FUNCTION FemModel::Echo {{{1*/
 void FemModel::Echo(void){
 
@@ -102,5 +107,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::DeepEcho {{{1*/
 void FemModel::DeepEcho(void){
 	
@@ -137,33 +143,46 @@
 
 }
+/*}}}*/
+/*FUNCTION FemModel::GetId {{{1*/
 int   FemModel::GetId(void){
 	ISSMERROR("not implemented yet!");
 }
-
+/*}}}*/
+/*FUNCTION FemModel::MyRank {{{1*/
 int   FemModel::MyRank(void){
 	ISSMERROR("not implemented yet!");
 }
-
+/*}}}*/
+/*FUNCTION FemModel::Marshall {{{1*/
 void  FemModel::Marshall(char** pmarshalled_dataset){
 	ISSMERROR("not implemented yet!");
 }
-
+/*}}}*/
+/*FUNCTION FemModel::MarshallSize {{{1*/
 int   FemModel::MarshallSize(void){
 	ISSMERROR("not implemented yet!");
 }
+/*}}}*/
+/*FUNCTION FemModel::GetName {{{1*/
 char*   FemModel::GetName(void){
 	ISSMERROR("not implemented yet!");
 }
+/*}}}*/
+/*FUNCTION FemModel::Demarshall {{{1*/
 void  FemModel::Demarshall(char** pmarshalled_dataset){
 	ISSMERROR("not implemented yet!");
 }
+/*}}}*/
+/*FUNCTION FemModel::Enum {{{1*/
 int   FemModel::Enum(void){
 	ISSMERROR("not implemented yet!");
 }
+/*}}}*/
+/*FUNCTION FemModel::copy {{{1*/
 Object* FemModel::copy(void){
 	ISSMERROR("not implemented yet!");
 }
-
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(double* pscalar,char* name) {{{1*/
 int FemModel::FindParam(double* pscalar,char* name){
 	
@@ -171,4 +190,6 @@
 
 }
+/*}}}*/
+/*FUNCTION FemModel::FindParam(int* pinteger,char* name) {{{1*/
 int FemModel::FindParam(int* pinteger,char* name){
 	
@@ -176,5 +197,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(char** pstring,char* name) {{{1*/
 int FemModel::FindParam(char** pstring,char* name){
 	
@@ -182,5 +204,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(char*** pstringarray,int* pM,char* name) {{{1*/
 int FemModel::FindParam(char*** pstringarray,int* pM,char* name){
 	
@@ -188,5 +211,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(double** pdoublearray,int* pM,int* pN,char* name) {{{1*/
 int FemModel::FindParam(double** pdoublearray,int* pM,int* pN,char* name){
 	
@@ -194,5 +218,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(Vec* pvec,char* name) {{{1*/
 int FemModel::FindParam(Vec* pvec,char* name){
 	
@@ -200,5 +225,6 @@
 
 }
-
+/*}}}*/
+/*FUNCTION FemModel::FindParam(Mat* pmat,char* name) {{{1*/
 int FemModel::FindParam(Mat* pmat,char* name){
 	
@@ -206,21 +232,80 @@
 
 }
-
+/*}}}*/
 
 /*access to internal data: */
-DataSet*            FemModel::get_elements(void){return elements;}
-DataSet*            FemModel::get_nodes(void){return nodes ;}
-DataSet*            FemModel::get_vertices(void){return vertices ;}
-DataSet*            FemModel::get_constraints(void){return constraints ;}
-DataSet*            FemModel::get_loads(void){return loads;}
-DataSet*            FemModel::get_materials(void){return materials;}
-DataSet*            FemModel::get_parameters(void){return parameters;}
-DofVec*                 FemModel::get_partition(void){return partition;}
-DofVec*                 FemModel::get_tpartition(void){return tpartition;}
-DofVec*                 FemModel::get_yg(void){return yg;}
-Mat                 FemModel::get_Rmg(void){return Rmg;}
-NodeSets*           FemModel::get_nodesets(void){return nodesets;}
-Vec                 FemModel::get_ys(void){return ys;}
-Vec                 FemModel::get_ys0(void){return ys0;}
-Mat                 FemModel::get_Gmn(void){return Gmn;}
-/*}}}*/
+/*FUNCTION FemModel::get_elements {{{1*/
+DataSet* FemModel::get_elements(void){
+	return elements;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_nodes {{{1*/
+DataSet* FemModel::get_nodes(void){
+	return nodes;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_vertices {{{1*/
+DataSet* FemModel::get_vertices(void){
+	return vertices;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_constraints {{{1*/
+DataSet* FemModel::get_constraints(void){
+	return constraints;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_loads {{{1*/
+DataSet* FemModel::get_loads(void){
+	return loads;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_materials {{{1*/
+DataSet* FemModel::get_materials(void){
+	return materials;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_parameters {{{1*/
+DataSet* FemModel::get_parameters(void){
+	return parameters;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_partition {{{1*/
+DofVec* FemModel::get_partition(void){
+	return partition;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_tpartition {{{1*/
+DofVec* FemModel::get_tpartition(void){
+	return tpartition;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_yg {{{1*/
+DofVec* FemModel::get_yg(void){
+	return yg;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_Rmg {{{1*/
+Mat FemModel::get_Rmg(void){
+	return Rmg;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_nodesets {{{1*/
+NodeSets* FemModel::get_nodesets(void){
+	return nodesets;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_ys {{{1*/
+Vec FemModel::get_ys(void){
+	return ys;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_ys0 {{{1*/
+Vec FemModel::get_ys0(void){
+	return ys0;
+}
+/*}}}*/
+/*FUNCTION FemModel::get_Gmn {{{1*/
+Mat FemModel::get_Gmn(void){
+	return Gmn;
+}
+/*}}}*/
Index: /issm/trunk/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk/src/c/shared/Elements/elements.h	(revision 3495)
+++ /issm/trunk/src/c/shared/Elements/elements.h	(revision 3496)
@@ -12,5 +12,16 @@
 double Paterson(double temperature);
 int GetVerticesCoordinates(double* xyz,  Node** nodes, int numgrids);
+inline void printarray(double* array,int lines,int cols){
+	int i,j;
+	printf("\n");
+	for (i=0;i<lines;i++){  
+		printf("   [ ");
+		for (j=0;j<cols;j++){
+			printf(" %10g ",array[i*cols+j]);
+		}  
+		printf(" ]\n");
+	}  
+	printf("\n");
+}
 
 #endif //ifndef _SHARED_ELEMENTS_H_
-
