Index: /issm/trunk/src/m/solutions/NewFemModel.m
===================================================================
--- /issm/trunk/src/m/solutions/NewFemModel.m	(revision 4196)
+++ /issm/trunk/src/m/solutions/NewFemModel.m	(revision 4197)
@@ -8,5 +8,4 @@
 %      femmodel=NewFemModel(md,solution_type,analysis_types,nummodels)
 %
-
 
    femmodel.solution_type=solution_type;
@@ -27,8 +26,10 @@
 
    %now, go through all analyses types and post-process datasets
-   for i=1:length(analysis_types),
+   for i=1:nummodels,
 	   
 	   analysis_type=femmodel.analysis_type_list(i);
 	   displaystring(md.verbose,'%s%s','   dealing with analysis type: ',EnumAsString(analysis_type));
+
+		femmodel=SetCurrentAnalysis(femmodel,analysis_type);
 
 	   displaystring(md.verbose,'%s','      generating degrees of freedofemmodel...');
@@ -56,7 +57,6 @@
 		displaystring(md.verbose,'%s','      configuring element and loads...');
 		[femmodel.elements,femmodel.loads,femmodel.nodes,femmodel.parameters] = ConfigureObjects( femmodel.elements, femmodel.loads, femmodel.nodes, femmodel.vertices,femmodel.materials,femmodel.parameters);
+		disp('file: NewFemModel.m line: 59');
+
 
    end
-
-	displaystring(md.verbose,'%s','   process parameters...');
-	femmodel.parameters=ProcessParams(femmodel.parameters);
Index: /issm/trunk/src/m/solutions/SetCurrentAnalysis.m
===================================================================
--- /issm/trunk/src/m/solutions/SetCurrentAnalysis.m	(revision 4196)
+++ /issm/trunk/src/m/solutions/SetCurrentAnalysis.m	(revision 4197)
@@ -11,6 +11,8 @@
 
 	%first, look for analysis: 
-	found=0;
+	found=-1;
 	for i=1:length(femmodel.analysis_type_list),
+		femmodel.analysis_type_list(i)
+		analysis_enum
 		if femmodel.analysis_type_list(i)==analysis_enum,
 			found=i;
@@ -19,8 +21,7 @@
 	end
 
-	if ~found==1,
+	if found==-1,
 		error('SetCurrentAnalysis error message: could not find analysis_type in list of FemModel analyses');
 	end
-
 	
 	%activate matrices and vectors: 
@@ -34,4 +35,4 @@
 	%Now, plug analysis_counter and analysis_type inside the parameters: 
 	%set counter and analyse_type
-	femmodel.parameters.AnalysisCounter=found;
+	femmodel.parameters.AnalysisCounter=found-1; %for c indices
 	femmodel.parameters.AnalysisType=analysis_enum;
