Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 468)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 469)
@@ -59,4 +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);
+	which diagnostic_core_linear
 	u_g_vert=diagnostic_core_linear(m_dv,inputs,'diagnostic','vert');
 
Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m	(revision 468)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core_linear.m	(revision 469)
@@ -1,7 +1,7 @@
-function u_g=diagnostic_core_linear(m,inputs,analysis_type)
+function u_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type)
 %DIAGNOSTIC_CORE_LINEAR - linear solution sequence
 %
 %   Usage:
-%      u_g=diagnostic_core_linear(m,inputs,analysis_type)
+%      u_g=diagnostic_core_linear(m,inputs,analysis_type,sub_analysis_type)
 
 	%stiffness and load generation only:
@@ -12,5 +12,5 @@
 	
 	%system matrices
-	[K_gg, p_g]=SystemMatrices(m.elements,m.nodes,m.loads,m.materials,m.parameters,inputs,analysis_type);
+	[K_gg, p_g]=SystemMatrices(m.elements,m.nodes,m.loads,m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
 	
 	%Reduce tangent matrix from g size to f size