Index: /issm/trunk/src/m/solutions/cielo/prognostic.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 1719)
+++ /issm/trunk/src/m/solutions/cielo/prognostic.m	(revision 1720)
@@ -24,5 +24,5 @@
 	inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofnodes);
 	inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofnodes);
-	inputs=add(inputs,'dt',models.p.parameters.dt,'double');
+	inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
 
 	displaystring(md.debug,'\n%s',['call computational core:']);
Index: /issm/trunk/src/m/solutions/cielo/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 1719)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 1720)
@@ -25,5 +25,5 @@
 	inputs=add(inputs,'velocity',models.t.parameters.u_g,'doublevec',3,models.t.parameters.numberofnodes);
 	inputs=add(inputs,'pressure',models.t.parameters.p_g,'doublevec',1,models.t.parameters.numberofnodes);
-	inputs=add(inputs,'dt',models.t.parameters.dt,'double');
+	inputs=add(inputs,'dt',models.t.parameters.dt*models.t.parameters.yts,'double');
 	
 	%call core
Index: /issm/trunk/src/m/solutions/cielo/transient2d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 1719)
+++ /issm/trunk/src/m/solutions/cielo/transient2d.m	(revision 1720)
@@ -44,5 +44,5 @@
 inputs=add(inputs,'melting',models.p.parameters.m_g,'doublevec',1,models.p.parameters.numberofnodes);
 inputs=add(inputs,'accumulation',models.p.parameters.a_g,'doublevec',1,models.p.parameters.numberofnodes);
-inputs=add(inputs,'dt',models.p.parameters.dt,'double');
+inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double'); %change time from yrs to sec
 
 % figure out number of dof: just for information purposes.
@@ -58,5 +58,5 @@
 while  time<finaltime+dt, %make sure we run up to finaltime.
 
-	disp(sprintf('\n%s%g%s%g%s%g\n','time [yr]: ',time/yts,'    iteration number: ',n,'/',floor(finaltime/dt)));
+	disp(sprintf('\n%s%g%s%g%s%g\n','time [yr]: ',time,'    iteration number: ',n,'/',floor(finaltime/dt)));
 
 	solution(n+1).step=n+1;
Index: /issm/trunk/src/m/solutions/cielo/transient3d.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 1719)
+++ /issm/trunk/src/m/solutions/cielo/transient3d.m	(revision 1720)
@@ -56,5 +56,5 @@
 inputs=add(inputs,'melting',results.m_g,'doublevec',1,models.p.parameters.numberofnodes);
 inputs=add(inputs,'accumulation',results.a_g,'doublevec',1,models.p.parameters.numberofnodes);
-inputs=add(inputs,'dt',models.p.parameters.dt,'double');
+inputs=add(inputs,'dt',models.p.parameters.dt*models.p.parameters.yts,'double');
 
 % figure out number of dof: just for information purposes.
@@ -70,5 +70,5 @@
 while  time<finaltime+dt, %make sure we run up to finaltime.
 
-	displaystring(md.debug,'\n%s%g%s%g%s%g\n','time [yr]: ',time/yts,'    iteration number: ',n,'/',floor(finaltime/dt));
+	displaystring(md.debug,'\n%s%g%s%g%s%g\n','time [yr]: ',time,'    iteration number: ',n,'/',floor(finaltime/dt));
 
 	results(n+1).step=n+1;