Index: /issm/trunk/src/c/parallel/diagnostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 1830)
+++ /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 1831)
@@ -13,5 +13,5 @@
 #include "../issm.h"
 
-void diagnostic_core(DataSet* results,FemModel* fems, ParameterInputs* inputs){
+void diagnostic_core(DataSet* results,Model* model, ParameterInputs* inputs){
 
 	extern int my_rank;
@@ -58,9 +58,9 @@
 
 	/*recover fem models: */
-	fem_dh=fems+0;
-	fem_dv=fems+1;
-	fem_ds=fems+2;
-	fem_dhu=fems+3;
-	fem_sl=fems+4;
+	fem_dh=model->DiagnosticHorizontal();
+	fem_dv=model->DiagnosticVertical();
+	fem_ds=model->DiagnosticStokes();
+	fem_dhu=model->DiagnosticHutter();
+	fem_sl=model->Slope();
 
 	//first recover parameters common to all solutions
@@ -186,5 +186,4 @@
 	result=new Result(results->Size()+1,0,1,"p_g",pg);
 	results->AddObject(result);
-
 	/*output if we have rifts: */
 	if(numrifts){