Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 325)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 326)
@@ -16,6 +16,6 @@
 
 	displaystring(debug,'\n%s',['computing surface slope (x and y derivatives)...']);
-	slopex=diagnostic_core_linear(m_ss,'surface_slope_compute_x',inputs);
-	slopey=diagnostic_core_linear(m_ss,'surface_slope_compute_y',inputs);
+	slopex=diagnostic_core_linear(m_ss,inputs,'surface_slope_compute_x');
+	slopey=diagnostic_core_linear(m_ss,inputs,'surface_slope_compute_y');
 
 	if dim==3,
@@ -30,5 +30,5 @@
 
 	displaystring(debug,'\n%s',['computing hutter velocities...']);
-	u_g=diagnostic_core_linear(m_dhu,'diagnostic_hutter',inputs);
+	u_g=diagnostic_core_linear(m_dhu,inputs,'diagnostic_hutter');
 
 	displaystring(debug,'\n%s',['computing pressure according to MacAyeal...']);
@@ -59,5 +59,5 @@
 	displaystring(debug,'\n%s',['computing vertical velocities...']);
 	inputs=add(inputs,'velocity',u_g_horiz,'doublevec',m_dh.parameters.numberofdofspernode,m_dh.parameters.numberofnodes);
-	u_g_vert=diagnostic_core_linear(m_dv,inputs);
+	u_g_vert=diagnostic_core_linear(m_dv,inputs,'diagnostic_vert');
 
 	displaystring(debug,'\n%s',['combining horizontal and vertical velocities...']);
@@ -75,6 +75,6 @@
 
 		displaystring(debug,'\n%s',['computing bed slope (x and y derivatives)...']);
-		slopex=diagnostic_core_linear(m_bs,'bed_slope_compute_x',inputs);
-		slopey=diagnostic_core_linear(m_bs,'bed_slope_compute_y',inputs);
+		slopex=diagnostic_core_linear(m_bs,inputs,'bed_slope_compute_x');
+		slopey=diagnostic_core_linear(m_bs,inputs,'bed_slope_compute_y');
 		slopex=SlopeExtrude(m_bs.elements,m_bs.nodes,m_bs.loads,m_bs.materials,slopex);
 		slopey=SlopeExtrude(m_bs.elements,m_bs.nodes,m_bs.loads,m_bs.materials,slopey);
Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m	(revision 325)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m	(revision 326)
@@ -1,7 +1,7 @@
-function u_g=diagnostic_core_linear(m,inputs)
+function u_g=diagnostic_core_linear(m,inputs,analysis_type)
 %DIAGNOSTIC_CORE_LINEAR - linear solution sequence
 %
 %   Usage:
-%      u_g=diagnostic_core_linear(m,inputs)
+%      u_g=diagnostic_core_linear(m,inputs,analysis_type)
 
 	%stiffness and load generation only: