Index: /issm/trunk/src/m/solutions/cielo/thermal_core.m
===================================================================
--- /issm/trunk/src/m/solutions/cielo/thermal_core.m	(revision 528)
+++ /issm/trunk/src/m/solutions/cielo/thermal_core.m	(revision 529)
@@ -52,5 +52,5 @@
 		%penalty constraints
 		displaystring(m.parameters.debug,'%s',['   penalty constraints']);
-		[m.loads,constraints_converged,num_unstable_constraints] =PenaltyConstraints( m.elements,m.nodes, m.loads, m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
+		[loads,constraints_converged,num_unstable_constraints] =PenaltyConstraints( m.elements,m.nodes, loads, m.materials,m.parameters,inputs,analysis_type,sub_analysis_type);
 	
 		if ~converged,