Index: /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 570)
+++ /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 571)
@@ -51,4 +51,5 @@
 	Reducematrixfromgtofx(&Kff,&Kfs,Kgg,fem->Gmn,fem->nodesets);
 
+
 	/*!Reduce load from g to f size: */
 	if (debug) _printf_("   reducing load from g to f set\n");
Index: /issm/trunk/src/c/parallel/thermal.cpp
===================================================================
--- /issm/trunk/src/c/parallel/thermal.cpp	(revision 570)
+++ /issm/trunk/src/c/parallel/thermal.cpp	(revision 571)
@@ -115,8 +115,4 @@
 		if(debug)_printf_("computing melting:\n");
 		diagnostic_core_linear(&m_g[0],&femmodels[1],inputs,MeltingAnalysisEnum(),SteadyAnalysisEnum());
-		PetscSynchronizedPrintf(MPI_COMM_WORLD,"melting_offset=%g\n",melting_offset);
-		PetscSynchronizedFlush(MPI_COMM_WORLD);
-		VecView(m_g[0],PETSC_VIEWER_STDOUT_WORLD);
-
 	}
 	else{