Index: /issm/trunk/src/m/solutions/jpl/NewFemModel.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/NewFemModel.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/NewFemModel.m	(revision 4119)
@@ -1,3 +1,3 @@
-function femmodel=NewFemModel(md,solution_type,analysis_types);
+function femmodel=NewFemModel(md,solution_type,analysis_types,nummodels);
 %NEWFEMMODEL - create a finite element model out of the matlab base \@model md. 
 %   For each analysis_type contained in analysis_types, create a set of nodes, constraints 
@@ -6,5 +6,5 @@
 %
 %   Usage:
-%      femmodel=NewFemModel(md,solution_type,analysis_types)
+%      femmodel=NewFemModel(md,solution_type,analysis_types,nummodels)
 %
 
@@ -23,5 +23,5 @@
 	   
 	   analysis_type=femmodel.analysis_type_list(i);
-	   displaystring(md.verbose,'%s%s','   dealing with analysis type: ',EnumAsString(analysis_type);
+	   displaystring(md.verbose,'%s%s','   dealing with analysis type: ',EnumAsString(analysis_type));
 
 	   displaystring(md.verbose,'%s','      generating degrees of freedofemmodel...');
Index: /issm/trunk/src/m/solutions/jpl/balancedthickness.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedthickness.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/balancedthickness.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/balancedthickness2.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedthickness2.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/balancedthickness2.m	(revision 4119)
@@ -9,5 +9,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/balancedvelocities.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/balancedvelocities.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/balancedvelocities.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/bedslope.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/bedslope.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/bedslope.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/diagnostic.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/diagnostic.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,5);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/prognostic.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/prognostic.m	(revision 4119)
@@ -13,5 +13,5 @@
 	
 	displaystring(md.verbose,'%s',['reading prognostic model data']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	% figure out number of dof: just for information purposes.
Index: /issm/trunk/src/m/solutions/jpl/prognostic2.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/prognostic2.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/prognostic2.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/steadystate.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/steadystate.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/steadystate.m	(revision 4119)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,7);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/surfaceslope.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/surfaceslope.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/surfaceslope.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,1);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/jpl/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/thermal.m	(revision 4118)
+++ /issm/trunk/src/m/solutions/jpl/thermal.m	(revision 4119)
@@ -12,5 +12,5 @@
 
 	displaystring(md.verbose,'%s',['reading thermal model data']);
-	femmodel=NewFemModel(md,solution_type,analysis_types);
+	femmodel=NewFemModel(md,solution_type,analysis_types,2);
 
 	%retrieve parameters