Index: /issm/trunk/src/m/solutions/cielo/ControlInitialization.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/ControlInitialization.m	(revision 1702)
+++ /issm/trunk/src/m/solutions/cielo/ControlInitialization.m	(revision 1703)
@@ -38,5 +38,5 @@
 %horizontal velocities
 displaystring(debug,'\n%s',['computing horizontal velocities...']);
-u_g=diagnostic_core_nonlinear(m_dh,inputs,'diagnostic','horiz');
+u_g=diagnostic_core_nonlinear(m_dh,inputs,DiagnosticAnalysisEnum(),HorizAnalysisEnum());
 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);
@@ -45,5 +45,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,'diagnostic','vert');
+u_g_vert=diagnostic_core_linear(m_dv,inputs,DiagnosticAnalysisEnum(),VertAnalysisEnum());
 
 %create 3d u_g
@@ -70,5 +70,5 @@
 %Compute Stokes velocities once to have a reasonably good ug in input
 displaystring(debug,'\n%s',['computing stokes velocities and pressure ...']);
-u_g=diagnostic_core_nonlinear(m_ds,inputs,'diagnostic','stokes');
+u_g=diagnostic_core_nonlinear(m_ds,inputs,DiagnosticAnalysisEnum(),StokesAnalysisEnum());
 inputs=add(inputs,'velocity',u_g,'doublevec',4,m_ds.parameters.numberofnodes);