Index: /issm/trunk/src/c/objects/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Pengrid.cpp	(revision 926)
+++ /issm/trunk/src/c/objects/Pengrid.cpp	(revision 927)
@@ -349,9 +349,4 @@
 	if (temperature<t_pmp){ //If T<Tpmp, there must be no melting. Therefore, melting should be  constrained to 0 when T<Tpmp, instead of using spcs, use penalties
 		Ke[0][0]=kmax*pow(10,penalty_offset);
-		if (id==65){
-			printf("Ke_terms[0] = %10.10lf\n",Ke[0]);
-			printf("penalty_offset(K) = %g\n",penalty_offset);
-			printf("kmax = %10.10lf\n",kmax);
-		}
 	}
 	
@@ -507,9 +502,4 @@
 		if (sub_analysis_type==SteadyAnalysisEnum()){
 			P_terms[0]=melting_offset*pow(10,penalty_offset)*(temperature-t_pmp);
-			if (id==65){
-				printf("penalty_offset(P) = %g\n",penalty_offset);
-				printf("melting_offset(P) = %10.10lf\n",melting_offset);
-				printf("P_terms[0] = %10.10lf\n",P_terms[0]);
-			}
 		}
 		else{