Index: /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 17463)
@@ -853,5 +853,4 @@
 	femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
 	if(islevelset){
-		_printf0_("   Updating active and non-active nodes for ThermalAnalysis \n");
 		SetActiveNodesLSMx(femmodel);
 	}
Index: /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/ExtrapolationAnalysis.cpp	(revision 17463)
@@ -12,4 +12,5 @@
 void ExtrapolationAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
 	//do nothing for now
+	return;
 }
 /*}}}*/
@@ -44,8 +45,10 @@
 void ExtrapolationAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/
 	// do nothing for now
+	return;
 }
 /*}}}*/
 void ExtrapolationAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/
-// 	do nothing for now
+	//	do nothing for now
+	return;
 }/*}}}*/
 
@@ -70,5 +73,4 @@
 		femmodel->RequestedOutputsx(&femmodel->results,&outputs[0],2);
 	}
-
 }/*}}}*/
 ElementVector* ExtrapolationAnalysis::CreateDVector(Element* element){/*{{{*/
@@ -348,4 +350,3 @@
 		this->SetConstraintsOnIce(element);
 	}
-
-}/*}}}*/
+}/*}}}*/
Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 17463)
@@ -11,5 +11,5 @@
 /*}}}*/
 void LevelsetAnalysis::UpdateParameters(Parameters* parameters,IoModel* iomodel,int solution_enum,int analysis_enum){/*{{{*/
-	// parameters->AddObject(iomodel->CopyConstantObject(FlowequationIsSIAEnum));
+	return;
 }
 /*}}}*/
@@ -34,9 +34,4 @@
 	iomodel->FetchDataToInput(elements,VyEnum);
 	iomodel->FetchDataToInput(elements,MasstransportCalvingrateEnum);
-	
-	if(iomodel->meshtype==Mesh3DEnum){
-		iomodel->FetchDataToInput(elements,MeshElementonbedEnum);
-		iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
-	}
 }
 /*}}}*/
@@ -49,17 +44,9 @@
 /*}}}*/
 void LevelsetAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/
-
-	/*Intermediary*/
-
-	/*Fetch parameters: */
-
+	return;
 }
 /*}}}*/
 void LevelsetAnalysis::CreateLoads(Loads* loads, IoModel* iomodel){/*{{{*/
-	
-	/*Intermediary*/
-
-	/*Fetch parameters: */
-
+	return;
 }/*}}}*/
 
@@ -73,5 +60,4 @@
 	/*parameters: */
 	bool save_results;
-	
 	femmodel->parameters->FindParam(&save_results,SaveResultsEnum);
 
@@ -525,3 +511,2 @@
 	return fabs(a[0]*b[1]-a[1]*b[0])/norm_b;
 }/*}}}*/
-
Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp	(revision 17463)
@@ -755,5 +755,4 @@
 	femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
 	if(islevelset){
-		_printf0_("   Updating active and non-active nodes for MasstransportAnalysis \n");
 		SetActiveNodesLSMx(femmodel);
 	}
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17463)
@@ -1038,5 +1038,4 @@
 	femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
 	if(islevelset){
-		_printf0_("   Updating active and non-active nodes for StressbalanceAnalysis \n");
 		SetActiveNodesLSMx(femmodel);
 	}
Index: /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp	(revision 17463)
@@ -580,5 +580,4 @@
 	femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
 	if(islevelset){
-		_printf0_("   Updating active and non-active nodes for StressbalanceSIAAnalysis \n");
 		SetActiveNodesLSMx(femmodel);
 	}
Index: /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/analyses/ThermalAnalysis.cpp	(revision 17463)
@@ -717,5 +717,4 @@
 	femmodel->parameters->FindParam(&islevelset,TransientIslevelsetEnum);
 	if(islevelset){
-		_printf0_("   Updating active and non-active nodes for ThermalAnalysis \n");
 		SetActiveNodesLSMx(femmodel);
 	}
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17463)
@@ -2912,6 +2912,8 @@
 	/*If the level set has always same sign, there is no ice front here*/
 	iszerols = false;
-	if(ls[0]*ls[1]<0. || ls[0]*ls[2]<0. || (ls[0]*ls[1]*ls[2]==0. && ls[0]*ls[1]+ls[0]*ls[2]+ls[1]*ls[2]<=0.)){
-		iszerols = true;
+	if(IsIceInElement()){
+		if(ls[0]*ls[1]<0. || ls[0]*ls[2]<0. || (ls[0]*ls[1]*ls[2]==0. && ls[0]*ls[1]+ls[0]*ls[2]+ls[1]*ls[2]<=0.)){
+			iszerols = true;
+		}
 	}
 	return iszerols;
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17463)
@@ -2347,6 +2347,8 @@
 	/*If the level set is awlays <0, there is no ice front here*/
 	iszerols= false;
-	if(ls[0]*ls[1]<0. || ls[0]*ls[2]<0. || (ls[0]*ls[1]*ls[2]==0. && ls[0]*ls[1]+ls[0]*ls[2]+ls[1]*ls[2]<=0.)){
-		iszerols = true;
+	if(IsIceInElement()){
+		if(ls[0]*ls[1]<0. || ls[0]*ls[2]<0. || (ls[0]*ls[1]*ls[2]==0. && ls[0]*ls[1]+ls[0]*ls[2]+ls[1]*ls[2]<=0.)){
+			iszerols = true;
+		}
 	}
 
Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17463)
@@ -407,5 +407,5 @@
 
 	/*start module: */
-	if(VerboseModule()) _printf0_("   Updating constraints for time: " << time << "\n");
+	if(VerboseModule()) _printf0_("   Updating constraints and active domain of analysis " << EnumToStringx(analysis_type)  << " for time: " << time << "\n");
 
 	Analysis* analysis= EnumToAnalysis(analysis_type);
Index: /issm/trunk-jpl/src/c/cores/transient_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/cores/transient_core.cpp	(revision 17463)
@@ -93,11 +93,4 @@
 	}
 
-// 	if(islevelset){
-// 		/* set distance on elements intersected by zero level set */
-// 		LevelsetAnalysis* lsfanalysis = new LevelsetAnalysis();
-// 		lsfanalysis->SetDistanceOnIntersectedElements(femmodel);
-// 		delete lsfanalysis;
-// 	}
-	
 	while(time < finaltime - (yts*DBL_EPSILON)){ //make sure we run up to finaltime.
 
Index: /issm/trunk-jpl/src/c/modules/SetActiveNodesLSMx/SetActiveNodesLSMx.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/SetActiveNodesLSMx/SetActiveNodesLSMx.cpp	(revision 17462)
+++ /issm/trunk-jpl/src/c/modules/SetActiveNodesLSMx/SetActiveNodesLSMx.cpp	(revision 17463)
@@ -49,4 +49,5 @@
 		}
 
+		/* if solving 2d problem on vertically extende mesh, solve on basal layer only*/
 		if(meshtype!=Mesh2DhorizontalEnum){
 			femmodel->parameters->FindParam(&analysis_type,AnalysisTypeEnum);
@@ -97,4 +98,5 @@
 	/* Intermediaries */
 	int i;
+
 	/*Initialize vector with number of vertices*/
 	int numvertices=femmodel->vertices->NumberOfVertices();