Index: /issm/trunk/src/c/parallel/thermal.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal.cpp	(revision 578)
+++ /issm/trunk/src/c/parallel/thermal.cpp	(revision 579)
@@ -120,4 +120,5 @@
 		nsteps=(int)(ndt/dt);
 		time=(double*)xmalloc((nsteps+1)*sizeof(double));
+		time[0]=0;
 
 		/*allocate t_g and m_g arrays: */
@@ -136,6 +137,6 @@
 
 		for(i=0;i<nsteps;i++){
-			if(debug)_printf_("time step: %i/%i\n",n,nsteps);
-			time[i]=n*dt;
+			if(debug)_printf_("time step: %i/%i\n",i+1,nsteps);
+			time[i+1]=(i+1)*dt;
 			
 			if(debug)_printf_("computing temperatures:\n");
@@ -146,5 +147,5 @@
 			inputs->Add("temperature",t_g[i+1],1,numberofnodes);
 			inputs->Add("melting_offset",melting_offset);
-			diagnostic_core_linear(&m_g[0],&femmodels[1],inputs,MeltingAnalysisEnum(),TransientAnalysisEnum());
+			diagnostic_core_linear(&m_g[i+1],&femmodels[1],inputs,MeltingAnalysisEnum(),TransientAnalysisEnum());
 		}
 	}