Index: /issm/trunk/src/m/solutions/cielo/diagnostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 103)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic.m	(revision 104)
@@ -14,6 +14,8 @@
 
 	if strcmpi(md.type,'3d'),
-		m_dbv=CreateFemModel(md,analysis_dbv);
-		m_dv=CreateFemModel(md,analysis_dv);
+		md.analysis_type=analysis_dbv; 
+		m_dbv=CreateFemModel(md);
+		md.analysis_type=analysis_dv; 
+		m_dv=CreateFemModel(md);
 	end
 
@@ -28,9 +30,6 @@
 	if strcmpi(md.type,'3d'),
 	
-		%swith u_g to workspace partitioning.
-		u_g=SwitchPartitioning(u_g,'workspace',part,tpart,[1 2]); 
-		
 		%extrude velocities for collapsed penta elements
-		u_g=CieloVelocityExtrude(md,u_g);
+		u_g=VelocityExtrude(md,u_g);
 
 		%Compute depth averaged velocity and add it to inputs
@@ -41,5 +40,4 @@
 		velocity_average=SwitchPartitioning(velocity_average,'cluster',part,tpart,[1 2]); 
 		
-
 		disp(sprintf('\n%s',['computing basal vertical velocities...']));
 		SetUset(m_dbv);