Index: /issm/trunk/src/m/solutions/cielo/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 548)
+++ /issm/trunk/src/m/solutions/cielo/thermal.m	(revision 549)
@@ -44,5 +44,4 @@
 		t_g=m_t.parameters.t_g;
 		melting_g=m_t.parameters.melting_g;
-
 		nsteps=m_t.parameters.ndt/m_t.parameters.dt;
 
@@ -50,8 +49,10 @@
 		soln.t_g=t_g;
 		soln.melting_g=melting_g;
+		soln.time=0;
 
 		for n=1:nsteps, 
 
 			displaystring(md.debug,'\n%s%i/%i\n','time step: ',n,nsteps);
+			soln(n+1).time=n*m_t.parameters.dt;
 			
 			displaystring(md.debug,'\n%s',['    computing temperatures...']);
@@ -67,15 +68,11 @@
 		
 		%Wrap up
-		solution_temperature=struct('temperature',{});
+		solution=struct('time',{},'temperature',{},'melting',{});
 		for n=1:nsteps+1,
-			solution_temperature(n).temperature=soln(n).t_g;
+			solution(n).temperature=soln(n).t_g;
+			solution(n).melting=soln(n).melting_g*md.yts; %in m/year
+			solution(n).time=soln(n).time/md.yts;         %in years
 		end
-		md.temperature=solution_temperature;
-		
-		solution_melting=struct('melting',{});
-		for n=1:nsteps+1,
-			solution_melting(n).melting=soln(n).melting_g*md.yts; %in m/a
-		end
-		md.melting=solution_melting;
+		md.thermaltransient_results=solution;
 	end