Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 25489)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 25490)
@@ -5068,4 +5068,7 @@
 	int         counter  = 0; //this is probably not acurate but we are only interested in positivity
 	Element*    basalelement=NULL;
+	int         domaintype;
+	this->parameters->FindParam(&domaintype,DomainTypeEnum);
+
 	for(int i=0;i<elements->Size();i++){
 		Element    *element  = xDynamicCast<Element*>(elements->GetObjectByOffset(i));
@@ -5271,5 +5274,4 @@
 		else{
 			for(int j=0;j<elements->Size();j++){
-
 				/*Get the right transient input*/
 				Element* element=xDynamicCast<Element*>(elements->GetObjectByOffset(j));