Index: /issm/trunk/src/c/ModelProcessorx/CreateParameters.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/CreateParameters.cpp	(revision 3773)
+++ /issm/trunk/src/c/ModelProcessorx/CreateParameters.cpp	(revision 3774)
@@ -30,4 +30,6 @@
 	if (strcmp(iomodel->meshtype,"2d")==0) parameters->AddObject(new IntParam(DimEnum,2));
 	else parameters->AddObject(new IntParam(DimEnum,3));
+	parameters->AddObject(new    IntParam(AnalysisTypeEnum,iomodel->analysis_type));
+	parameters->AddObject(new    IntParam(SubAnalysisTypeEnum,iomodel->sub_analysis_type));
 	parameters->AddObject(new StringParam(OutputFileNameEnum,iomodel->outputfilename));
 	parameters->AddObject(new   BoolParam(IsHutterEnum,iomodel->ishutter));
Index: /issm/trunk/src/c/parallel/prognostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/prognostic_core.cpp	(revision 3773)
+++ /issm/trunk/src/c/parallel/prognostic_core.cpp	(revision 3774)
@@ -34,5 +34,5 @@
 
 	/*recover fem model: */
-	fem_p=model->GetFormulation(PrognosticAnalysisEnum);
+	fem_p=model->GetFormulation(PrognosticAnalysisEnum); ISSMASSERT(fem_p);
 
 	//first recover parameters common to all solutions