Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 13908)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 13909)
@@ -342,5 +342,5 @@
 	int  fsize,ssize,flocalsize,slocalsize;
 	int  connectivity, numberofdofspernode;
-	int  analysis_type,configuration_type;
+	int  configuration_type;
 	int  m,n,M,N;
 	int *d_nnz = NULL;
@@ -357,5 +357,4 @@
 
 	/*retrieve parameters: */
-	this->parameters->FindParam(&analysis_type,AnalysisTypeEnum);
 	this->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
 	this->parameters->FindParam(&connectivity,MeshAverageVertexConnectivityEnum);
@@ -364,6 +363,6 @@
 	fsize      = this->nodes->NumberOfDofs(configuration_type,FsetEnum);
 	ssize      = this->nodes->NumberOfDofs(configuration_type,SsetEnum);
-	flocalsize = this->nodes->NumberOfDofsLocal(analysis_type,FsetEnum);
-	slocalsize = this->nodes->NumberOfDofsLocal(analysis_type,SsetEnum);
+	flocalsize = this->nodes->NumberOfDofsLocal(configuration_type,FsetEnum);
+	slocalsize = this->nodes->NumberOfDofsLocal(configuration_type,SsetEnum);
 
 	numberofdofspernode=this->nodes->MaxNumDofs(configuration_type,GsetEnum);
@@ -408,5 +407,5 @@
 	/*Intermediary*/
 	int      i,j,k,index,offset,count;
-	int      analysis_type,configuration_type;
+	int      configuration_type;
 	int      d_nz,o_nz;
 	Element *element      = NULL;
@@ -421,14 +420,13 @@
 
 	/*retrive parameters: */
-	this->parameters->FindParam(&analysis_type,AnalysisTypeEnum);
 	this->parameters->FindParam(&configuration_type,ConfigurationTypeEnum);
 
 	/*Get vector size and number of nodes*/
-	int numnodes            = nodes->NumberOfNodes(analysis_type);
+	int numnodes            = nodes->NumberOfNodes(configuration_type);
 	int numberofdofspernode = nodes->MaxNumDofs(configuration_type,GsetEnum);
-	int M                   = nodes->NumberOfDofs(analysis_type,set1enum);
-	int N                   = nodes->NumberOfDofs(analysis_type,set2enum);
-	int m                   = nodes->NumberOfDofsLocal(analysis_type,set1enum);
-	int n                   = nodes->NumberOfDofsLocal(analysis_type,set2enum);
+	int M                   = nodes->NumberOfDofs(configuration_type,set1enum);
+	int N                   = nodes->NumberOfDofs(configuration_type,set2enum);
+	int m                   = nodes->NumberOfDofsLocal(configuration_type,set1enum);
+	int n                   = nodes->NumberOfDofsLocal(configuration_type,set2enum);
 	int numnodesperobject   = elements->MaxNumNodes();
 
@@ -471,5 +469,5 @@
 	for(i=0;i<nodes->Size();i++){
 		Node* node=dynamic_cast<Node*>(nodes->GetObjectByOffset(i));
-		if(node->InAnalysis(analysis_type)){
+		if(node->InAnalysis(configuration_type)){
 
 			/*Reinitialize flags to 0*/
@@ -509,5 +507,5 @@
 		for(i=0;i<nodes->Size();i++){
 			Node* node=dynamic_cast<Node*>(nodes->GetObjectByOffset(i));
-			if(node->InAnalysis(analysis_type) && !node->IsClone()){
+			if(node->InAnalysis(configuration_type) && !node->IsClone()){
 				for(j=0;j<node->indexing.fsize;j++){
 					_assert_(count<m);
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 13908)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Penta.cpp	(revision 13909)
@@ -904,4 +904,6 @@
 /*FUNCTION Penta::GetNumberOfNodes{{{*/
 int Penta::GetNumberOfNodes(void){
+
+	if(this->nodes==NULL) return 0;
 
 	switch(this->element_type){
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.cpp	(revision 13908)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.cpp	(revision 13909)
@@ -70,5 +70,4 @@
 void PentaHook::SetHookNodes(int* node_ids,int analysis_counter){
 	this->hnodes[analysis_counter]= new Hook(node_ids,6);
-
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.h	(revision 13908)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/PentaHook.h	(revision 13909)
@@ -13,18 +13,18 @@
 
 	public: 
-		int   numanalyses; //number of analysis types
-		Hook** hnodes; // 6 nodes for each analysis type
-		Hook*  hmaterial; // 1 ice material
-		Hook*  hmatpar; // 1 material parameter
-		Hook*  hneighbors; // 2 elements, first down, second up
+		int    numanalyses;   //number of analysis types
+		Hook **hnodes;        // 6 nodes for each analysis type
+		Hook  *hmaterial;     // 1 ice material
+		Hook  *hmatpar;       // 1 material parameter
+		Hook  *hneighbors;    // 2 elements, first down, second up
 
-		/*FUNCTION constructors, destructors {{{*/
+		/*constructors, destructors*/
 		PentaHook();
 		PentaHook(int in_numanalyses,int material_id, IoModel* iomodel);
 		~PentaHook();
+
 		void SetHookNodes(int* node_ids,int analysis_counter);
 		void SpawnTriaHook(TriaHook* triahook,int* indices);
 		void InitHookNeighbors(int* element_ids);
-		/*}}}*/
 };
 
Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 13908)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 13909)
@@ -1139,4 +1139,6 @@
 /*FUNCTION Tria::GetNumberOfNodes{{{*/
 int Tria::GetNumberOfNodes(void){
+
+	if(this->nodes==NULL) return 0;
 
 	switch(this->element_type){