Index: /issm/trunk/src/m/solutions/jpl/ModelUpdateInputsFromVector.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/ModelUpdateInputsFromVector.m	(revision 3852)
+++ /issm/trunk/src/m/solutions/jpl/ModelUpdateInputsFromVector.m	(revision 3853)
@@ -24,5 +24,7 @@
 	field=modelfields(i); field=field{1}; model=models.(field);
 
-	[model.elements,model.nodes,model.vertices,model.loads,model.materials,model.parameters] = UpdateInputsFromVector(model.elements,model.nodes,model.vertices,model.loads,model.materials,model.parameters,vector,enum, typeenum);
+	if isstruct(model), %there is an analysis_type model
+		[model.elements,model.nodes,model.vertices,model.loads,model.materials,model.parameters] = UpdateInputsFromVector(model.elements,model.nodes,model.vertices,model.loads,model.materials,model.parameters,vector,enum, typeenum);
+	end
 
 end
Index: /issm/trunk/src/m/solutions/jpl/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/diagnostic_core.m	(revision 3852)
+++ /issm/trunk/src/m/solutions/jpl/diagnostic_core.m	(revision 3853)
@@ -72,5 +72,5 @@
 	u_g_horiz=FieldExtrude(m_dh.elements,m_dh.nodes,m_dh.vertices,m_dh.loads,m_dh.materials,m_dh.parameters,u_g,'velocity',1);
 
-	[vx,vy]=SplitSolutionVector(ug,numberofnodes,numberofdofspernode_dh);
+	[vx,vy]=SplitSolutionVector(u_g_horiz,m_dh.parameters.NumberOfNodes,m_dh.parameters.NumberOfDofsPerNode);
 	ModelUpdateInputsFromVector(model,vx,VxEnum,VertexEnum);
 	ModelUpdateInputsFromVector(model,vy,VyEnum,VertexEnum);