Index: /issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp	(revision 574)
+++ /issm/trunk/src/c/ModelProcessorx/CreateDataSets.cpp	(revision 575)
@@ -73,7 +73,19 @@
 	else if (strcmp(model->analysis_type,"thermal")==0){
 
+		PetscSynchronizedPrintf(MPI_COMM_WORLD,"ok1\n");
+		PetscSynchronizedFlush(MPI_COMM_WORLD);
+
 		CreateElementsNodesAndMaterialsThermal(pelements,pnodes,pmaterials, model,model_handle);
+		PetscSynchronizedPrintf(MPI_COMM_WORLD,"ok2\n");
+		PetscSynchronizedFlush(MPI_COMM_WORLD);
+
 		CreateConstraintsThermal(pconstraints,model,model_handle);
+		PetscSynchronizedPrintf(MPI_COMM_WORLD,"ok3\n");
+		PetscSynchronizedFlush(MPI_COMM_WORLD);
+
 		CreateLoadsThermal(ploads,model,model_handle);
+		PetscSynchronizedPrintf(MPI_COMM_WORLD,"ok4\n");
+		PetscSynchronizedFlush(MPI_COMM_WORLD);
+
 		CreateParametersThermal(pparameters,model,model_handle);
 					
Index: /issm/trunk/src/c/ModelProcessorx/Model.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Model.cpp	(revision 574)
+++ /issm/trunk/src/c/ModelProcessorx/Model.cpp	(revision 575)
@@ -58,4 +58,6 @@
 	model->vz=NULL;
 	model->pressure=NULL;
+	model->temperature=NULL;
+	model->melting=NULL;
 	model->elementonbed=NULL;
 	model->elementonsurface=NULL;
@@ -213,4 +215,6 @@
 	xfree((void**)&model->vz);
 	xfree((void**)&model->pressure);
+	xfree((void**)&model->temperature);
+	xfree((void**)&model->melting);
 	xfree((void**)&model->drag);
 	xfree((void**)&model->p);
Index: /issm/trunk/src/c/ModelProcessorx/Model.h
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Model.h	(revision 574)
+++ /issm/trunk/src/c/ModelProcessorx/Model.h	(revision 575)
@@ -47,4 +47,5 @@
 	double* vz;
 	double* pressure;
+	double* temperature;
 
 	/*observations: */
@@ -220,6 +221,6 @@
 	void	CreateElementsNodesAndMaterialsMelting(DataSet** pelements,DataSet** pnodes, DataSet** pmaterials, Model* model,ConstDataHandle model_handle);
 	void	CreateConstraintsMelting(DataSet** pconstraints,Model* model,ConstDataHandle model_handle);
-	void    CreateLoadsMelting(DataSet** ploads, Model* model, ConstDataHandle model_handle);
-	void    CreateParametersMelting(DataSet** pparameters,Model* model,ConstDataHandle model_handle);
+	void  CreateLoadsMelting(DataSet** ploads, Model* model, ConstDataHandle model_handle);
+	void  CreateParametersMelting(DataSet** pparameters,Model* model,ConstDataHandle model_handle);
 
 
Index: /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp	(revision 574)
+++ /issm/trunk/src/c/ModelProcessorx/Thermal/CreateLoadsThermal.cpp	(revision 575)
@@ -85,4 +85,2 @@
 
 }
-
-