Index: /issm/trunk/src/m/solutions/thermal_core_step.m
===================================================================
--- /issm/trunk/src/m/solutions/thermal_core_step.m	(revision 4468)
+++ /issm/trunk/src/m/solutions/thermal_core_step.m	(revision 4469)
@@ -1,7 +1,7 @@
-function femmodel=thermal_core_step(femmodel)
+function femmodel=thermal_core_step(femmodel,step,time)
 %THERMAL_CORE_STEP - core of the thermal solution for one step 
 %
 %   Usage:
-%      femmodel=thermal_core_step(femmodel)
+%      femmodel=thermal_core_step(femmodel,step,time)
 
 	%recover parameters common to all solutions
@@ -10,8 +10,8 @@
 	displaystring(verbose,'\n%s',['computing temperature...']);
 	femmodel=SetCurrentConfiguration(femmodel,ThermalAnalysisEnum);
-	femmodel=solver_thermal_nonlineat(femmodel);
+	femmodel=solver_thermal_nonlinear(femmodel);
 
 	displaystring(verbose,'\n%s',['compute melting...']);
-	femmodel=SetCurrentConfiguraton(femmodel,MeltingAnalysisEnum);
+	femmodel=SetCurrentConfiguration(femmodel,MeltingAnalysisEnum);
 	femmodel=solver_linear(femmodel);
 	
Index: /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m
===================================================================
--- /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m	(revision 4468)
+++ /issm/trunk/src/m/solvers/solver_thermal_nonlinear.m	(revision 4469)
@@ -22,8 +22,8 @@
 		if count==1
 			displaystring(femmodel.parameters.Verbose,'%s',['   system matrices']);
-			[K_gg_nopenalty, p_g_nopenalty]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,analysis_type,sub_analysis_type);
+			[K_gg_nopenalty, p_g_nopenalty]=SystemMatrices(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 		end
 		displaystring(femmodel.parameters.Verbose,'%s',['   penalty system matrices']);
-		[K_gg , p_g, melting_offset]=PenaltySystemMatrices(K_gg_nopenalty,p_g_nopenalty,femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,analysis_type,sub_analysis_type);
+		[K_gg , p_g, melting_offset]=PenaltySystemMatrices(K_gg_nopenalty,p_g_nopenalty,femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
 
 		%Reduce tangent matrix from g size to f size
@@ -44,5 +44,5 @@
 		%penalty constraints
 		displaystring(femmodel.parameters.Verbose,'%s',['   penalty constraints']);
-		[femmodel.loads,constraints_converged,num_unstable_constraints] =PenaltyConstraints(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads, femmodel.materials,femmodel.parameters,analysis_type,sub_analysis_type);
+		[femmodel.loads,constraints_converged,num_unstable_constraints] =PenaltyConstraints(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads, femmodel.materials,femmodel.parameters);
 	
 		if ~converged,
@@ -57,3 +57,6 @@
 	end
 
+	%add melting_offset  into inputs
+	[femmodel.elements femmodel.loads]=InputUpdateFromConstant(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,melting_offset,MeltingOffsetEnum);
 end
+