Index: /issm/trunk/src/m/solutions/cielo/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 854)
+++ /issm/trunk/src/m/solutions/cielo/diagnostic_core.m	(revision 855)
@@ -21,6 +21,6 @@
 	if dim==3,
 		displaystring(debug,'\n%s',['extruding slopes in 3d...']);
-		slopex=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopex,"slopex",0);
-		slopey=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopey,"slopey",0);
+		slopex=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopex,'slopex',0);
+		slopey=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopey,'slopey',0);
 	end
 
@@ -55,5 +55,5 @@
 
 	displaystring(debug,'\n%s',['extruding horizontal velocities...']);
-	u_g_horiz=FieldExtrude(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,u_g,"velocity",1);
+	u_g_horiz=FieldExtrude(m_dh.elements,m_dh.nodes,m_dh.loads,m_dh.materials,u_g,'velocity',1);
 
 	displaystring(debug,'\n%s',['computing vertical velocities...']);
@@ -77,6 +77,6 @@
 		slopex=diagnostic_core_linear(m_sl,inputs,'slope_compute','bedx');
 		slopey=diagnostic_core_linear(m_sl,inputs,'slope_compute','bedy');
-		slopex=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopex,"slopex",0);
-		slopey=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopey,"slopey",0);
+		slopex=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopex,'slopex',0);
+		slopey=FieldExtrude(m_sl.elements,m_sl.nodes,m_sl.loads,m_sl.materials,slopey,'slopey',0);
 
 		inputs=add(inputs,'bedslopex',slopex,'doublevec',m_sl.parameters.numberofdofspernode,m_sl.parameters.numberofnodes);