Index: /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 23170)
+++ /issm/trunk-jpl/src/c/analyses/LevelsetAnalysis.cpp	(revision 23171)
@@ -12,5 +12,7 @@
 
 void LevelsetAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/
-	IoModelToConstraintsx(constraints,iomodel,"md.levelset.spclevelset",LevelsetAnalysisEnum,P1Enum);
+	int finiteelement;
+	iomodel->FindConstant(&finiteelement,"md.levelset.fe");
+	IoModelToConstraintsx(constraints,iomodel,"md.levelset.spclevelset",LevelsetAnalysisEnum,finiteelement);
 }
 /*}}}*/
@@ -19,5 +21,6 @@
 }/*}}}*/
 void LevelsetAnalysis::CreateNodes(Nodes* nodes,IoModel* iomodel){/*{{{*/
-	int finiteelement=P1Enum;
+	int finiteelement;
+	iomodel->FindConstant(&finiteelement,"md.levelset.fe");
 	if(iomodel->domaintype!=Domain2DhorizontalEnum) iomodel->FetchData(2,"md.mesh.vertexonbase","md.mesh.vertexonsurface");
 	::CreateNodes(nodes,iomodel,LevelsetAnalysisEnum,finiteelement);
@@ -32,5 +35,6 @@
 
 	/*Finite element type*/
-	int finiteelement = P1Enum;
+	int finiteelement;
+	iomodel->FindConstant(&finiteelement,"md.levelset.fe");
 
 	/*Update elements: */