Index: /issm/trunk/src/m/solutions/ice/thermal.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/thermal.m	(revision 951)
+++ /issm/trunk/src/m/solutions/ice/thermal.m	(revision 952)
@@ -60,5 +60,5 @@
 	melting(1:6:gridset.gsize)=md.melting/md.yts; %in m/s
 	soln.t_g=temperature;
-	soln.melting_g=melting;
+	soln.m_g=melting;
 	soln.time=0;
 
@@ -77,6 +77,7 @@
 		disp('   computing melting...');
 		gridset=m_m.gridset;
-		inputs=struct('pressure',pressure,'temperature',soln(n).t_g,'melting_offset',melting_offset,'dt',md.dt);
-		soln(n+1).melting_g=melting_core(m_m,analysis_type,inputs);
+		inputs=struct('pressure',pressure,'temperature',soln(n+1).t_g,'melting_offset',melting_offset,'dt',md.dt);
+		soln(n+1).m_g=melting_core(m_m,analysis_type,inputs);
+		any(soln(n+1).m_g)
 		
 	end
@@ -90,5 +91,5 @@
 		solution(n).step=n;
 		solution(n).temperature=soln(n).t_g(indx);
-		solution(n).melting=soln(n).melting_g(indx)*md.yts; %in m/a
+		solution(n).melting=soln(n).m_g(indx)*md.yts; %in m/a
 		solution(n).time=soln(n).time/md.yts; %in year
 	end