Index: /issm/trunk/src/m/solutions/balancedthickness.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/balancedthickness.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/balancedthickness2.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedthickness2.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/balancedthickness2.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/balancedvelocities.m
===================================================================
--- /issm/trunk/src/m/solutions/balancedvelocities.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/balancedvelocities.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/bedslope.m
===================================================================
--- /issm/trunk/src/m/solutions/bedslope.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/bedslope.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/prognostic.m	(revision 4296)
@@ -14,5 +14,5 @@
 	
 	displaystring(md.verbose,'%s',['reading prognostic model data']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	% figure out number of dof: just for information purposes.
Index: /issm/trunk/src/m/solutions/prognostic2.m
===================================================================
--- /issm/trunk/src/m/solutions/prognostic2.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/prognostic2.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/steadystate.m
===================================================================
--- /issm/trunk/src/m/solutions/steadystate.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/steadystate.m	(revision 4296)
@@ -9,10 +9,10 @@
 	t1=clock;
 
-	numanalyses=1;
+	numanalyses=7;
 	analyses=[DiagnosticHorizAnalysisEnum;DiagnosticVertAnalysisEnum;DiagnosticStokesAnalysisEnum;DiagnosticHutterAnalysisEnum;SlopeAnalysisEnum;ThermalAnalysisEnum;MeltingAnalysisEnum];
 	solution_type=SteadyStateSolutionEnum;
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,7);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/surfaceslope.m
===================================================================
--- /issm/trunk/src/m/solutions/surfaceslope.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/surfaceslope.m	(revision 4296)
@@ -13,5 +13,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,1);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/thermal.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/thermal.m	(revision 4296)
@@ -8,10 +8,10 @@
 	t1=clock;
 	
-	numanalyses=1;
+	numanalyses=2;
 	analyses=[ThermalAnalysisEnum;MeltingAnalysisEnum];
 	solution_type=ThermalSolutionEnum;
 
 	displaystring(md.verbose,'%s',['reading thermal model data']);
-	femmodel=NewFemModel(md,solution_type,analyses,2);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/transient2d.m
===================================================================
--- /issm/trunk/src/m/solutions/transient2d.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/transient2d.m	(revision 4296)
@@ -10,10 +10,10 @@
 	t1=clock;
 
-	numanalyses=1;
+	numanalyses=2;
 	analyses=[DiagnosticHorizAnalysisEnum;PrognosticAnalysisEnum];
 	solution_type=Transient2DSolutionEnum;
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,2);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters
Index: /issm/trunk/src/m/solutions/transient3d.m
===================================================================
--- /issm/trunk/src/m/solutions/transient3d.m	(revision 4295)
+++ /issm/trunk/src/m/solutions/transient3d.m	(revision 4296)
@@ -14,5 +14,5 @@
 
 	displaystring(md.verbose,'%s',['create finite element model']);
-	femmodel=NewFemModel(md,solution_type,analyses,8);
+	femmodel=NewFemModel(md,solution_type,analyses,numanalyses);
 
 	%retrieve parameters